Showing 1–33 of 33 results for author: Marsman, M

How to verify the precision of density-functional-theory implementations via reproducible and universal workflows

Authors: Emanuele Bosoni, Louis Beal, Marnik Bercx, Peter Blaha, Stefan Blügel, Jens Bröder, Martin Callsen, Stefaan Cottenier, Augustin Degomme, Vladimir Dikan, Kristjan Eimre, Espen Flage-Larsen, Marco Fornari, Alberto Garcia, Luigi Genovese, Matteo Giantomassi, Sebastiaan P. Huber, Henning Janssen, Georg Kastlunger, Matthias Krack, Georg Kresse, Thomas D. Kühne, Kurt Lejaeghere, Georg K. H. Madsen, Martijn Marsman , et al. (20 additional authors not shown)

Abstract: In the past decades many density-functional theory methods and codes adopting periodic boundary conditions have been developed and are now extensively used in condensed matter physics and materials science research. Only in 2016, however, their precision (i.e., to which extent properties computed with different codes agree among each other) was systematically assessed on elemental crystals: a firs… ▽ More In the past decades many density-functional theory methods and codes adopting periodic boundary conditions have been developed and are now extensively used in condensed matter physics and materials science research. Only in 2016, however, their precision (i.e., to which extent properties computed with different codes agree among each other) was systematically assessed on elemental crystals: a first crucial step to evaluate the reliability of such computations. We discuss here general recommendations for verification studies aiming at further testing precision and transferability of density-functional-theory computational approaches and codes. We illustrate such recommendations using a greatly expanded protocol covering the whole periodic table from Z=1 to 96 and characterizing 10 prototypical cubic compounds for each element: 4 unaries and 6 oxides, spanning a wide range of coordination numbers and oxidation states. The primary outcome is a reference dataset of 960 equations of state cross-checked between two all-electron codes, then used to verify and improve nine pseudopotential-based approaches. Such effort is facilitated by deploying AiiDA common workflows that perform automatic input parameter selection, provide identical input/output interfaces across codes, and ensure full reproducibility. Finally, we discuss the extent to which the current results for total energies can be reused for different goals (e.g., obtaining formation energies). △ Less

Submitted 26 May, 2023; originally announced May 2023.

Comments: Main text: 23 pages, 4 figures. Supplementary: 68 pages. Nature Review Physics 2023

Journal ref: Nat. Rev. Phys. 6, 45 (2024)

Zero-point Renormalization of the Band Gap of Semiconductors and Insulators Using the PAW Method

Authors: Manuel Engel, Henrique Miranda, Laurent Chaput, Atsushi Togo, Carla Verdi, Martijn Marsman, Georg Kresse

Abstract: We evaluate the zero-point renormalization (ZPR) due to electron-phonon interactions of 28 solids using the projector-augmented-wave (PAW) method. The calculations cover diamond, many zincblende semiconductors, rock-salt and wurtzite oxides, as well as silicate and titania. Particular care is taken to include long-range electrostatic interactions via a generalized Fröhlich model, as discussed in P… ▽ More We evaluate the zero-point renormalization (ZPR) due to electron-phonon interactions of 28 solids using the projector-augmented-wave (PAW) method. The calculations cover diamond, many zincblende semiconductors, rock-salt and wurtzite oxides, as well as silicate and titania. Particular care is taken to include long-range electrostatic interactions via a generalized Fröhlich model, as discussed in Phys. Rev. Lett. 115, 176401 (2015) and Phys. Rev. B 92, 054307 (2015). The data are compared to recent calculations, npj Computational Materials 6, 167 (2020), and generally very good agreement is found. We discuss in detail the evaluation of the electron-phonon matrix elements within the PAW method. We show that two distinct versions can be obtained depending on when the atomic derivatives are taken. If the PAW transformation is applied before taking derivatives with respect to the ionic positions, equations similar to the ones conventionally used in pseudopotential codes are obtained. If the PAW transformation is used after taking the derivatives, the full-potential spirit is largely maintained. We show that both variants yield very similar ZPRs for selected materials when the rigid-ion approximation is employed. In practice, we find however that the pseudo version converges more rapidly with respect to the number of included unoccupied states. △ Less

Submitted 28 June, 2022; v1 submitted 9 May, 2022; originally announced May 2022.

Comments: Added pseudopotential information; revised Sec. II B to avoid possible misconceptions; revised discussion in Sec. III A; added diamond calculation with harder pseudopotential; corrected ZPR comparisons for AlSb and GaP due to inconsistent band-gap transitions; conclusion unchanged

Relations between networks, regression, partial correlation, and latent variable model

Authors: Lourens Waldorp, Maarten Marsman

Abstract: The Gaussian graphical model (GGM) has become a popular tool for analyzing networks of psychological variables. In a recent paper in this journal, Forbes, Wright, Markon, and Krueger (FWMK) voiced the concern that GGMs that are estimated from partial correlations wrongfully remove the variance that is shared by its constituents. If true, this concern has grave consequences for the application of G… ▽ More The Gaussian graphical model (GGM) has become a popular tool for analyzing networks of psychological variables. In a recent paper in this journal, Forbes, Wright, Markon, and Krueger (FWMK) voiced the concern that GGMs that are estimated from partial correlations wrongfully remove the variance that is shared by its constituents. If true, this concern has grave consequences for the application of GGMs. Indeed, if partial correlations only capture the unique covariances, then the data that come from a unidimensional latent variable model ULVM should be associated with an empty network (no edges), as there are no unique covariances in a ULVM. We know that this cannot be true, which suggests that FWMK are missing something with their claim. We introduce a connection between the ULVM and the GGM and use that connection to prove that we find a fully-connected and not an empty network associated with a ULVM. We then use the relation between GGMs and linear regression to show that the partial correlation indeed does not remove the common variance. △ Less

Submitted 21 July, 2020; originally announced July 2020.

Comments: 21 pages 2 figures

MSC Class: 62P15 Primary)

Reliability of decisions based on tests: Fourier analysis of Boolean decision functions

Authors: Lourens Waldorp, Maarten Marsman, Denny Borsboom

Abstract: Items in a test are often used as a basis for making decisions and such tests are therefore required to have good psychometric properties, like unidimensionality. In many cases the sum score is used in combination with a threshold to decide between pass or fail, for instance. Here we consider whether such a decision function is appropriate, without a latent variable model, and which properties of… ▽ More Items in a test are often used as a basis for making decisions and such tests are therefore required to have good psychometric properties, like unidimensionality. In many cases the sum score is used in combination with a threshold to decide between pass or fail, for instance. Here we consider whether such a decision function is appropriate, without a latent variable model, and which properties of a decision function are desirable. We consider reliability (stability) of the decision function, i.e., does the decision change upon perturbations, or changes in a fraction of the outcomes of the items (measurement error). We are concerned with questions of whether the sum score is the best way to aggregate the items, and if so why. We use ideas from test theory, social choice theory, graphical models, computer science and probability theory to answer these questions. We conclude that a weighted sum score has desirable properties that (i) fit with test theory and is observable (similar to a condition like conditional association), (ii) has the property that a decision is stable (reliable), and (iii) satisfies Rousseau's criterion that the input should match the decision. We use Fourier analysis of Boolean functions to investigate whether a decision function is stable and to figure out which (set of) items has proportionally too large an influence on the decision. To apply these techniques we invoke ideas from graphical models and use a pseudo-likelihood factorisation of the probability distribution. △ Less

Submitted 11 July, 2020; originally announced July 2020.

Comments: 41 pages, 4 figures

MSC Class: 62P15 (Primary) 94D10 (Secondary)

A Generalization of the Savage-Dickey Density Ratio for Testing Equality and Order Constrained Hypotheses

Authors: J. Mulder, E. -J. Wagenmakers, M. Marsman

Abstract: The Savage-Dickey density ratio is a specific expression of the Bayes factor when testing a precise (equality constrained) hypothesis against an unrestricted alternative. The expression greatly simplifies the computation of the Bayes factor at the cost of assuming a specific form of the prior under the precise hypothesis as a function of the unrestricted prior. A generalization was proposed by Ver… ▽ More The Savage-Dickey density ratio is a specific expression of the Bayes factor when testing a precise (equality constrained) hypothesis against an unrestricted alternative. The expression greatly simplifies the computation of the Bayes factor at the cost of assuming a specific form of the prior under the precise hypothesis as a function of the unrestricted prior. A generalization was proposed by Verdinelli and Wasserman (1995) such that the priors can be freely specified under both hypotheses while keeping the computational advantage. This paper presents an extension of this generalization when the hypothesis has equality as well as order constraints on the parameters of interest. The methodology is used for a constrained multivariate t test using the JZS Bayes factor and a constrained hypothesis test under the multinomial model. △ Less

Submitted 14 July, 2020; v1 submitted 21 April, 2020; originally announced April 2020.

Comments: 23 pages, 1 figure

Electron-Phonon Interactions Using the PAW Method and Wannier Functions

Authors: Manuel Engel, Martijn Marsman, Cesare Franchini, Georg Kresse

Abstract: We present an ab-initio density-functional-theory approach for calculating electron-phonon interactions within the projector augmented-wave method. The required electron-phonon matrix elements are defined as the second derivative of the one-electron energies in the PAW method. As the PAW method leads to a generalized eigenvalue problem, the resulting electron-phonon matrix elements lack some symme… ▽ More We present an ab-initio density-functional-theory approach for calculating electron-phonon interactions within the projector augmented-wave method. The required electron-phonon matrix elements are defined as the second derivative of the one-electron energies in the PAW method. As the PAW method leads to a generalized eigenvalue problem, the resulting electron-phonon matrix elements lack some symmetries that are usually present for simple eigenvalue problems and all-electron formulations. We discuss the relation between our definition of the electron-phonon matrix element and other formulations. To allow for efficient evaluation of physical properties, we introduce a Wannier-interpolation scheme, again adapted to generalized eigenvalue problems. To explore the method's numerical characteristics, the temperature-dependent band-gap renormalization of diamond is calculated and compared with previous publications. Furthermore, we apply the method to selected binary compounds and show that the obtained zero-point renormalizations agree well with other values found in literature and experiments. △ Less

Submitted 23 January, 2020; originally announced January 2020.

Journal ref: Phys. Rev. B 101, 184302 (2020)

Intervention in undirected Ising graphs and the partition function

Authors: Lourens Waldorp, Maarten Marsman

Abstract: Undirected graphical models have many applications in such areas as machine learning, image processing, and, recently, psychology. Psychopathology in particular has received a lot of attention, where symptoms of disorders are assumed to influence each other. One of the most relevant questions practically is on which symptom (node) to intervene to have the most impact. Interventions in undirected g… ▽ More Undirected graphical models have many applications in such areas as machine learning, image processing, and, recently, psychology. Psychopathology in particular has received a lot of attention, where symptoms of disorders are assumed to influence each other. One of the most relevant questions practically is on which symptom (node) to intervene to have the most impact. Interventions in undirected graphical models is equal to conditioning, and so we have available the machinery with the Ising model to determine the best strategy to intervene. In order to perform such calculations the partition function is required, which is computationally difficult. Here we use a Curie-Weiss approach to approximate the partition function in applications of interventions. We show that when the connection weights in the graph are equal within each clique then we obtain exactly the correct partition function. And if the weights vary according to a sub-Gaussian distribution, then the approximation is exponentially close to the correct one. We confirm these results with simulations. △ Less

Submitted 27 May, 2019; originally announced May 2019.

Comments: Preprint for original paper

The physics of (ir)rational choice

Authors: Joost Kruis, Gunter Maris, Maarten Marsman, Dylan Molenaar, Maria Bolsinova, Han L. J. van der Maas

Abstract: Even though classic theories and models of discrete choice pose man as a rational being, it has been shown extensively that people persistently violate rationality in their actual choices. Recent models of decision-making take these violations often (partially) into account, however, a unified framework has not been established. Here we propose such a framework, inspired by the Ising model from st… ▽ More Even though classic theories and models of discrete choice pose man as a rational being, it has been shown extensively that people persistently violate rationality in their actual choices. Recent models of decision-making take these violations often (partially) into account, however, a unified framework has not been established. Here we propose such a framework, inspired by the Ising model from statistical physics, and show that representing choice problems as a graph, together with a simple choice process, allows us to explain both rational decisions as well as violations of rationality. △ Less

Submitted 18 April, 2019; originally announced April 2019.

Assessing model-dielectric-dependent hybrid functionals on the antiferromagnetic transition-metal monoxides MnO, FeO, CoO, and NiO

Authors: Peitao Liu, Cesare Franchini, Martijn Marsman, Georg Kresse

Abstract: Recently, two nonempirical hybrid functionals, dielectric-dependent range-separated hybrid functional based on the Coulomb-attenuating method (DD-RSH-CAM) and doubly screened hybrid functional (DSH), have been suggested by [Chen et al, Phys. Rev. Mater. 2, 073803 (2018)] and [Cui et al, J. Phys. Chem. Lett. 9, 2338 (2018)], respectively. These two hybrid functionals are both based on a common mode… ▽ More Recently, two nonempirical hybrid functionals, dielectric-dependent range-separated hybrid functional based on the Coulomb-attenuating method (DD-RSH-CAM) and doubly screened hybrid functional (DSH), have been suggested by [Chen et al, Phys. Rev. Mater. 2, 073803 (2018)] and [Cui et al, J. Phys. Chem. Lett. 9, 2338 (2018)], respectively. These two hybrid functionals are both based on a common model dielectric function approach, but differ in the way how to non-empirically obtain the range-separation parameter. By retaining the full short-range Fock exchange and a fraction of the long-range Fock exchange that equals the inverse of the dielectric constant, both DD-RSH-CAM and DSH turn out to perform very well in predicting the band gaps for a large variety of semiconductors and insulators. Here, we assess how these two hybrid functionals perform on challenging antiferromagnetic transition-metal monoxides MnO, FeO, CoO, and NiO by comparing them to other conventional hybrid functionals and the $GW$ method. We find that single-shot DD0-RSH-CAM and DSH0 improve the band gaps towards experiments as compared to conventional hybrid functionals. The magnetic moments are slightly increased, but the predicted dielectric constants are decreased. The valence band density of states (DOS) predicted by DD0-RSH-CAM and DSH0 are as satisfactory as HSE03 in comparison to experimental spectra, however, the conduction band DOS are shifted to higher energies by about 2 eV compared to HSE03. Self-consistent DD-RSH-CAM and DSH deteriorate the results with a significant overestimation of band gaps. △ Less

Submitted 20 August, 2019; v1 submitted 13 February, 2019; originally announced February 2019.

Comments: 11 pages, 5 figures

Journal ref: J. Phys.: Condens. Matter 32 (2019) 015502

Interpreting the Ising Model: The Input Matters

Authors: Jonas Haslbeck, Sacha Epskamp, Maarten Marsman, Lourens Waldorp

Abstract: The Ising model is a model for pairwise interactions between binary variables that has become popular in the psychological sciences. It has been first introduced as a theoretical model for the alignment between positive (+1) and negative (-1) atom spins. In many psychological applications, however, the Ising model is defined on the domain $\{0,1\}$ instead of the classical domain $\{-1,1\}$. While… ▽ More The Ising model is a model for pairwise interactions between binary variables that has become popular in the psychological sciences. It has been first introduced as a theoretical model for the alignment between positive (+1) and negative (-1) atom spins. In many psychological applications, however, the Ising model is defined on the domain $\{0,1\}$ instead of the classical domain $\{-1,1\}$. While it is possible to transform the parameters of a given Ising model in one domain to obtain a statistically equivalent model in the other domain, the parameters in the two versions of the Ising model lend themselves to different interpretations and imply different dynamics, when studying the Ising model as a dynamical system. In this tutorial paper, we provide an accessible discussion of the interpretation of threshold and interaction parameters in the two domains and show how the dynamics of the Ising model depends on the choice of domain. Finally, we provide a transformation that allows to transform the parameters in an Ising model in one domain into a statistically equivalent Ising model in the other domain. △ Less

Submitted 13 March, 2020; v1 submitted 7 November, 2018; originally announced November 2018.

Combine Statistical Thinking With Open Scientific Practice: A Protocol of a Bayesian Research Project

Authors: Alexandra Sarafoglou, Anna van der Heijden, Tim Draws, Joran Cornelisse, Eric-Jan Wagenmakers, Maarten Marsman

Abstract: Current developments in the statistics community suggest that modern statistics education should be structured holistically, that is, by allowing students to work with real data and to answer concrete statistical questions, but also by educating them about alternative frameworks, such as Bayesian inference. In this article, we describe how we incorporated such a holistic structure in a Bayesian re… ▽ More Current developments in the statistics community suggest that modern statistics education should be structured holistically, that is, by allowing students to work with real data and to answer concrete statistical questions, but also by educating them about alternative frameworks, such as Bayesian inference. In this article, we describe how we incorporated such a holistic structure in a Bayesian research project on ordered binomial probabilities. The project was conducted with a group of three undergraduate psychology students who had basic knowledge of Bayesian statistics and programming, but lacked formal mathematical training. The research project aimed to (1) convey the basic mathematical concepts of Bayesian inference; (2) have students experience the entire empirical cycle including collection, analysis, and interpretation of data and (3) teach students open science practices. △ Less

Submitted 6 January, 2022; v1 submitted 17 October, 2018; originally announced October 2018.

An Asymptotically Efficient Metropolis-Hastings Sampler for Bayesian Inference in Large-Scale Educational Measuremen

Authors: Timo Bechger, Gunter Maris, Maarten Marsman

Abstract: This paper discusses a Metropolis-Hastings algorithm developed by \citeA{MarsmanIsing}. The algorithm is derived from first principles, and it is proven that the algorithm becomes more efficient with more data and meets the growing demands of large scale educational measurement. This paper discusses a Metropolis-Hastings algorithm developed by \citeA{MarsmanIsing}. The algorithm is derived from first principles, and it is proven that the algorithm becomes more efficient with more data and meets the growing demands of large scale educational measurement. △ Less

Submitted 12 August, 2018; originally announced August 2018.

Logistic regression and Ising networks: prediction and estimation when violating lasso assumptions

Authors: Lourens Waldorp, Maarten Marsman, Gunter Maris

Abstract: The Ising model was originally developed to model magnetisation of solids in statistical physics. As a network of binary variables with the probability of becoming 'active' depending only on direct neighbours, the Ising model appears appropriate for many other processes. For instance, it was recently applied in psychology to model co-occurrences of mental disorders. It has been shown that the conn… ▽ More The Ising model was originally developed to model magnetisation of solids in statistical physics. As a network of binary variables with the probability of becoming 'active' depending only on direct neighbours, the Ising model appears appropriate for many other processes. For instance, it was recently applied in psychology to model co-occurrences of mental disorders. It has been shown that the connections between the variables (nodes) in the Ising network can be estimated with a series of logistic regressions. This naturally leads to questions of how well such a model predicts new observations and how well parameters of the Ising model can be estimated using logistic regressions. Here we focus on the high-dimensional setting with more parameters than observations and consider violations of assumptions of the lasso. In particular, we determine the consequences for both prediction and estimation when the sparsity and restricted eigenvalue assumptions are not satisfied. We explain by using the idea of connected copies (extreme multicollinearity) the fact that prediction becomes better when either sparsity or multicollinearity is not satisfied. We illustrate these results with simulations. △ Less

Submitted 28 July, 2018; originally announced July 2018.

Comments: to appear, Behaviormetrika, 2018

Charge self-consistent many-body corrections using optimized projected localized orbitals

Authors: Malte Schüler, Oleg E. Peil, Gernot J. Kraberger, Ronald Pordzik, Martijn Marsman, Georg Kresse, Tim O. Wehling, Markus Aichhorn

Abstract: In order for methods combining ab initio density-functional theory and many-body techniques to become routinely used, a flexible, fast, and easy-to-use implementation is crucial. We present an implementation of a general charge self-consistent scheme based on projected localized orbitals in the projector augmented wave framework in the Vienna Ab Initio Simulation Package (VASP). We give a detailed… ▽ More In order for methods combining ab initio density-functional theory and many-body techniques to become routinely used, a flexible, fast, and easy-to-use implementation is crucial. We present an implementation of a general charge self-consistent scheme based on projected localized orbitals in the projector augmented wave framework in the Vienna Ab Initio Simulation Package (VASP). We give a detailed description on how the projectors are optimally chosen and how the total energy is calculated. We benchmark our implementation in combination with dynamical mean-field theory: first we study the charge-transfer insulator NiO using a Hartree-Fock approach to solve the many-body Hamiltonian. We address the advantages of the optimized against non-optimized projectors and furthermore find that charge self-consistency decreases the dependence of the spectral function - especially the gap - on the double counting. Second, using continuous-time quantum Monte Carlo we study a monolayer of SrVO$_3$, where strong orbital polarization occurs due to the reduced dimensionality. Using total-energy calculation for structure determination, we find that electronic correlations have a non-negligible influence on the position of the apical oxygens, and therefore on the thickness of the single SrVO$_3$ layer. △ Less

Submitted 5 November, 2018; v1 submitted 5 April, 2018; originally announced April 2018.

Comments: 11 pages, 6 figures

Journal ref: J. Phys.: Condens. Matter 30, 475901 (2018)

Bayesian Rank-Based Hypothesis Testing for the Rank Sum Test, the Signed Rank Test, and Spearman's $ρ$

Authors: Johnny van Doorn, Alexander Ly, Maarten Marsman, Eric-Jan Wagenmakers

Abstract: Bayesian inference for rank-order problems is frustrated by the absence of an explicit likelihood function. This hurdle can be overcome by assuming a latent normal representation that is consistent with the ordinal information in the data: the observed ranks are conceptualized as an impoverished reflection of an underlying continuous scale, and inference concerns the parameters that govern the lat… ▽ More Bayesian inference for rank-order problems is frustrated by the absence of an explicit likelihood function. This hurdle can be overcome by assuming a latent normal representation that is consistent with the ordinal information in the data: the observed ranks are conceptualized as an impoverished reflection of an underlying continuous scale, and inference concerns the parameters that govern the latent representation. We apply this generic data-augmentation method to obtain Bayes factors for three popular rank-based tests: the rank sum test, the signed rank test, and Spearman's $ρ_s$. △ Less

Submitted 17 May, 2019; v1 submitted 19 December, 2017; originally announced December 2017.

Comments: 33 pages, 12 figures

A Tutorial on Fisher Information

Authors: Alexander Ly, Maarten Marsman, Josine Verhagen, Raoul Grasman, Eric-Jan Wagenmakers

Abstract: In many statistical applications that concern mathematical psychologists, the concept of Fisher information plays an important role. In this tutorial we clarify the concept of Fisher information as it manifests itself across three different statistical paradigms. First, in the frequentist paradigm, Fisher information is used to construct hypothesis tests and confidence intervals using maximum like… ▽ More In many statistical applications that concern mathematical psychologists, the concept of Fisher information plays an important role. In this tutorial we clarify the concept of Fisher information as it manifests itself across three different statistical paradigms. First, in the frequentist paradigm, Fisher information is used to construct hypothesis tests and confidence intervals using maximum likelihood estimators; second, in the Bayesian paradigm, Fisher information is used to define a default prior; lastly, in the minimum description length paradigm, Fisher information is used to measure model complexity. △ Less

Submitted 17 October, 2017; v1 submitted 2 May, 2017; originally announced May 2017.

MSC Class: 62-01; 62B10 (Primary); 62F03; 62F12; 62F15; 62B10 (Secondary)

Embedding for bulk systems using localized atomic orbitals

Authors: F. Libisch, M. Marsman, J. Burgdörfer, G. Kresse

Abstract: We present an embedding approach for semiconductors and insulators based on or- bital rotations in the space of occupied Kohn-Sham orbitals. We have implemented our approach in the popular VASP software package. We demonstrate its power for defect structures in silicon and polaron formation in titania, two challenging cases for conventional Kohn-Sham density functional theory. We present an embedding approach for semiconductors and insulators based on or- bital rotations in the space of occupied Kohn-Sham orbitals. We have implemented our approach in the popular VASP software package. We demonstrate its power for defect structures in silicon and polaron formation in titania, two challenging cases for conventional Kohn-Sham density functional theory. △ Less

Submitted 29 March, 2017; originally announced March 2017.

Comments: 2 figures

Journal ref: J. Chem. Phys. 147, 034110 (2017)

A Tutorial on Bridge Sampling

Authors: Quentin F. Gronau, Alexandra Sarafoglou, Dora Matzke, Alexander Ly, Udo Boehm, Maarten Marsman, David S. Leslie, Jonathan J. Forster, Eric-Jan Wagenmakers, Helen Steingroever

Abstract: The marginal likelihood plays an important role in many areas of Bayesian statistics such as parameter estimation, model comparison, and model averaging. In most applications, however, the marginal likelihood is not analytically tractable and must be approximated using numerical methods. Here we provide a tutorial on bridge sampling (Bennett, 1976; Meng & Wong, 1996), a reliable and relatively str… ▽ More The marginal likelihood plays an important role in many areas of Bayesian statistics such as parameter estimation, model comparison, and model averaging. In most applications, however, the marginal likelihood is not analytically tractable and must be approximated using numerical methods. Here we provide a tutorial on bridge sampling (Bennett, 1976; Meng & Wong, 1996), a reliable and relatively straightforward sampling method that allows researchers to obtain the marginal likelihood for models of varying complexity. First, we introduce bridge sampling and three related sampling methods using the beta-binomial model as a running example. We then apply bridge sampling to estimate the marginal likelihood for the Expectancy Valence (EV) model---a popular model for reinforcement learning. Our results indicate that bridge sampling provides accurate estimates for both a single participant and a hierarchical version of the EV model. We conclude that bridge sampling is an attractive method for mathematical psychologists who typically aim to approximate the marginal likelihood for a limited set of possibly high-dimensional models. △ Less

Submitted 11 October, 2017; v1 submitted 17 March, 2017; originally announced March 2017.

Bayesian Estimation of Kendall's tau Using a Latent Normal Approach

Authors: Johnny van Doorn, Alexander Ly, Maarten Marsman, Eric-Jan Wagenmakers

Abstract: The rank-based association between two variables can be modeled by introducing a latent normal level to ordinal data. We demonstrate how this approach yields Bayesian inference for Kendall's rank correlation coefficient, improving on a recent Bayesian solution from asymptotic properties of the test statistic. The rank-based association between two variables can be modeled by introducing a latent normal level to ordinal data. We demonstrate how this approach yields Bayesian inference for Kendall's rank correlation coefficient, improving on a recent Bayesian solution from asymptotic properties of the test statistic. △ Less

Submitted 24 May, 2018; v1 submitted 6 March, 2017; originally announced March 2017.

Comments: 13 pages, 2 figures

NMR shieldings from density functional perturbation theory: GIPAW versus all-electron calculations

Authors: G. A. de Wijs, R. Laskowski, P. Blaha, R. W. A. Havenith, G. Kresse, M. Marsman

Abstract: We present a benchmark of the density functional linear response calculation of NMR shieldings within the Gauge-Including Projector-Augmented-Wave method against all-electron Augmented-Plane-Wave$+$local-orbital and uncontracted Gaussian basis set results for NMR shieldings in molecular and solid state systems. In general, excellent agreement between the aforementioned methods is obtained. Scalar… ▽ More We present a benchmark of the density functional linear response calculation of NMR shieldings within the Gauge-Including Projector-Augmented-Wave method against all-electron Augmented-Plane-Wave$+$local-orbital and uncontracted Gaussian basis set results for NMR shieldings in molecular and solid state systems. In general, excellent agreement between the aforementioned methods is obtained. Scalar relativistic effects are shown to be quite large for nuclei in molecules in the deshielded limit. The small component makes up a substantial part of the relativistic corrections. △ Less

Submitted 13 September, 2016; originally announced September 2016.

Comments: 3 figures, supplementary material included

Estimating psychopathological networks: be careful what you wish for

Authors: Sacha Epskamp, Joost Kruis, Maarten Marsman

Abstract: Network models, in which psychopathological disorders are conceptualized as a complex interplay of psychological and biological components, have become increasingly popular in the recent psychopathological literature. These network models often contain significant numbers of unknown parameters, yet the sample sizes available in psychological research are limited. As such, general assumptions about… ▽ More Network models, in which psychopathological disorders are conceptualized as a complex interplay of psychological and biological components, have become increasingly popular in the recent psychopathological literature. These network models often contain significant numbers of unknown parameters, yet the sample sizes available in psychological research are limited. As such, general assumptions about the true network are introduced to reduce the number of free parameters. Incorporating these assumptions, however, means that the resulting network will lead to reflect the particular structure assumed by the estimation method---a crucial and often ignored aspect of psychopathological networks. For example, observing a sparse structure and simultaneously assuming a sparse structure does not imply that the true model is, in fact, sparse. To illustrate this point, we discuss recent literature and show the effect of the assumption of sparsity in three simulation studies. △ Less

Submitted 11 September, 2017; v1 submitted 26 April, 2016; originally announced April 2016.

Comments: Published in PlosOne

Analytic Posteriors for Pearson's Correlation Coefficient

Authors: Alexander Ly, Maarten Marsman, Eric-Jan Wagenmakers

Abstract: Pearson's correlation is one of the most common measures of linear dependence. Recently, Bernardo (2015) introduced a flexible class of priors to study this measure in a Bayesian setting. For this large class of priors we show that the (marginal) posterior for Pearson's correlation coefficient and all of the posterior moments are analytic. Our results are available in the open-source software pack… ▽ More Pearson's correlation is one of the most common measures of linear dependence. Recently, Bernardo (2015) introduced a flexible class of priors to study this measure in a Bayesian setting. For this large class of priors we show that the (marginal) posterior for Pearson's correlation coefficient and all of the posterior moments are analytic. Our results are available in the open-source software package JASP. △ Less

Submitted 28 April, 2017; v1 submitted 5 October, 2015; originally announced October 2015.

MSC Class: 62H20; 62E15; 62F15

Anisotropic Magnetic Couplings and Structure-Driven Canted to Collinear Transitions in Spin-orbit Coupled Sr2IrO4

Authors: Peitao Liu, Sergii Khmelevskyi, Bongjae Kim, Martijn Marsman, Dianzhong Li, Xing-Qiu Chen, D. D. Sarma, Georg Kresse, Cesare Franchini

Abstract: We put forward a scheme to study the anisotropic magnetic couplings in Sr2IrO4 by mapping fully relativistic constrained noncollinear density functional theory including an on-site Hubbard U correction onto a general spin model Hamiltonian. This procedure allows for the simultaneous account and direct control of the lattice, spin and orbital interactions within a fully ab initio scheme. We compute… ▽ More We put forward a scheme to study the anisotropic magnetic couplings in Sr2IrO4 by mapping fully relativistic constrained noncollinear density functional theory including an on-site Hubbard U correction onto a general spin model Hamiltonian. This procedure allows for the simultaneous account and direct control of the lattice, spin and orbital interactions within a fully ab initio scheme. We compute the isotropic, single site anisotropy and Dzyaloshinskii-Moriya (DM) coupling parameters, and clarify that the origin of the canted magnetic state in Sr2IrO4 arises from the interplay between structural distortions and the competition between isotropic exchange and DM interactions. A complete magnetic phase diagram with respect to the tetragonal distortion and the rotation of IrO6 octahedra is constructed, revealing the presence of two types of canted to collinear magnetic transitions: a spin-flop transition with increasing tetragonal distortion and a complete quenching of the basal weak ferromagnetic moment below a critical octahedral rotation. △ Less

Submitted 24 March, 2015; v1 submitted 23 March, 2015; originally announced March 2015.

Journal ref: Phys. Rev. B 92, 054428 (2015)

The formation of the positive, fixed charge at c-Si(111)/a-Si$_3$N$_{3.5}$:H interfaces

Authors: L. E. Hintzsche, C. M. Fang, M. Marsman, M. W. P. E. Lamers, A. W. Weeber, G. Kresse

Abstract: Modern electronic devices are unthinkable without the well-controlled formation of interfaces at heterostructures. These often involve at least one amorphous material. Modeling such interfaces poses a significant challenge, since a meaningful result can only be expected by using huge models or by drawing from many statistically independent samples. Here we report on the results of high throughput… ▽ More Modern electronic devices are unthinkable without the well-controlled formation of interfaces at heterostructures. These often involve at least one amorphous material. Modeling such interfaces poses a significant challenge, since a meaningful result can only be expected by using huge models or by drawing from many statistically independent samples. Here we report on the results of high throughput calculations for interfaces between crystalline silicon (c-Si) and amorphous silicon nitride (a-Si$_3$N$_{3.5}$:H), which are omnipresent in commercially available solar cells. The findings reconcile only partly understood key features. At the interface, threefold coordinated Si atoms are present. These are caused by the structural mismatch between the amorphous and crystalline part. The local Fermi level of undoped c-Si lies well below that of a-SiN:H. To align the Fermi levels in the device, charge is transferred from the a-SiN:H part to the c-Si part resulting in an abundance of positively charged, threefold coordinated Si atoms at the interface. This explains the existence of a positive, fixed charge at the interface that repels holes. △ Less

Submitted 15 January, 2015; originally announced January 2015.

Journal ref: Phys. Rev. Applied 3, 064005 (2015)

Electronic structure and energy level schemes of RE3+:LaSi3N5 and RE2+:LaSi3N5-xOx phosphors (RE= Ce, Pr, Nd, Pm, Sm, Eu) from first principles

Authors: Ismail A. M. Ibrahim, Zoltan Lences, Pavol Sajgalik, Lubomir Benco, Martijn Marsman

Abstract: First principles calculations of rare earth (RE)-doped LaSi3N5 host lattice are performed to obtain the electronic structure, the band gap (BG), and the character of electronic transitions. Doping with both trivalent and bivalent RE cations is inspected. RE 4f states form two bands of occupied and unoccupied states separated by ~5 eV. In RE3+-doped compounds 4f states are shifted by ~6 eV to more… ▽ More First principles calculations of rare earth (RE)-doped LaSi3N5 host lattice are performed to obtain the electronic structure, the band gap (BG), and the character of electronic transitions. Doping with both trivalent and bivalent RE cations is inspected. RE 4f states form two bands of occupied and unoccupied states separated by ~5 eV. In RE3+-doped compounds 4f states are shifted by ~6 eV to more negative energies compared with RE2+-compounds. This stabilization causes that RE3+ 4f bands are in a different position relative to the valence band and the conduction band than RE2+ 4f bands and therefore different electronic transitions apply. BG of RE3+-compounds decreases from ~4.6 eV (Ce) to ~0.5 eV (Eu). Except for Ce3+, exhibiting the 4f-5d transition, other RE3+-compounds show the charge transfer of the p - 4f character. BG of RE2+-compounds increases from ~0.80 eV (Ce, Pr) to ~0.95 eV (Nd, Pm), ~1.43 eV (Sm), and ~3.28 eV (Eu) and the electronic transition is of the 4f-5d character. The energy level scheme constructed from ab initio calculated electronic structures agrees well with the experimental energy level diagram. The agreement demonstrates the reliability of the hybrid functional HSE06 to describe correctly bands of nonbonding RE 4f electrons. △ Less

Submitted 20 October, 2014; originally announced October 2014.

Comments: 23 pages, 5 figures, no tables

MSC Class: 82D20 ACM Class: J.2

Magnetic structure map for face-centered tetragonal iron: appearance of a new collinear spin structure

Authors: D. Reith, R. Podloucky, M. Marsman, P. O. Bedolla-Velazquez, P. Mohn

Abstract: For face-centered cubic (fcc) and tetragonal (fct) iron a large number of magnetic configurations as a function of crystal structural parameters were studied by means of density functional theory. The stability of magnetic structures was defined by the magnetic re-orientation energy $Δ{E}^i_\text{reor}$ as the difference of the total energy of configuration $i$ and that of the fcc ferromagnetic st… ▽ More For face-centered cubic (fcc) and tetragonal (fct) iron a large number of magnetic configurations as a function of crystal structural parameters were studied by means of density functional theory. The stability of magnetic structures was defined by the magnetic re-orientation energy $Δ{E}^i_\text{reor}$ as the difference of the total energy of configuration $i$ and that of the fcc ferromagnetic state. The Cluster Expansion technique was applied to several volumes deriving $Δ{E}_\text{reor}$ for more than 90.000 collinear spin structures for each volume. A variety of structures with promisingly low $Δ{E}_\text{reor}$ were tetragonally distorted according to a two-dimensional mesh defined by volume per atom and $c/a$ ratio of the distortion. At each of the points on this mesh $ΔE_\text{reor}$ of all collinear structure were compared to results for non-collinear spin spirals (SS) which were calculated for a grid of propagation directions. The lowest $Δ{E}_\text{reor}$ of all investigated spin structures then defined the magnetic structure map spanned by volume per atom and $c/a$ ratio. Three local minima were identified and for each of the minima SS were calculated on a fine grid of propagation vectors. For the volume per atom of 10.6 Å$^3$ and the distortion range $0.94 \le c/a \le 1.01$ we found a new, surprisingly simple collinear spin structure with four atoms per fct unit cell to be the most stable one. This structure was called AFM/NM because it consists of two atoms with anti-ferromagnetically ordered local moments of $\pm 1.8 μ_\text{B}$ and of two nonmagnetic atoms with zero local moment. It seems that the newly detected AFM/NM structure explains a variety of puzzling experimental results. △ Less

Submitted 20 December, 2013; v1 submitted 16 December, 2013; originally announced December 2013.

Maximally localized Wannier functions in LaMnO3 within PBE+U, hybrid functionals, and partially self-consistent GW: an efficient route to construct ab-initio tight-binding parameters for e_g perovskites

Authors: C. Franchini, R. Kovacik, M. Marsman, S. Sathyanarayana Murthy, J. He, C. Ederer, G. Kresse

Abstract: Using the newly developed VASP2WANNIER90 interface we have constructed maximally localized Wannier functions (MLWFs) for the e_g states of the prototypical Jahn-Teller magnetic perovskite LaMnO3 at different levels of approximation for the exchange-correlation kernel. These include conventional density functional theory (DFT) with and without additional on-site Hubbard U term, hybrid-DFT, and part… ▽ More Using the newly developed VASP2WANNIER90 interface we have constructed maximally localized Wannier functions (MLWFs) for the e_g states of the prototypical Jahn-Teller magnetic perovskite LaMnO3 at different levels of approximation for the exchange-correlation kernel. These include conventional density functional theory (DFT) with and without additional on-site Hubbard U term, hybrid-DFT, and partially self-consistent GW. By suitably mapping the MLWFs onto an effective e_g tight-binding (TB) Hamiltonian we have computed a complete set of TB parameters which should serve as guidance for more elaborate treatments of correlation effects in effective Hamiltonian-based approaches. The method-dependent changes of the calculated TB parameters and their interplay with the electron-electron (el-el) interaction term are discussed and interpreted. We discuss two alternative model parameterizations: one in which the effects of the el-el interaction are implicitly incorporated in the otherwise "noninteracting" TB parameters, and a second where we include an explicit mean-field el-el interaction term in the TB Hamiltonian. Both models yield a set of tabulated TB parameters which provide the band dispersion in excellent agreement with the underlying ab initio and MLWF bands. △ Less

Submitted 14 May, 2012; v1 submitted 7 November, 2011; originally announced November 2011.

Comments: 30 pages, 7 figures

Journal ref: J. Phys.: Condens. Matter 24 (2012) 235602

Proximity of Iron Pnictide Superconductors to a Quantum Tricritical Point

Authors: G. Giovannetti, C. Ortix, M. Marsman, M. Capone, J. van den Brink, J. Lorenzana

Abstract: We determine the nature of the magnetic quantum critical point in the doped LaFeAsO using a set of constrained density functional calculations that provide ab initio coefficients for a Landau order parameter analysis. The system turns out to be remarkably close to a quantum tricritical point, where the nature of the phase transition changes from first to second order. We compare with the effective… ▽ More We determine the nature of the magnetic quantum critical point in the doped LaFeAsO using a set of constrained density functional calculations that provide ab initio coefficients for a Landau order parameter analysis. The system turns out to be remarkably close to a quantum tricritical point, where the nature of the phase transition changes from first to second order. We compare with the effective field theory and discuss the experimental consequences. △ Less

Submitted 31 August, 2010; originally announced September 2010.

Comments: 4 pages, 4 figures

Journal ref: Nat. Commun. 2, 298 (2011)

The multiferroic phase of DyFeO$_{3}$:an ab--initio study

Authors: Alessandro Stroppa, Martijn Marsman, Georg Kresse, Silvia Picozzi

Abstract: By performing accurate ab-initio density functional theory calculations, we study the role of $4f$ electrons in stabilizing the magnetic-field-induced ferroelectric state of DyFeO$_{3}$. We confirm that the ferroelectric polarization is driven by an exchange-strictive mechanism, working between adjacent spin-polarized Fe and Dy layers, as suggested by Y. Tokunaga [Phys. Rev. Lett, \textbf{101},… ▽ More By performing accurate ab-initio density functional theory calculations, we study the role of $4f$ electrons in stabilizing the magnetic-field-induced ferroelectric state of DyFeO$_{3}$. We confirm that the ferroelectric polarization is driven by an exchange-strictive mechanism, working between adjacent spin-polarized Fe and Dy layers, as suggested by Y. Tokunaga [Phys. Rev. Lett, \textbf{101}, 097205 (2008)]. A careful electronic structure analysis suggests that coupling between Dy and Fe spin sublattices is mediated by Dy-$d$ and O-$2p$ hybridization. Our results are robust with respect to the different computational schemes used for $d$ and $f$ localized states, such as the DFT+$U$ method, the Heyd-Scuseria-Ernzerhof (HSE) hybrid functional and the GW approach. Our findings indicate that the interaction between the $f$ and $d$ sublattice might be used to tailor ferroelectric and magnetic properties of multiferroic compounds. △ Less

Submitted 3 March, 2010; v1 submitted 17 December, 2009; originally announced December 2009.

Comments: 6 pages, 4 figures-Revised version

Determining the Anisotropic Exchange Coupling of CrO_2 via First-Principles Density Functional Theory Calculations

Authors: H. Sims, S. J. Oset, W. H. Butler, James M. MacLaren, Martijn Marsman

Abstract: We report a study of the anisotropic exchange interactions in bulk CrO_2 calculated from first principles within density functional theory. We determine the exchange coupling energies, using both the experimental lattice parameters and those obtained within DFT, within a modified Heisenberg model Hamiltonian in two ways. We employ a supercell method in which certain spins within a cell are rotated… ▽ More We report a study of the anisotropic exchange interactions in bulk CrO_2 calculated from first principles within density functional theory. We determine the exchange coupling energies, using both the experimental lattice parameters and those obtained within DFT, within a modified Heisenberg model Hamiltonian in two ways. We employ a supercell method in which certain spins within a cell are rotated and the energy dependence is calculated and a spin-spiral method that modifies the periodic boundary conditions of the problem to allow for an overall rotation of the spins between unit cells. Using the results from each of these methods, we calculate the spin-wave stiffness constant D from the exchange energies using the magnon dispersion relation. We employ a Monte Carlo method to determine the DFT-predicted Curie temperature from these calculated energies and compare with accepted values. Finally, we offer an evaluation of the accuracy of the DFT-based methods and suggest implications of the competing ferro- and antiferromagnetic interactions. △ Less

Submitted 30 January, 2013; v1 submitted 24 October, 2009; originally announced October 2009.

Comments: 10 pages, 13 figures

Structural, vibrational and quasiparticle properties of the Peierls semiconductor $\rm BaBiO_3$: a hybrid functional and self-consistent GW+vertex-corrections study

Authors: C. Franchini, A. Sanna, M. Marsman, G. Kresse

Abstract: $\rm BaBiO_3$ is characterized by a charge disproportionation with half of the Bi atoms possessing a valence 3+ and half a valence 5+. Because of selfinteraction errors, local and semi-local density functionals fail to describe the charge disproportionation quantitatively, yielding a too small structural distortion and no band gap. Using hybrid functionals we obtain a satisfactory description of… ▽ More $\rm BaBiO_3$ is characterized by a charge disproportionation with half of the Bi atoms possessing a valence 3+ and half a valence 5+. Because of selfinteraction errors, local and semi-local density functionals fail to describe the charge disproportionation quantitatively, yielding a too small structural distortion and no band gap. Using hybrid functionals we obtain a satisfactory description of the structural, electronic, optical, and vibrational properties of $\rm BaBiO_3$. The results obtained using GW (Green's function G and screened Coulomb potential W) based schemes on top of hybrid functionals, including fully selfconsistent GW calculations with vertex corrections in the dielectric screening, qualitatively confirm the HSE picture but a systematic overestimation of the bandgap by about 0.4 eV is observed. △ Less

Submitted 18 November, 2009; v1 submitted 5 March, 2008; originally announced March 2008.

Comments: 5 pages, 5 figures

Energies of ions in water and nanopores within Density Functional Theory

Authors: Kevin Leung, Martijn Marsman

Abstract: Accurate calculations of electrostatic potentials and treatment of substrate polarizability are critical for predicting the permeation of ions inside water-filled nanopores. The {\it ab initio} molecular dynamics method (AIMD), based on Density Functional Theory (DFT), accounts for the polarizability of materials, water, and solutes, and it should be the method of choice for predicting accurate… ▽ More Accurate calculations of electrostatic potentials and treatment of substrate polarizability are critical for predicting the permeation of ions inside water-filled nanopores. The {\it ab initio} molecular dynamics method (AIMD), based on Density Functional Theory (DFT), accounts for the polarizability of materials, water, and solutes, and it should be the method of choice for predicting accurate electrostatic energies of ions. In practice, DFT coupled with the use of periodic boundary conditions in a charged system leads to large energy shifts. Results obtained using different DFT packages may vary because of the way pseudopotentials and long-range electrostatics are implemented. Using maximally localized Wannier functions, we apply robust corrections that yield relatively unambiguous ion energies in select molecular and aqueous systems and inside carbon nanotubes. Large binding energies are predicted for ions in metallic carbon nanotube arrays, while with consistent definitions Na$^+$ and Cl$^-$ energies are found to exhibit asymmetries comparable with those computed using non-polarizable water force fields. △ Less

Submitted 2 November, 2007; originally announced November 2007.

Comments: 12 pages, 3 figures

Journal ref: Journal of Chemical Physics Vol. 127, No. 15, pages 154722 (19 October 2007)

A new polymorphic material? Structural degeneracy of ZrMn_2

Authors: Xing-Qiu Chen, W. Wolf, R. Podloucky, P. Rogl, M. Marsman

Abstract: Based on density functional calculations, we propose that ZrMn_2 is a polymorphic material. We predict that at low temperatures the cubic C15, and the hexagonal C14 and C36 structures of the Laves phase compound ZrMn_2 are nearly equally stable within 0.3 kJmol^{-1} or 30 K. This degeneracy occurs when the Mn atoms magnetize spontaneously in a ferromagnetic arrangement forming the states of lowe… ▽ More Based on density functional calculations, we propose that ZrMn_2 is a polymorphic material. We predict that at low temperatures the cubic C15, and the hexagonal C14 and C36 structures of the Laves phase compound ZrMn_2 are nearly equally stable within 0.3 kJmol^{-1} or 30 K. This degeneracy occurs when the Mn atoms magnetize spontaneously in a ferromagnetic arrangement forming the states of lowest energy. From the temperature dependent free energies at T approx 160K we predict a transition from the most stable C15 to the C14 structure, which is the experimentally observed structure at elevated temperatures. △ Less

Submitted 18 March, 2004; originally announced March 2004.

Comments: 4 pages, 3 figures

Journal ref: EuroPhys. Lett, 2004, Vol.67, No.5, pp.807-813