-
A Bayesian factor analysis model for high-dimensional microbiome count data
Authors:
Ismaïla Ba,
Maxime Turgeon,
Simona Veniamin,
Juan Joel,
Richard Miller,
Morag Graham,
Christine Bonner,
Charles N. Bernstein,
Douglas L. Arnold,
Amit Bar-Or,
Ruth Ann Marrie,
Julia O'Mahony,
E. Ann Yeh,
Brenda Banwell,
Emmanuelle Waubant,
Natalie Knox,
Gary Van Domselaar,
Ali I. Mirza,
Heather Armstrong,
Saman Muthukumarana,
Kevin McGregor
Abstract:
Dimension reduction techniques are among the most essential analytical tools in the analysis of high-dimensional data. Generalized principal component analysis (PCA) is an extension to standard PCA that has been widely used to identify low-dimensional features in high-dimensional discrete data, such as binary, multi-category and count data. For microbiome count data in particular, the multinomial…
▽ More
Dimension reduction techniques are among the most essential analytical tools in the analysis of high-dimensional data. Generalized principal component analysis (PCA) is an extension to standard PCA that has been widely used to identify low-dimensional features in high-dimensional discrete data, such as binary, multi-category and count data. For microbiome count data in particular, the multinomial PCA is a natural counterpart of the standard PCA. However, this technique fails to account for the excessive number of zero values, which is frequently observed in microbiome count data. To allow for sparsity, zero-inflated multivariate distributions can be used. We propose a zero-inflated probabilistic PCA model for latent factor analysis. The proposed model is a fully Bayesian factor analysis technique that is appropriate for microbiome count data analysis. In addition, we use the mean-field-type variational family to approximate the marginal likelihood and develop a classification variational approximation algorithm to fit the model. We demonstrate the efficiency of our procedure for predictions based on the latent factors and the model parameters through simulation experiments, showcasing its superiority over competing methods. This efficiency is further illustrated with two real microbiome count datasets. The method is implemented in R.
△ Less
Submitted 23 April, 2024; v1 submitted 3 April, 2024;
originally announced April 2024.
-
A foundation model for atomistic materials chemistry
Authors:
Ilyes Batatia,
Philipp Benner,
Yuan Chiang,
Alin M. Elena,
Dávid P. Kovács,
Janosh Riebesell,
Xavier R. Advincula,
Mark Asta,
Matthew Avaylon,
William J. Baldwin,
Fabian Berger,
Noam Bernstein,
Arghya Bhowmik,
Samuel M. Blau,
Vlad Cărare,
James P. Darby,
Sandip De,
Flaviano Della Pia,
Volker L. Deringer,
Rokas Elijošius,
Zakariya El-Machachi,
Fabio Falcioni,
Edvin Fako,
Andrea C. Ferrari,
Annalena Genreith-Schriever
, et al. (51 additional authors not shown)
Abstract:
Machine-learned force fields have transformed the atomistic modelling of materials by enabling simulations of ab initio quality on unprecedented time and length scales. However, they are currently limited by: (i) the significant computational and human effort that must go into development and validation of potentials for each particular system of interest; and (ii) a general lack of transferabilit…
▽ More
Machine-learned force fields have transformed the atomistic modelling of materials by enabling simulations of ab initio quality on unprecedented time and length scales. However, they are currently limited by: (i) the significant computational and human effort that must go into development and validation of potentials for each particular system of interest; and (ii) a general lack of transferability from one chemical system to the next. Here, using the state-of-the-art MACE architecture we introduce a single general-purpose ML model, trained on a public database of 150k inorganic crystals, that is capable of running stable molecular dynamics on molecules and materials. We demonstrate the power of the MACE-MP-0 model - and its qualitative and at times quantitative accuracy - on a diverse set problems in the physical sciences, including the properties of solids, liquids, gases, chemical reactions, interfaces and even the dynamics of a small protein. The model can be applied out of the box and as a starting or "foundation model" for any atomistic system of interest and is thus a step towards democratising the revolution of ML force fields by lowering the barriers to entry.
△ Less
Submitted 1 March, 2024; v1 submitted 29 December, 2023;
originally announced January 2024.
-
ACEpotentials.jl: A Julia Implementation of the Atomic Cluster Expansion
Authors:
William C. Witt,
Cas van der Oord,
Elena Gelžinytė,
Teemu Järvinen,
Andres Ross,
James P. Darby,
Cheuk Hin Ho,
William J. Baldwin,
Matthias Sachs,
James Kermode,
Noam Bernstein,
Gábor Csányi,
Christoph Ortner
Abstract:
We introduce ACEpotentials.jl, a Julia-language software package that constructs interatomic potentials from quantum mechanical reference data using the Atomic Cluster Expansion (Drautz, 2019). As the latter provides a complete description of atomic environments, including invariance to overall translation and rotation as well as permutation of like atoms, the resulting potentials are systematical…
▽ More
We introduce ACEpotentials.jl, a Julia-language software package that constructs interatomic potentials from quantum mechanical reference data using the Atomic Cluster Expansion (Drautz, 2019). As the latter provides a complete description of atomic environments, including invariance to overall translation and rotation as well as permutation of like atoms, the resulting potentials are systematically improvable and data efficient. Furthermore, the descriptor's expressiveness enables use of a linear model, facilitating rapid evaluation and straightforward application of Bayesian techniques for active learning. We summarize the capabilities of ACEpotentials.jl and demonstrate its strengths (simplicity, interpretability, robustness, performance) on a selection of prototypical atomistic modelling workflows.
△ Less
Submitted 7 September, 2023; v1 submitted 6 September, 2023;
originally announced September 2023.
-
wfl Python Toolkit for Creating Machine Learning Interatomic Potentials and Related Atomistic Simulation Workflows
Authors:
Elena Gelžinytė,
Simon Wengert,
Tamás K. Stenczel,
Hendrik H. Heenen,
Karsten Reuter,
Gábor Csányi,
Noam Bernstein
Abstract:
Predictive atomistic simulations are increasingly employed for data intensive high throughput studies that take advantage of constantly growing computational resources. To handle the sheer number of individual calculations that are needed in such studies, workflow management packages for atomistic simulations have been developed for a rapidly growing user base. These packages are predominantly des…
▽ More
Predictive atomistic simulations are increasingly employed for data intensive high throughput studies that take advantage of constantly growing computational resources. To handle the sheer number of individual calculations that are needed in such studies, workflow management packages for atomistic simulations have been developed for a rapidly growing user base. These packages are predominantly designed to handle computationally heavy ab initio calculations, usually with a focus on data provenance and reproducibility. However, in related simulation communities, e.g. the developers of machine learning interatomic potentials (MLIPs), the computational requirements are somewhat different: the types, sizes, and numbers of computational tasks are more diverse, and therefore require additional ways of parallelization and local or remote execution for optimal efficiency. In this work, we present the atomistic simulation and MLIP fitting workflow management package wfl and Python remote execution package ExPyRe to meet these requirements. With wfl and ExPyRe, versatile Atomic Simulation Environment based workflows that perform diverse procedures can be written. This capability is based on a low-level developer-oriented framework, which can be utilized to construct high level functionality for user-friendly programs. Such high level capabilities to automate machine learning interatomic potential fitting procedures are already incorporated in wfl, which we use to showcase its capabilities in this work. We believe that wfl fills an important niche in several growing simulation communities and will aid the development of efficient custom computational tasks.
△ Less
Submitted 1 August, 2023; v1 submitted 20 June, 2023;
originally announced June 2023.
-
The Ferris ferromagnetic resonance technique: principles and applications
Authors:
Amit Rothschild,
Benjamin Assouline,
Nadav Am Shalom,
Nirel Bernstein,
Goni Daniel,
Gil Cohen,
Amir Capua
Abstract:
Measurements of ferromagnetic resonance (FMR) are pivotal to modern magnetism and spintronics. Recently, we reported on the Ferris FMR technique, which relies on large-amplitude modulation of the externally applied magnetic field. It was shown to benefit from high sensitivity while being broadband. The Ferris FMR also expanded the resonance linewidth such that the sensitivity to spin currents was…
▽ More
Measurements of ferromagnetic resonance (FMR) are pivotal to modern magnetism and spintronics. Recently, we reported on the Ferris FMR technique, which relies on large-amplitude modulation of the externally applied magnetic field. It was shown to benefit from high sensitivity while being broadband. The Ferris FMR also expanded the resonance linewidth such that the sensitivity to spin currents was enhanced as well. Eventually, the spin Hall angle (θ_SH) was measurable even in wafer-level measurements that require low current densities to reduce the Joule heating. Despite the various advantages, analysis of the Ferris FMR response is limited to numerical modeling where the linewidth depends on multiple factors such as the field modulation profile and the magnetization saturation. Here, we describe in detail the basic principles of operation of the Ferris FMR and discuss its applicability and engineering considerations. We demonstrated these principles in a measurement of the orbital Hall effect taking place in Cu, using an Au layer as the orbital to spin current converter. This illustrates the potential of the Ferris FMR for the future development of spintronics technology.
△ Less
Submitted 31 May, 2023;
originally announced June 2023.
-
Efficient generation of spin currents by the Orbital Hall effect in pure Cu and Al and their measurement by a Ferris-wheel ferromagnetic resonance technique at the wafer level
Authors:
Amit Rothschild,
Nadav Am-Shalom,
Nirel Bernstein,
Mayan Meron,
Tal David,
Benjamin Assouline,
Elichai Frohlich,
Jiewen Xiao,
Binghai Yan,
Amir Capua
Abstract:
We present a new ferromagnetic resonance (FMR) method that we term the Ferris FMR. It is wideband, has significantly higher sensitivity as compared to conventional FMR systems, and measures the absorption line rather than its derivative. It is based on large-amplitude modulation of the externally applied magnetic field that effectively magnifies signatures of the spin-transfer torque making its me…
▽ More
We present a new ferromagnetic resonance (FMR) method that we term the Ferris FMR. It is wideband, has significantly higher sensitivity as compared to conventional FMR systems, and measures the absorption line rather than its derivative. It is based on large-amplitude modulation of the externally applied magnetic field that effectively magnifies signatures of the spin-transfer torque making its measurement possible even at the wafer-level. Using the Ferris FMR, we report on the generation of spin currents from the orbital Hall effect taking place in pure Cu and Al. To this end, we use the spin-orbit coupling of a thin Pt layer introduced at the interface that converts the orbital current to a measurable spin current. While Cu reveals a large effective spin Hall angle exceeding that of Pt, Al possesses an orbital Hall effect of opposite polarity in agreement with the theoretical predictions. Our results demonstrate additional spin- and orbit- functionality for two important metals in the semiconductor industry beyond their primary use as interconnects with all the advantages in power, scaling, and cost.
△ Less
Submitted 21 August, 2022;
originally announced August 2022.
-
Nested sampling for physical scientists
Authors:
Greg Ashton,
Noam Bernstein,
Johannes Buchner,
Xi Chen,
Gábor Csányi,
Andrew Fowlie,
Farhan Feroz,
Matthew Griffiths,
Will Handley,
Michael Habeck,
Edward Higson,
Michael Hobson,
Anthony Lasenby,
David Parkinson,
Livia B. Pártay,
Matthew Pitkin,
Doris Schneider,
Joshua S. Speagle,
Leah South,
John Veitch,
Philipp Wacker,
David J. Wales,
David Yallup
Abstract:
We review Skilling's nested sampling (NS) algorithm for Bayesian inference and more broadly multi-dimensional integration. After recapitulating the principles of NS, we survey developments in implementing efficient NS algorithms in practice in high-dimensions, including methods for sampling from the so-called constrained prior. We outline the ways in which NS may be applied and describe the applic…
▽ More
We review Skilling's nested sampling (NS) algorithm for Bayesian inference and more broadly multi-dimensional integration. After recapitulating the principles of NS, we survey developments in implementing efficient NS algorithms in practice in high-dimensions, including methods for sampling from the so-called constrained prior. We outline the ways in which NS may be applied and describe the application of NS in three scientific fields in which the algorithm has proved to be useful: cosmology, gravitational-wave astronomy, and materials science. We close by making recommendations for best practice when using NS and by summarizing potential limitations and optimizations of NS.
△ Less
Submitted 31 May, 2022;
originally announced May 2022.
-
Spin-current amplification by stimulated spin pumping
Authors:
Benjamin Assouline,
Marina Brik,
Nirel Bernstein,
Amir Capua
Abstract:
Stimulated emission is a process during which an atomic system gives away energy to create a coherent photon. It is fundamental to the operation of the optical amplifier. Here we propose two mechanisms for amplifying AC spin currents in a solid-state magnetic medium by a stimulated spin pumping process. The first is synchronous and consists of phase-locked pulses that perturb a precessing magnetic…
▽ More
Stimulated emission is a process during which an atomic system gives away energy to create a coherent photon. It is fundamental to the operation of the optical amplifier. Here we propose two mechanisms for amplifying AC spin currents in a solid-state magnetic medium by a stimulated spin pumping process. The first is synchronous and consists of phase-locked pulses that perturb a precessing magnetic moment. The second is asynchronous and is driven by DC spin currents. The amplification relies on a non-adiabatic interaction taking place in a ferromagnetic medium in which the magnetic moment emits spin angular momentum in the form of spin current before equilibrating with the environment. The pumped spin current amplifies or absorbs the injected AC spin current mimicking the operation of the optical gain medium as readily seen from the gain saturation profiles. The mechanisms we propose are a first step towards a realistic spin current amplifier.
△ Less
Submitted 5 September, 2021;
originally announced September 2021.
-
Determination of the spin Hall angle by the inverse spin Hall effect, device level ferromagnetic resonance, and spin torque ferromagnetic resonance: a comparison of methods
Authors:
Ranen Ben-Shalom,
Nirel Bernstein,
See-Hun Yang,
Amir Capua
Abstract:
The spin torque ferromagnetic resonance (STFMR) is one of the popular methods for measurement of the spin Hall angle (SHA). However, in order to accurately determine SHA from STFMR measurements, the acquired data must be carefully analyzed: The resonance linewidth should be determined to an accuracy of a fraction of an Oe, while the dynamical interaction leading to the measured response consists o…
▽ More
The spin torque ferromagnetic resonance (STFMR) is one of the popular methods for measurement of the spin Hall angle (SHA). However, in order to accurately determine SHA from STFMR measurements, the acquired data must be carefully analyzed: The resonance linewidth should be determined to an accuracy of a fraction of an Oe, while the dynamical interaction leading to the measured response consists of the conventional field-induced ferromagnetic resonance (FMR), spin-torque induced FMR, and of the inverse spin Hall effect (ISHE). Additionally, the signal often deteriorates when DC current is passed through the device. In this work we compare the STFMR method with two other FMR-based methods that are used to extract SHA. The first is a device-level FMR and the second is based on the ISHE. We identify artefacts that are caused by the noise floor of the instrumentation that make the measurement of SHA illusive even when the signal to noise ratio seems to be reasonable. Additionally, we estimate a 10% error in SHA that results from neglecting the magnetic anisotropies as in conventional measurements. Overall, we find the STFMR to be the most robust of the three methods despite the complexity of the interaction taking place therein. The conclusions of our work lead to a more accurate determination of SHA and will assist in the search of novel materials for energy efficient spin-based applications.
△ Less
Submitted 19 May, 2021;
originally announced May 2021.
-
Microwave and spin transfer torque driven coherent control in ferromagnets
Authors:
Marina Brik,
Nirel Bernstein,
Amir Capua
Abstract:
Coherent control is a method used to manipulate the state of matter using oscillatory electromagnetic radiation which relies on the non-adiabatic interaction. It is commonly applied in quantum processing applications. This technique is interesting in the context of ferromagnetic materials because of the ability to combine it with spintronics for the purpose of fundamental spin transport research,…
▽ More
Coherent control is a method used to manipulate the state of matter using oscillatory electromagnetic radiation which relies on the non-adiabatic interaction. It is commonly applied in quantum processing applications. This technique is interesting in the context of ferromagnetic materials because of the ability to combine it with spintronics for the purpose of fundamental spin transport research, low-power information processing, and potentially future quantum bit (Qubit) applications. In this work we address the theoretical grounds of coherent manipulation in practical ferromagnetic systems. We study electromagnetic radiation driven interaction that is enhanced in the presence of spin polarized currents and map the conditions that allow coherent manipulation for which Rabi oscillations take place. The role of the magnetic anisotropy field is shown to act as an additional oscillatory driving field. We discuss the Gilbert losses in the context of effective coherence decay rates and show that it is possible to control these rates by application of a static spin current. The case of coherent manipulation using oscillatory spin currents that is free of radiation is discussed as well. Our work paves the way towards spin current amplification as well as radiation-free coherent control schemes that may potentially lead to novel Qubits that are robust and scalable.
△ Less
Submitted 22 September, 2020; v1 submitted 30 August, 2020;
originally announced August 2020.
-
Structural transitions in dense disordered silicon from quantum-accurate ultra-large-scale simulations
Authors:
Volker L. Deringer,
Noam Bernstein,
Gábor Csányi,
Mark Wilson,
David A. Drabold,
Stephen R. Elliott
Abstract:
Structurally disordered materials continue to pose fundamental questions, including that of how different disordered phases ("polyamorphs") can coexist and transform from one to another. As a widely studied case, amorphous silicon (a-Si) forms a fourfold-coordinated, covalent random network at ambient conditions, but much higher-coordinated, metallic-like phases under pressure. However, a detailed…
▽ More
Structurally disordered materials continue to pose fundamental questions, including that of how different disordered phases ("polyamorphs") can coexist and transform from one to another. As a widely studied case, amorphous silicon (a-Si) forms a fourfold-coordinated, covalent random network at ambient conditions, but much higher-coordinated, metallic-like phases under pressure. However, a detailed mechanistic understanding of the liquid-amorphous and amorphous-amorphous transitions in silicon has been lacking, due to intrinsic limitations of even the most advanced experimental and computational techniques. Here, we show how machine-learning (ML)-driven simulations can break through this long-standing barrier, affording a comprehensive, quantum-accurate, and fully atomistic description of all relevant liquid and amorphous phases of silicon. Combining a model system size of 100,000 atoms (ten-nanometre length scale) with a prediction accuracy of a few meV per atom, our simulations reveal a remarkable, three-step transformation sequence for a-Si under increasing external pressure. First, up to 10-11 GPa, polyamorphic low- and high-density amorphous (LDA and HDA) regions are found to coexist, rather than appearing sequentially. Then, we observe a structural collapse into a distinct, very-high-density amorphous (VHDA) phase at 12-13 GPa, reminiscent of the dense liquid but being formed at a much lower temperature. Finally, our simulations indicate the transient nature of this VHDA phase: it rapidly nucleates crystallites at 13-16 GPa, ultimately leading to the formation of a poly-crystalline, simple-hexagonal structure, consistent with experiments but not seen in earlier simulations.
△ Less
Submitted 16 December, 2019;
originally announced December 2019.
-
Impact of biaxial and uniaxial strain on V$_2$O$_3$
Authors:
Darshana Wickramaratne,
Noam Bernstein,
I. I. Mazin
Abstract:
Using first-principles calculations we determine the role of compressive and tensile uniaxial and equibiaxial strain on the structural, electronic and magnetic properties of V$_2$O$_3$. We find that compressive strain increases the energy cost to transition from the high-temperature paramagnetic metallic phase to the low-temperature antiferromagnetic insulating phase. This shift in the energy diff…
▽ More
Using first-principles calculations we determine the role of compressive and tensile uniaxial and equibiaxial strain on the structural, electronic and magnetic properties of V$_2$O$_3$. We find that compressive strain increases the energy cost to transition from the high-temperature paramagnetic metallic phase to the low-temperature antiferromagnetic insulating phase. This shift in the energy difference can be explained by changes in the V-V bond lengths that are antiferromagnetically aligned in the low temperature structure. The insights that we have obtained provide a microscopic explanation for the shifts in the metal-insulator transition temperature that have been observed in experiments of V$_2$O$_3$ films grown on different substrates.
△ Less
Submitted 29 September, 2019;
originally announced September 2019.
-
Machine-learned Interatomic Potentials for Alloys and Alloy Phase Diagrams
Authors:
Conrad W. Rosenbrock,
Konstantin Gubaev,
Alexander V. Shapeev,
Livia B. Pártay,
Noam Bernstein,
Gábor Csányi,
Gus L. W. Hart
Abstract:
We introduce machine-learned potentials for Ag-Pd to describe the energy of alloy configurations over a wide range of compositions. We compare two different approaches. Moment tensor potentials (MTP) are polynomial-like functions of interatomic distances and angles. The Gaussian Approximation Potential (GAP) framework uses kernel regression, and we use the Smooth Overlap of Atomic Positions (SOAP)…
▽ More
We introduce machine-learned potentials for Ag-Pd to describe the energy of alloy configurations over a wide range of compositions. We compare two different approaches. Moment tensor potentials (MTP) are polynomial-like functions of interatomic distances and angles. The Gaussian Approximation Potential (GAP) framework uses kernel regression, and we use the Smooth Overlap of Atomic Positions (SOAP) representation of atomic neighbourhoods that consists of a complete set of rotational and permutational invariants provided by the power spectrum of the spherical Fourier transform of the neighbour density. Both types of potentials give excellent accuracy for a wide range of compositions and rival the accuracy of cluster expansion, a benchmark for this system. While both models are able to describe small deformations away from the lattice positions, SOAP-GAP excels at transferability as shown by sensible transformation paths between configurations, and MTP allows, due to its lower computational cost, the calculation of compositional phase diagrams. Given the fact that both methods perform as well as cluster expansion would but yield off-lattice models, we expect them to open new avenues in computational materials modeling for alloys.
△ Less
Submitted 9 July, 2019; v1 submitted 18 June, 2019;
originally announced June 2019.
-
De novo exploration and self-guided learning of potential-energy surfaces
Authors:
Noam Bernstein,
Gábor Csányi,
Volker L. Deringer
Abstract:
Interatomic potential models based on machine learning (ML) are rapidly developing as tools for materials simulations. However, because of their flexibility, they require large fitting databases that are normally created with substantial manual selection and tuning of reference configurations. Here, we show that ML potentials can be built in a largely automated fashion, exploring and fitting poten…
▽ More
Interatomic potential models based on machine learning (ML) are rapidly developing as tools for materials simulations. However, because of their flexibility, they require large fitting databases that are normally created with substantial manual selection and tuning of reference configurations. Here, we show that ML potentials can be built in a largely automated fashion, exploring and fitting potential-energy surfaces from the beginning (de novo) within one and the same protocol. The key enabling step is the use of a configuration-averaged kernel metric that allows one to select the few most relevant structures at each step. The resulting potentials are accurate and robust for the wide range of configurations that occur during structure searching, despite only requiring a relatively small number of single-point DFT calculations on small unit cells. We apply the method to materials with diverse chemical nature and coordination environments, marking a milestone toward the more routine application of ML potentials in physics, chemistry, and materials science.
△ Less
Submitted 24 May, 2019;
originally announced May 2019.
-
The role of defects in the metal-insulator transition in VO$_{2}$ and V$_2$O$_3$
Authors:
Darshana Wickramaratne,
Noam Bernstein,
I. I. Mazin
Abstract:
The vanadates VO$_2$ and V$_2$O$_3$ are prototypical examples of strongly correlated materials that exhibit a metal-insulator transition. While the phase transitions in these materials have been studied extensively, there is a limited understanding of how the properties of these materials are affected by the presence of defects and doping. In this study we investigate the impact of native point de…
▽ More
The vanadates VO$_2$ and V$_2$O$_3$ are prototypical examples of strongly correlated materials that exhibit a metal-insulator transition. While the phase transitions in these materials have been studied extensively, there is a limited understanding of how the properties of these materials are affected by the presence of defects and doping. In this study we investigate the impact of native point defects in the form of Frenkel defects on the structural, magnetic and electronic properties of VO$_2$ and V$_2$O$_3$, using first-principles calculations. In VO$_2$ the vanadium Frenkel pairs lead to a non-trivial insulating state. The unpaired vanadium interstitial bonds to a single dimer, which leads to a trimer that has one singlet state and one localized single-electron $S=1/2$ state. The unpaired broken dimer created by the vanadium vacancy also has a localized $S=1/2$ state. Thus, the insulating state is created by the singlet dimers, the trimer and the two localized $S=1/2$ states. Oxygen Frenkel pairs, on the other hand, lead to a metallic state in VO$_2$, but are expected to be present in much lower concentrations. In contrast, the Frenkel defects in V$_2$O$_3$ do not directly suppress the insulating character of the material. However, the disorder created by defects in V$_2$O$_3$ alters the local magnetic moments and in turn reduces the energy cost of a transition between the insulating and conducting phases of the material. We also find self-trapped small polarons in V$_2$O$_3$, which has implications for transport properties in the insulating phase.
△ Less
Submitted 12 February, 2019;
originally announced February 2019.
-
Quantifying Chemical Structure and Atomic Energies in Amorphous Silicon Networks
Authors:
Noam Bernstein,
Bishal Bhattarai,
Gábor Csányi,
David A. Drabold,
Stephen R. Elliott,
Volker L. Deringer
Abstract:
Amorphous materials are coming within reach of realistic computer simulations, but new approaches are needed to fully understand their intricate atomic structures. Here, we show how machine-learning (ML)-based techniques can give new, quantitative chemical insight into the atomic-scale structure of amorphous silicon (a-Si). Based on a similarity function ("kernel"), we define a structural metric t…
▽ More
Amorphous materials are coming within reach of realistic computer simulations, but new approaches are needed to fully understand their intricate atomic structures. Here, we show how machine-learning (ML)-based techniques can give new, quantitative chemical insight into the atomic-scale structure of amorphous silicon (a-Si). Based on a similarity function ("kernel"), we define a structural metric that unifies the description of nearest- and next-nearest-neighbor environments in the amorphous state. We apply this to an ensemble of a-Si networks, generated in melt-quench simulations with an ML-based interatomic potential, in which we tailor the degree of ordering by varying the quench rates down to $10^{10}$ K/s (leading to a structural model that is lower in energy than the established WWW network). We then show how "machine-learned" atomic energies permit a chemical interpretation, associating coordination defects in a-Si with distinct energetic stability regions. The approach is straightforward and inexpensive to apply to arbitrary structural models, and it is therefore expected to have more general significance for developing a quantitative understanding of the amorphous state.
△ Less
Submitted 27 November, 2018;
originally announced November 2018.
-
Machine learning a general purpose interatomic potential for silicon
Authors:
Albert P. Bartok,
James Kermode,
Noam Bernstein,
Gabor Csanyi
Abstract:
The success of first principles electronic structure calculation for predictive modeling in chemistry, solid state physics, and materials science is constrained by the limitations on simulated length and time scales due to computational cost and its scaling. Techniques based on machine learning ideas for interpolating the Born-Oppenheimer potential energy surface without explicitly describing elec…
▽ More
The success of first principles electronic structure calculation for predictive modeling in chemistry, solid state physics, and materials science is constrained by the limitations on simulated length and time scales due to computational cost and its scaling. Techniques based on machine learning ideas for interpolating the Born-Oppenheimer potential energy surface without explicitly describing electrons have recently shown great promise, but accurately and efficiently fitting the physically relevant space of configurations has remained a challenging goal. Here we present a Gaussian Approximation Potential for silicon that achieves this milestone, accurately reproducing density functional theory reference results for a wide range of observable properties, including crystal, liquid, and amorphous bulk phases, as well as point, line, and plane defects. We demonstrate that this new potential enables calculations that would be extremely expensive with a first principles electronic structure method, such as finite temperature phase boundary lines, self-diffusivity in the liquid, formation of the amorphous by slow quench, and dynamic brittle fracture. We show that the uncertainty quantification inherent to the Gaussian process regression framework gives a qualitative estimate of the potential's accuracy for a given atomic configuration. The success of this model shows that it is indeed possible to create a useful machine-learning-based interatomic potential that comprehensively describes a material, and serves as a template for the development of such models in the future.
△ Less
Submitted 3 May, 2018;
originally announced May 2018.
-
Realistic atomistic structure of amorphous silicon from machine-learning-driven molecular dynamics
Authors:
Volker L. Deringer,
Noam Bernstein,
Albert P. Bartók,
Matthew J. Cliffe,
Rachel N. Kerber,
Lauren E. Marbella,
Clare P. Grey,
Stephen R. Elliott,
Gábor Csányi
Abstract:
Amorphous silicon (a-Si) is a widely studied non-crystalline material, and yet the subtle details of its atomistic structure are still unclear. Here, we show that accurate structural models of a-Si can be obtained by harnessing the power of machine-learning algorithms to create interatomic potentials. Our best a-Si network is obtained by cooling from the melt in molecular-dynamics simulations, at…
▽ More
Amorphous silicon (a-Si) is a widely studied non-crystalline material, and yet the subtle details of its atomistic structure are still unclear. Here, we show that accurate structural models of a-Si can be obtained by harnessing the power of machine-learning algorithms to create interatomic potentials. Our best a-Si network is obtained by cooling from the melt in molecular-dynamics simulations, at a rate of 10$^{11}$ K/s (that is, on the 10 ns timescale). This structure shows a defect concentration of below 2% and agrees with experiments regarding excess energies, diffraction data, as well as $^{29}$Si solid-state NMR chemical shifts. We show that this level of quality is impossible to achieve with faster quench simulations. We then generate a 4,096-atom system which correctly reproduces the magnitude of the first sharp diffraction peak (FSDP) in the structure factor, achieving the closest agreement with experiments to date. Our study demonstrates the broader impact of machine-learning interatomic potentials for elucidating accurate structures and properties of amorphous functional materials.
△ Less
Submitted 7 March, 2018;
originally announced March 2018.
-
Uncovering the mechanism of the impurity-selective Mott transition in paramagnetic V$_{2}$O$_{3}$
Authors:
Frank Lechermann,
Noam Bernstein,
I. I. Mazin,
Roser Valentí
Abstract:
While the phase diagrams of the one- and multi-orbital Hubbard model have been well studied, the physics of real Mott insulators is often much richer, material dependent, and poorly understood. In the prototype Mott insulator V$_{2}$O$_{3}$, chemical pressure was initially believed to explain why the paramagnetic-metal to antiferromagnetic-insulator transition temperature is lowered by Ti doping w…
▽ More
While the phase diagrams of the one- and multi-orbital Hubbard model have been well studied, the physics of real Mott insulators is often much richer, material dependent, and poorly understood. In the prototype Mott insulator V$_{2}$O$_{3}$, chemical pressure was initially believed to explain why the paramagnetic-metal to antiferromagnetic-insulator transition temperature is lowered by Ti doping while Cr doping strengthens correlations, eventually rendering the high-temperature phase paramagnetic insulating. However, this scenario has been recently shown both experimentally and theoretically to be untenable. Based on full structural optimization, we demonstrate via the charge self-consistent combination of density functional theory and dynamical mean-field theory that changes in the V$_{2}$O$_{3}$ phase diagram are driven by defect-induced local symmetry breakings resulting from dramatically different couplings of Cr and Ti dopants to the host system. This finding emphasizes the high sensitivity of the Mott metal-insulator transition to the local environment and the importance of accurately accounting for the one-electron Hamiltonian, since correlations crucially respond to it.
△ Less
Submitted 23 May, 2018; v1 submitted 26 January, 2018;
originally announced January 2018.
-
Constant-pressure nested sampling with atomistic dynamics
Authors:
Robert J. N. Baldock,
Noam Bernstein,
K. Michael Salerno,
Lívia B. Pártay,
Gábor Csányi
Abstract:
The nested sampling algorithm has been shown to be a general method for calculating the pressure-temperature-composition phase diagrams of materials. While the previous implementation used single-particle Monte Carlo moves, these are inefficient for condensed systems with general interactions where single-particle moves cannot be evaluated faster than the energy of the whole system. Here we enhanc…
▽ More
The nested sampling algorithm has been shown to be a general method for calculating the pressure-temperature-composition phase diagrams of materials. While the previous implementation used single-particle Monte Carlo moves, these are inefficient for condensed systems with general interactions where single-particle moves cannot be evaluated faster than the energy of the whole system. Here we enhance the method by using all-particle moves: either Galilean Monte Carlo or a total enthalpy Hamiltonian Monte Carlo algorithm, introduced in this paper. We show that these algorithms enable the determination of phase transition temperatures with equivalent accuracy to the previous method at $1/N$ of the cost for an $N$-particle system with general interactions, or at equal cost when single particle moves can be done in $1/N$ of the cost of a full $N$-particle energy evaluation.
△ Less
Submitted 16 November, 2017; v1 submitted 30 October, 2017;
originally announced October 2017.
-
Bright triplet excitons in lead halide perovskites
Authors:
Michael A. Becker,
Roman Vaxenburg,
Georgian Nedelcu,
Peter C. Sercel,
Andrew Shabaev,
Michael J. Mehl,
John G. Michopoulos,
Samuel G. Lambrakos,
Noam Bernstein,
John L. Lyons,
Thilo Stöferle,
Rainer F. Mahrt,
Maksym V. Kovalenko,
David J. Norris,
Gabriele Rainò,
Alexander L. Efros
Abstract:
Nanostructured semiconductors emit light from electronic states known as excitons[1]. According to Hund's rules[2], the lowest energy exciton in organic materials should be a poorly emitting triplet state. Analogously, the lowest exciton level in all known inorganic semiconductors is believed to be optically inactive. These 'dark' excitons (into which the system can relax) hinder light-emitting de…
▽ More
Nanostructured semiconductors emit light from electronic states known as excitons[1]. According to Hund's rules[2], the lowest energy exciton in organic materials should be a poorly emitting triplet state. Analogously, the lowest exciton level in all known inorganic semiconductors is believed to be optically inactive. These 'dark' excitons (into which the system can relax) hinder light-emitting devices based on semiconductor nanostructures. While strategies to diminish their influence have been developed[3-5], no materials have been identified in which the lowest exciton is bright. Here we show that the lowest exciton in quasi-cubic lead halide perovskites is optically active. We first use the effective-mass model and group theory to explore this possibility, which can occur when the strong spin-orbit coupling in the perovskite conduction band is combined with the Rashba effect [6-10]. We then apply our model to CsPbX3 (X=Cl,Br,I) nanocrystals[11], for which we measure size- and composition-dependent fluorescence at the single-nanocrystal level. The bright character of the lowest exciton immediately explains the anomalous photon-emission rates of these materials, which emit 20 and 1,000 times faster[12] than any other semiconductor nanocrystal at room[13-16] and cryogenic[17] temperatures, respectively. The bright exciton is further confirmed by detailed analysis of the fine structure in low-temperature fluorescence spectra. For semiconductor nanocrystals[18], which are already used in lighting[19,20], lasers[21,22], and displays[23], these optically active excitons can lead to materials with brighter emission and enhanced absorption. More generally, our results provide criteria for identifying other semiconductors exhibiting bright excitons with potentially broad implications for optoelectronic devices.
△ Less
Submitted 10 July, 2017;
originally announced July 2017.
-
Machine Learning Unifies the Modelling of Materials and Molecules
Authors:
Albert P. Bartok,
Sandip De,
Carl Poelking,
Noam Bernstein,
James Kermode,
Gabor Csanyi,
Michele Ceriotti
Abstract:
Determining the stability of molecules and condensed phases is the cornerstone of atomistic modelling, underpinning our understanding of chemical and materials properties and transformations. Here we show that a machine learning model, based on a local description of chemical environments and Bayesian statistical learning, provides a unified framework to predict atomic-scale properties. It capture…
▽ More
Determining the stability of molecules and condensed phases is the cornerstone of atomistic modelling, underpinning our understanding of chemical and materials properties and transformations. Here we show that a machine learning model, based on a local description of chemical environments and Bayesian statistical learning, provides a unified framework to predict atomic-scale properties. It captures the quantum mechanical effects governing the complex surface reconstructions of silicon, predicts the stability of different classes of molecules with chemical accuracy, and distinguishes active and inactive protein ligands with more than 99% reliability. The universality and the systematic nature of our framework provides new insight into the potential energy surface of materials and molecules.
△ Less
Submitted 15 December, 2017; v1 submitted 1 June, 2017;
originally announced June 2017.
-
Modulus-Pressure Equation for Confined Fluids
Authors:
Gennady Y. Gor,
Daniel W. Siderius,
Vincent K. Shen,
Noam Bernstein
Abstract:
Ultrasonic experiments allow one to measure the elastic modulus of bulk solid or fluid samples. Recently such experiments have been carried out on fluid-saturated nanoporous glass to probe the modulus of a confined fluid. In our previous work [J. Chem. Phys., (2015) 143, 194506], using Monte Carlo simulations we showed that the elastic modulus $K$ of a fluid confined in a mesopore is a function of…
▽ More
Ultrasonic experiments allow one to measure the elastic modulus of bulk solid or fluid samples. Recently such experiments have been carried out on fluid-saturated nanoporous glass to probe the modulus of a confined fluid. In our previous work [J. Chem. Phys., (2015) 143, 194506], using Monte Carlo simulations we showed that the elastic modulus $K$ of a fluid confined in a mesopore is a function of the pore size. Here we focus on modulus-pressure dependence $K(P)$, which is linear for bulk materials, a relation known as the Tait-Murnaghan equation. Using transition-matrix Monte Carlo simulations we calculated the elastic modulus of bulk argon as a function of pressure and argon confined in silica mesopores as a function of Laplace pressure. Our calculations show that while the elastic modulus is strongly affected by confinement and temperature, the slope of the modulus versus pressure is not. Moreover, the calculated slope is in a good agreement with the reference data for bulk argon and experimental data for confined argon derived from ultrasonic experiments. We propose to use the value of the slope of $K(P)$ to estimate the elastic moduli of an unknown porous medium.
△ Less
Submitted 23 August, 2016;
originally announced August 2016.
-
A universal preconditioner for simulating condensed phase materials
Authors:
David Packwood,
James Kermode,
Letif Mones,
Noam Bernstein,
John Woolley,
Nicholas Gould,
Christoph Ortner,
Gabor Csanyi
Abstract:
We introduce a universal sparse preconditioner that accelerates geometry optimisation and saddle point search tasks that are common in the atomic scale simulation of materials. Our preconditioner is based on the neighbourhood structure and we demonstrate the gain in computational efficiency in a wide range of materials that include metals, insulators and molecular solids. The simple structure of t…
▽ More
We introduce a universal sparse preconditioner that accelerates geometry optimisation and saddle point search tasks that are common in the atomic scale simulation of materials. Our preconditioner is based on the neighbourhood structure and we demonstrate the gain in computational efficiency in a wide range of materials that include metals, insulators and molecular solids. The simple structure of the preconditioner means that the gains can be realised in practice not only when using expensive electronic structure models but also for fast empirical potentials. Even for relatively small systems of a few hundred atoms, we observe speedups of a factor of two or more, and the gain grows with system size. An open source Python implementation within the Atomic Simulation Environment is available, offering interfaces to a wide range of atomistic codes.
△ Less
Submitted 9 April, 2016;
originally announced April 2016.
-
Relation Between Pore Size and the Compressibility of a Confined Fluid
Authors:
Gennady Y. Gor,
Daniel W. Siderius,
Christopher J. Rasmussen,
William P. Krekelberg,
Vincent K. Shen,
Noam Bernstein
Abstract:
When a fluid is confined to a nanopore, its thermodynamic properties differ from the properties of a bulk fluid, so measuring such properties of the confined fluid can provide information about the pore sizes. Here we report a simple relation between the pore size and isothermal compressibility of argon confined in these pores. Compressibility is calculated from the fluctuations of the number of p…
▽ More
When a fluid is confined to a nanopore, its thermodynamic properties differ from the properties of a bulk fluid, so measuring such properties of the confined fluid can provide information about the pore sizes. Here we report a simple relation between the pore size and isothermal compressibility of argon confined in these pores. Compressibility is calculated from the fluctuations of the number of particles in the grand canonical ensemble using two different simulation techniques: conventional grand-canonical Monte Carlo and grand-canonical ensemble transition-matrix Monte Carlo. Our results provide a theoretical framework for extracting the information on the pore sizes of fluid-saturated samples by measuring the compressibility from ultrasonic experiments.
△ Less
Submitted 3 November, 2015;
originally announced November 2015.
-
Elastic Response of Mesoporous Silicon to Capillary Pressures in the Pores
Authors:
Gennady Y. Gor,
Luca Bertinetti,
Noam Bernstein,
Peter Fratzl,
Patrick Huber
Abstract:
We study water adsorption-induced deformation of a monolithic, mesoporous silicon membrane traversed by independent channels of $\sim$8 nm diameter. We focus on the elastic constant associated with the Laplace pressure-induced deformation of the membrane upon capillary condensation, i.e. the pore-load modulus. We perform finite-element method (FEM) simulations of the adsorption-induced deformation…
▽ More
We study water adsorption-induced deformation of a monolithic, mesoporous silicon membrane traversed by independent channels of $\sim$8 nm diameter. We focus on the elastic constant associated with the Laplace pressure-induced deformation of the membrane upon capillary condensation, i.e. the pore-load modulus. We perform finite-element method (FEM) simulations of the adsorption-induced deformation of hexagonal and square lattices of cylindrical pores representing the membrane. We find that the pore-load modulus weakly depends on the geometrical arrangement of pores, and can be expressed as a function of porosity. We propose an analytical model which relates the pore-load modulus to the porosity and to the elastic properties of bulk silicon (Young's modulus and Poisson's ratio), and provides an excellent agreement with FEM results. We find good agreement between our experimental data and the predictions of the analytical model, with the Young's modulus of the pore walls slightly lower than the bulk value. This model is applicable to a large class of materials with morphologies similar to mesoporous silicon. Moreover, our findings suggest that liquid condensation experiments allow one to elegantly access the elastic constants of a mesoporous medium.
△ Less
Submitted 23 July, 2015;
originally announced July 2015.
-
What superconducts in sulfur hydrides under pressure, and why
Authors:
N. Bernstein,
C. S. Hellberg,
M. D. Johannes,
I. I. Mazin,
M. J. Mehl
Abstract:
The recent discovery of superconductivity at 190~K in highly compressed H$_{2}$S is spectacular not only because it sets a record high critical temperature, but because it does so in a material that appears to be, and we argue here that it is, a conventional strong-coupling BCS superconductor. Intriguingly, superconductivity in the observed pressure and temperature range was predicted theoreticall…
▽ More
The recent discovery of superconductivity at 190~K in highly compressed H$_{2}$S is spectacular not only because it sets a record high critical temperature, but because it does so in a material that appears to be, and we argue here that it is, a conventional strong-coupling BCS superconductor. Intriguingly, superconductivity in the observed pressure and temperature range was predicted theoretically in a similar compound H$_{3}$S. Several important questions about this remarkable result, however, are left unanswered: (1) Does the stoichiometry of the superconducting compound differ from the nominal composition, and could it be the predicted H$_{3}$S compound? (2) Is the physical origin of the anomalously high critical temperature related only to the high H phonon frequencies, or does strong electron-ion coupling play a role? We show that at experimentally relevant pressures H$_2$S is unstable, decomposing into H$_3$S and S, and that H$_3$S has a record high $T_c$ due to its covalent bonds driven metallic. The main reason for this extraordinarily high $T_c$ in H$_3$S as compared with MgB$_2$, another compound with a similar superconductivity mechanism, is the high vibrational frequency of the much lighter H atoms.
△ Less
Submitted 31 December, 2014;
originally announced January 2015.
-
The Adaptive Buffered Force QM/MM method in the CP2K and AMBER software packages
Authors:
Letif Mones,
Andrew Jones,
Andreas W. Götz,
Teodoro Laino,
Ross C. Walker,
Ben Leimkuhler,
Gábor Csányi,
Noam Bernstein
Abstract:
The implementation and validation of the adaptive buffered force QM/MM method in two popular packages, CP2K and AMBER are presented. The implementations build on the existing QM/MM functionality in each code, extending it to allow for redefinition of the QM and MM regions during the simulation and reducing QM-MM interface errors by discarding forces near the boundary according to the buffered forc…
▽ More
The implementation and validation of the adaptive buffered force QM/MM method in two popular packages, CP2K and AMBER are presented. The implementations build on the existing QM/MM functionality in each code, extending it to allow for redefinition of the QM and MM regions during the simulation and reducing QM-MM interface errors by discarding forces near the boundary according to the buffered force-mixing approach. New adaptive thermostats, needed by force-mixing methods, are also implemented. Different variants of the method are benchmarked by simulating the structure of bulk water, water autoprotolysis in the presence of zinc and dimethyl-phosphate hydrolysis using various semiempirical Hamiltonians and density functional theory as the QM model. It is shown that with suitable parameters, based on force convergence tests, the adaptive buffered-force QM/MM scheme can provide an accurate approximation of the structure in the dynamical QM region matching the corresponding fully QM simulations, as well as reproducing the correct energetics in all cases. Adaptive unbuffered force-mixing and adaptive conventional QM/MM methods also provide reasonable results for some systems, but are more likely to suffer from instabilities and inaccuracies.
△ Less
Submitted 18 September, 2014;
originally announced September 2014.
-
Free energy surface reconstruction from umbrella samples using Gaussian process regression II: Multiple collective variables
Authors:
Noam Bernstein,
Thomas Stecher,
Gábor Csányi
Abstract:
We demonstrate how a prior assumption of smoothness can be used to enhance the reconstruction of free energy profiles from multiple umbrella sampling simulations using the Bayesian Gaussian process regression approach. The method we derive allows the concurrent use of histograms and free energy gradients and can easily be extended to include further data. In Part I we review the necessary theory a…
▽ More
We demonstrate how a prior assumption of smoothness can be used to enhance the reconstruction of free energy profiles from multiple umbrella sampling simulations using the Bayesian Gaussian process regression approach. The method we derive allows the concurrent use of histograms and free energy gradients and can easily be extended to include further data. In Part I we review the necessary theory and test the method for one collective variable. We demonstrate improved performance with respect to the weighted histogram analysis method and obtain meaningful error bars without any significant additional computation. In Part II we consider the case of multiple collective variables and compare to a reconstruction using least squares fitting of radial basis functions. We find substantial improvements in the regimes of spatially sparse data or short sampling trajectories. A software implementation is made available on www.libatoms.org.
△ Less
Submitted 24 July, 2014; v1 submitted 16 December, 2013;
originally announced December 2013.
-
Free energy surface reconstruction from umbrella samples using Gaussian process regression
Authors:
Thomas Stecher,
Noam Bernstein,
Gábor Csányi
Abstract:
We demonstrate how the Gaussian process regression approach can be used to efficiently reconstruct free energy surfaces from umbrella sampling simulations. By making a prior assumption of smoothness and taking account of the sampling noise in a consistent fashion, we achieve a significant improvement in accuracy over the state of the art in two or more dimensions or, equivalently, a significant co…
▽ More
We demonstrate how the Gaussian process regression approach can be used to efficiently reconstruct free energy surfaces from umbrella sampling simulations. By making a prior assumption of smoothness and taking account of the sampling noise in a consistent fashion, we achieve a significant improvement in accuracy over the state of the art in two or more dimensions or, equivalently, a significant cost reduction to obtain the free energy surface within a prescribed tolerance in both regimes of spatially sparse data and short sampling trajectories. Stemming from its Bayesian interpretation the method provides meaningful error bars without significant additional computation. A software implementation is made available on www.libatoms.org.
△ Less
Submitted 24 July, 2014; v1 submitted 16 December, 2013;
originally announced December 2013.
-
Origin of the structural phase transition in Li7La3Zr2O12
Authors:
N. Bernstein,
M. D. Johannes,
Khang Hoang
Abstract:
Garnet-type Li7La3Zr2O12 (LLZO) is a solid electrolyte material with a low-conductivity tetragonal and a high-conductivity cubic phase. Using density-functional theory and variable cell shape molecular dynamics simulations, we show that the tetragonal phase stability is dependent on a simultaneous ordering of the Li ions on the Li sublattice and a volume-preserving tetragonal distortion that relie…
▽ More
Garnet-type Li7La3Zr2O12 (LLZO) is a solid electrolyte material with a low-conductivity tetragonal and a high-conductivity cubic phase. Using density-functional theory and variable cell shape molecular dynamics simulations, we show that the tetragonal phase stability is dependent on a simultaneous ordering of the Li ions on the Li sublattice and a volume-preserving tetragonal distortion that relieves internal structural strain. Supervalent doping introduces vacancies into the Li sublattice, increasing the overall entropy and reducing the free energy gain from ordering, eventually stabilizing the cubic phase. We show that the critical temperature for cubic phase stability is lowered as Li vacancy concentration (dopant level) is raised and that an activated hop of Li ions from one crystallographic site to another always accompanies the transition. By identifying the relevant mechanism and critical concentrations for achieving the high conductivity phase, this work shows how targeted synthesis could be used to improve electrolytic performance.
△ Less
Submitted 4 September, 2012;
originally announced September 2012.
-
Incommensurate Transverse Anisotropy Induced by Disorder and Spin-Orbit-Vibron Coupling in Mn12-acetate
Authors:
Kyungwha Park,
Mark R. Pederson,
Tunna Baruah,
Noam Bernstein,
Jens Kortus,
Steven L. Richardson,
Enrique del Barco,
Andrew D. Kent,
Steve Hill,
Naresh S. Dalal
Abstract:
It has been shown within density-functional theory that in Mn$_{12}$-acetate there are effects due to disorder by solvent molecules and a coupling between vibrational and electronic degrees of freedom. We calculate the in-plane principal axes of the second-order anisotropy caused by the second effect and compare them with those of the fourth-order anisotropy due to the first effect. We find that…
▽ More
It has been shown within density-functional theory that in Mn$_{12}$-acetate there are effects due to disorder by solvent molecules and a coupling between vibrational and electronic degrees of freedom. We calculate the in-plane principal axes of the second-order anisotropy caused by the second effect and compare them with those of the fourth-order anisotropy due to the first effect. We find that the two types of the principal axes are not commensurate with each other, which results in a complete quenching of the tunnel-splitting oscillation as a function of an applied transverse field.
△ Less
Submitted 8 September, 2004;
originally announced September 2004.
-
Tight-binding study of structure and vibrations of amorphous silicon
Authors:
J. L. Feldman,
N. Bernstein,
D. A. Papaconstantopoulos,
M. J. Mehl
Abstract:
We present a tight-binding calculation that, for the first time, accurately describes the structural, vibrational and elastic properties of amorphous silicon. We compute the interatomic force constants and find an unphysical feature of the Stillinger-Weber empirical potential that correlates with a much noted error in the radial distribution function associated with that potential. We also find…
▽ More
We present a tight-binding calculation that, for the first time, accurately describes the structural, vibrational and elastic properties of amorphous silicon. We compute the interatomic force constants and find an unphysical feature of the Stillinger-Weber empirical potential that correlates with a much noted error in the radial distribution function associated with that potential. We also find that the intrinsic first peak of the radial distribution function is asymmetric, contrary to usual assumptions made in the analysis of diffraction data. We use our results for the normal mode frequencies and polarization vectors to obtain the zero-point broadening effect on the radial distribution function, enabling us to directly compare theory and a high resolution x-ray diffraction experiment.
△ Less
Submitted 14 May, 2004;
originally announced May 2004.
-
A First Principles Density-Functional Calculation of the Electronic and Vibrational Structure of the Key Melanin Monomers
Authors:
B. J. Powell,
T. Baruah,
N. Bernstein,
K. Brake,
Ross H. McKenzie,
P. Meredith,
M. R. Pederson
Abstract:
We report first principles density functional calculations for hydroquinone (HQ), indolequinone (IQ) and semiquinone (SQ). These molecules are believed to be the basic building blocks of the eumelanins, a class of bio-macromolecules with important biological functions (including photoprotection) and with potential for certain bioengineering applications. We have used the DeltaSCF (difference of…
▽ More
We report first principles density functional calculations for hydroquinone (HQ), indolequinone (IQ) and semiquinone (SQ). These molecules are believed to be the basic building blocks of the eumelanins, a class of bio-macromolecules with important biological functions (including photoprotection) and with potential for certain bioengineering applications. We have used the DeltaSCF (difference of self consistent fields) method to study the energy gap between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO), Delta_HL. We show that Delta_HL is similar in IQ and SQ but approximately twice as large in HQ. This may have important implications for our understanding of the observed broad band optical absorption of the eumelanins. The possibility of using this difference in Delta_HL to molecularly engineer the electronic properties of eumelanins is discussed. We calculate the infrared and Raman spectra of the three redox forms from first principles. Each of the molecules have significantly different infrared and Raman signatures, and so these spectra could be used in situ to non-destructively identify the monomeric content of macromolecules. It is hoped that this may be a helpful analytical tool in determining the structure of eumelanin macromolecules and hence in helping to determine the structure-property-function relationships that control the behaviour of the eumelanins.
△ Less
Submitted 22 January, 2004;
originally announced January 2004.
-
Second-order transverse magnetic anisotropy induced by disorders in the single-molecule magnet Mn12
Authors:
Kyungwha Park,
Tunna Baruah,
Noam Bernstein,
Mark R. Pederson
Abstract:
For the single-molecule magnet Mn12, Cornia et al. recently proposed that solvent molecules may cause the quantum tunneling that requires a lower symmetry than S4. However, magnetic quantum tunneling and electron paramagnetic resonance experiments suggested that the proposed theory may not correspond to the measurements. In this regard, we consider positional disorder induced by the solvent mole…
▽ More
For the single-molecule magnet Mn12, Cornia et al. recently proposed that solvent molecules may cause the quantum tunneling that requires a lower symmetry than S4. However, magnetic quantum tunneling and electron paramagnetic resonance experiments suggested that the proposed theory may not correspond to the measurements. In this regard, we consider positional disorder induced by the solvent molecules and orientational disorder by the methyl groups of a Mn12 molecule. We calculate, within density-functional theory, the second-order transverse magnetic anisotropy parameter E and an easy-axis tilting angle induced by the positional disorder and the E value by the orientational disorder. We also calculate the local magnetic anisotropy and the local easy axis for each inequivalent Mn site in different environments to investigate their effects on the global E value. We find that the hydrogen bonding between a Mn12 molecule and the solvent molecule is crucial to obtain a substantial E value and that the E value increases upon geometry relaxations. Our calculations on relaxed geometries show that the largest calculated E value is 0.016 K and that the largest tilting angle is 0.5 degrees. Our largest E value is comparable to experimental results but larger than Cornia et al.'s.
△ Less
Submitted 10 December, 2003;
originally announced December 2003.
-
Electronic, magnetic, and vibrational properties of the molecular magnet Mn4 monomer and dimer
Authors:
Kyungwha Park,
Mark R. Pederson,
Noam Bernstein
Abstract:
A new type of the single-molecule magnet [Mn_4 O_3 Cl_4 (O_2 CEt)_3(py)_3] forms dimers. Recent magnetic hysteresis measurements on this single-molecular magnet revealed interesting phenomena: an absence of quantum tunneling at zero magnetic field and tunneling before magnetic field reversal. This is attributed to a significant antiferromagnetic exchange interaction between different monomers. T…
▽ More
A new type of the single-molecule magnet [Mn_4 O_3 Cl_4 (O_2 CEt)_3(py)_3] forms dimers. Recent magnetic hysteresis measurements on this single-molecular magnet revealed interesting phenomena: an absence of quantum tunneling at zero magnetic field and tunneling before magnetic field reversal. This is attributed to a significant antiferromagnetic exchange interaction between different monomers. To investigate this system, we calculate the electronic structure, magnetic properties, intramolecular and intermolecular exchange interactions using density-functional theory within the generalized-gradient approximation. Our calculations agree with experiment. We also calculate vibrational infrared absorption and Raman scattering intensities for the monomer which can be tested experimentally.
△ Less
Submitted 7 July, 2003;
originally announced July 2003.
-
Magnetic ordering, electronic structure and magnetic anisotropy energy in the high-spin Mn$_{10}$ single molecule magnet
Authors:
Jens Kortus,
Tunna Baruah,
N. Bernstein,
Mark R. Pederson
Abstract:
We report the electronic structure and magnetic ordering of the single molecule magnet [Mn$_{10}$O$_{4}$(2,2'-biphenoxide)$_{4}$Br$_{12}$]$^{4-}$ based on first-principles all-electron density-functional calculations. We find that two of the ten core Mn atoms are coupled antiferromagnetically to the remaining eight, resulting in a ferrimagnetic ground state with total spin S=13. The calculated m…
▽ More
We report the electronic structure and magnetic ordering of the single molecule magnet [Mn$_{10}$O$_{4}$(2,2'-biphenoxide)$_{4}$Br$_{12}$]$^{4-}$ based on first-principles all-electron density-functional calculations. We find that two of the ten core Mn atoms are coupled antiferromagnetically to the remaining eight, resulting in a ferrimagnetic ground state with total spin S=13. The calculated magnetic anisotropy barrier is found to be 9 K in good agreement with experiment. The presence of the Br anions impact the electronic structure and therefore the magnetic properties of the 10 Mn atoms. However, the electric field due to the negative charges has no significant effect on the magnetic anisotropy.
△ Less
Submitted 22 April, 2002;
originally announced April 2002.
-
Vibrationally induced fourth-order magnetic anisotropy and tunnel splittings in Mn_{12}
Authors:
Mark R. Pederson,
Noam Bernstein,
Jens Kortus
Abstract:
From density-functional-theory (DFT) based methods we calculate the vibrational spectrum of the Mn_{12}O_{12}(COOH)_{16}(H_2 O)_4 molecular magnet. Calculated infrared intensities are in accord with experimental studies. There have been no ab initio attempts at determining which interactions account for the fourth-order anisotropy. We show that vibrationally induced distortions of the molecule c…
▽ More
From density-functional-theory (DFT) based methods we calculate the vibrational spectrum of the Mn_{12}O_{12}(COOH)_{16}(H_2 O)_4 molecular magnet. Calculated infrared intensities are in accord with experimental studies. There have been no ab initio attempts at determining which interactions account for the fourth-order anisotropy. We show that vibrationally induced distortions of the molecule contribute to the fourth-order anisotropy Hamiltonian and that the magnitude and sign of the effect (-6.2 K) is in good agreement with all experiments. Vibrationally induced tunnel splittings in isotopically pure and natural samples are predicted.
△ Less
Submitted 18 July, 2002; v1 submitted 19 January, 2002;
originally announced January 2002.
-
The Amorphous-Crystal Interface in Silicon: a Tight-Binding Simulation
Authors:
N. Bernstein,
M. J. Aziz,
E. Kaxiras
Abstract:
The structural features of the interface between the cystalline and amorphous phases of Si solid are studied in simulations based on a combination of empirical interatomic potentials and a nonorthogonal tight-binding model. The tight-binding Hamiltonian was created and tested for the types of structures and distortions anticipated to occur at this interface. The simulations indicate the presence…
▽ More
The structural features of the interface between the cystalline and amorphous phases of Si solid are studied in simulations based on a combination of empirical interatomic potentials and a nonorthogonal tight-binding model. The tight-binding Hamiltonian was created and tested for the types of structures and distortions anticipated to occur at this interface. The simulations indicate the presence of a number of interesting features near the interface. The features that may lead to crystallization upon heating include <110> chains with some defects, most prominently dimers similar to those on the Si(001) 2x1 reconstructed free surface. Within the amorphous region order is lost over very short distances. By examining six different samples with two interfaces each, we find the energy of the amorphous-crystal interface to be 0.49 +/- 0.05 J/m^2
△ Less
Submitted 7 May, 1998;
originally announced May 1998.