-
Extreme Metastability of Diamond and its Transformation to BC8 Post-Diamond Phase of Carbon
Authors:
Kien Nguyen-Cong,
Jonathan T. Willman,
Joseph M. Gonzalez,
Ashley S. Williams,
Anatoly B. Belonoshko,
Stan G. Moore,
Aidan P. Thompson,
Mitchell A. Wood,
Jon H. Eggert,
Marius Millot,
Luis A. Zepeda-Ruiz,
Ivan I. Oleynik
Abstract:
Diamond possesses exceptional physical properties due to its remarkably strong carbon-carbon bonding, leading to significant resilience to structural transformations at very high pressures and temperatures. Despite several experimental attempts, synthesis and recovery of the theoretically predicted post-diamond BC8 phase remains elusive. Through quantum accurate, multi-million atom molecular dynam…
▽ More
Diamond possesses exceptional physical properties due to its remarkably strong carbon-carbon bonding, leading to significant resilience to structural transformations at very high pressures and temperatures. Despite several experimental attempts, synthesis and recovery of the theoretically predicted post-diamond BC8 phase remains elusive. Through quantum accurate, multi-million atom molecular dynamics (MD) simulations, we have uncovered the extreme metastability of diamond at very high pressures, significantly exceeding its range of thermodynamic stability. We predict the post-diamond BC8 phase to be experimentally accessible only within a narrow high pressure-temperature region of the carbon phase diagram. The diamond to BC8 transformation proceeds through pre-melting followed by BC8 nucleation and growth in the metastable carbon liquid. We propose a double-shock compression pathway to achieve BC8 synthesis, which is currently being explored in theory-inspired experiments at the National Ignition Facility.
△ Less
Submitted 22 January, 2024; v1 submitted 16 January, 2024;
originally announced January 2024.
-
Dynamic Formation of Preferentially Lattice Oriented, Self Trapped Hydrogen Clusters
Authors:
M. A. Cusentino,
E. L. Sikorski,
M. J. McCarthy,
A. P. Thompson,
M. A. Wood
Abstract:
A series of MD and DFT simulations were performed to investigate hydrogen self-clustering and retention in tungsten. Using a newly develop machine learned interatomic potential, spontaneous formation of hydrogen platelets was observed after implanting low-energy hydrogen into tungsten at high fluxes and temperatures. The platelets formed along low miller index orientations and neighboring tetrahed…
▽ More
A series of MD and DFT simulations were performed to investigate hydrogen self-clustering and retention in tungsten. Using a newly develop machine learned interatomic potential, spontaneous formation of hydrogen platelets was observed after implanting low-energy hydrogen into tungsten at high fluxes and temperatures. The platelets formed along low miller index orientations and neighboring tetrahedral and octahedral sites and could grow to over 50 atoms in size. High temperatures above 600 K and high hydrogen concentrations were needed to observe significant platelet formation. A critical platelet size of six hydrogen atoms was needed for long term stability. Platelets smaller than this were found to be thermally unstable within a few nanoseconds. To verify these observations, characteristic platelets from the MD simulations were simulated using large-scale DFT. DFT corroborated the MD results in that large platelets were also found to be dynamically stable for five or more hydrogen atoms. The LDOS from the DFT simulated platelets indicated that hydrogen atoms, particularly at the periphery of the platelet, were found to be at least as stable as hydrogen atoms in bulk tungsten. In addition, electrons were found to be localized around hydrogen atoms in the platelet itself and that hydrogen atoms up to 4.2 Angstrom away within the platelet were found to be bonded suggesting that the hydrogen atoms are interacting across longer distances than previously suggested. These results reveal a self-clustering mechanisms for hydrogen within tungsten in the absence of radiation induced or microstructural defects that could be a precursor to blistering and potentially explain the experimentally observed high hydrogen retention particularly in the near surface region.
△ Less
Submitted 19 July, 2023;
originally announced July 2023.
-
The Page curve from the entanglement membrane
Authors:
Mike Blake,
Anthony P. Thompson
Abstract:
We study entanglement dynamics in toy models of black hole information built out of chaotic many-body quantum systems, by utilising a coarse-grained description of entanglement dynamics in such systems known as the `entanglement membrane'. We show that in these models the Page curve associated to the entropy of Hawking radiation arises from a transition in the entanglement membrane around the Page…
▽ More
We study entanglement dynamics in toy models of black hole information built out of chaotic many-body quantum systems, by utilising a coarse-grained description of entanglement dynamics in such systems known as the `entanglement membrane'. We show that in these models the Page curve associated to the entropy of Hawking radiation arises from a transition in the entanglement membrane around the Page time, in an analogous manner to the change in quantum extremal surfaces that leads to the Page curve in semi-classical gravity. We also use the entanglement membrane prescription to study the Hayden-Preskill protocol, and demonstrate how information initially encoded in the black hole is rapidly transferred to the radiation around the Page time. Our results relate recent developments in black hole information to generic features of entanglement dynamics in chaotic many-body quantum systems.
△ Less
Submitted 2 November, 2023; v1 submitted 22 June, 2023;
originally announced June 2023.
-
Exploring Model Complexity in Machine Learned Potentials for Simulated Properties
Authors:
Andrew Rohskopf,
James Goff,
Dionysios Sema,
Kiarash Gordiz,
Ngoc Cuong Nguyen,
Asegun Henry,
Aidan P. Thompson,
Mitchell A. Wood
Abstract:
Machine learning (ML) enables the development of interatomic potentials that promise the accuracy of first principles methods while retaining the low cost and parallel efficiency of empirical potentials. While ML potentials traditionally use atom-centered descriptors as inputs, different models such as linear regression and neural networks can map these descriptors to atomic energies and forces. T…
▽ More
Machine learning (ML) enables the development of interatomic potentials that promise the accuracy of first principles methods while retaining the low cost and parallel efficiency of empirical potentials. While ML potentials traditionally use atom-centered descriptors as inputs, different models such as linear regression and neural networks can map these descriptors to atomic energies and forces. This begs the question: what is the improvement in accuracy due to model complexity irrespective of choice of descriptors? We curate three datasets to investigate this question in terms of ab initio energy and force errors: (1) solid and liquid silicon, (2) gallium nitride, and (3) the superionic conductor LGPS. We further investigate how these errors affect simulated properties with these models and verify if the improvement in fitting errors corresponds to measurable improvement in property prediction. Since linear and nonlinear regression models have different advantages and disadvantages, the results presented herein help researchers choose models for their particular application. By assessing different models, we observe correlations between fitting quantity (e.g. atomic force) error and simulated property error with respect to ab initio values. Such observations can be repeated by other researchers to determine the level of accuracy, and hence model complexity, needed for their particular systems of interest.
△ Less
Submitted 4 June, 2023;
originally announced June 2023.
-
Atomic Representations of Local and Global Chemistry in Complex Alloys
Authors:
Megan J. McCarthy,
Jacob Startt,
Rémi Dingreville,
Aidan P. Thompson,
Mitchell A. Wood
Abstract:
The exceptional properties observed in complex concentrated alloys (CCAs) arise from the interplay between crystalline order and chemical disorder at the atomic scale, complicating a unique determination of properties. In contrast to conventional alloys, CCA properties emerge as distributions due to varying local chemical environments and the specific scale of measurement. Currently there are few…
▽ More
The exceptional properties observed in complex concentrated alloys (CCAs) arise from the interplay between crystalline order and chemical disorder at the atomic scale, complicating a unique determination of properties. In contrast to conventional alloys, CCA properties emerge as distributions due to varying local chemical environments and the specific scale of measurement. Currently there are few ways to quantitatively define, track, and compare local alloy compositions (versus a global label, i.e. equiatomic) contained in a CCA. Molecular dynamics is used here to build descriptive metrics that connect a global alloy composition to the diverse local alloy compositions that define it. A machine-learned interatomic potential for MoNbTaTi is developed and we use these metrics to investigate how property distributions change with excursions in global-local composition space. Short-range order is examined through the lens of local chemistry for the equiatomic composition, demonstrating stark changes in vacancy formation energy with local chemistry evolution.
△ Less
Submitted 18 April, 2023; v1 submitted 7 March, 2023;
originally announced March 2023.
-
Exact results for sheared polar active suspensions with variable liquid crystalline order
Authors:
Aurore Loisy,
Anthony P. Thompson,
Jens Eggers,
Tanniemola B. Liverpool
Abstract:
We consider a confined sheared active polar liquid crystal with a uniform orientation and study the effect of variations in the magnitude of polarization. Restricting our analysis to one-dimensional geometries, we demonstrate that with asymmetric boundary conditions, this system is characterized, macroscopically, by a linear shear stress vs. shear strain relationship that does not pass through the…
▽ More
We consider a confined sheared active polar liquid crystal with a uniform orientation and study the effect of variations in the magnitude of polarization. Restricting our analysis to one-dimensional geometries, we demonstrate that with asymmetric boundary conditions, this system is characterized, macroscopically, by a linear shear stress vs. shear strain relationship that does not pass through the origin: At a zero strain rate, the fluid sustains a non-zero stress. Analytic solutions for the polarization, density, and velocity fields are derived for asymptotically large or small systems and are shown by comparison with precise numerical solutions to be good approximations for finite-size systems.
△ Less
Submitted 11 December, 2022;
originally announced December 2022.
-
Machine Learned Interatomic Potential for Dispersion Strengthened Plasma Facing Components
Authors:
E. L. Sikorski,
M. A. Cusentino,
M. J. McCarthy,
J. Tranchida,
M. A. Wood,
A. P. Thompson
Abstract:
Tungsten (W) is a material of choice for the divertor material due to its high melting temperature, thermal conductivity, and sputtering threshold. However, W has a very high brittle-to-ductile transition temperature and at fusion reactor temperatures ($\geq$1000K) may undergo recrystallization and grain growth. Dispersion-strengthening W with zirconium carbide (ZrC) can improve ductility and limi…
▽ More
Tungsten (W) is a material of choice for the divertor material due to its high melting temperature, thermal conductivity, and sputtering threshold. However, W has a very high brittle-to-ductile transition temperature and at fusion reactor temperatures ($\geq$1000K) may undergo recrystallization and grain growth. Dispersion-strengthening W with zirconium carbide (ZrC) can improve ductility and limit grain growth, but much of the effects of the dispersoids on microstructural evolution and thermomechanical properties at high temperature are still unknown. We present a machine learned Spectral Neighbor Analysis Potential (SNAP) for W-ZrC that can now be used to study these materials. In order to construct a potential suitable for large-scale atomistic simulations at fusion reactor temperatures, it is necessary to train on ab initio data generated for a diverse set of structures, chemical environments, and temperatures. Further accuracy and stability tests of the potential were achieved using objective functions for both material properties and high temperature stability. Validation of lattice parameters, surface energies, bulk moduli, and thermal expansion is confirmed on the optimized potential. Tensile tests of W/ZrC bicrystals show that while the W(110)-ZrC(111) C-terminated bicrystal has the highest ultimate tensile strength (UTS) at room temperature, observed strength decreases with increasing temperature. At 2500K, the terminating C layer diffuses into the W, resulting in a weaker W-Zr interface. Meanwhile, the W(110)-ZrC(111) Zr-terminated bicrystal has the highest UTS at 2500K.
△ Less
Submitted 2 December, 2022;
originally announced December 2022.
-
Chaos and pole-skipping in a simply spinning plasma
Authors:
Markus A. G. Amano,
Mike Blake,
Casey Cartwright,
Matthias Kaminski,
Anthony P. Thompson
Abstract:
We study the relationship between many-body quantum chaos and energy dynamics in holographic quantum field theory states dual to the simply-spinning Myers-Perry-AdS$_5$ black hole. The enhanced symmetry of such black holes allows us to provide a thorough examination of the phenomenon of pole-skipping, that is significantly simpler than a previous analysis of quantum field theory states dual to the…
▽ More
We study the relationship between many-body quantum chaos and energy dynamics in holographic quantum field theory states dual to the simply-spinning Myers-Perry-AdS$_5$ black hole. The enhanced symmetry of such black holes allows us to provide a thorough examination of the phenomenon of pole-skipping, that is significantly simpler than a previous analysis of quantum field theory states dual to the Kerr-AdS$_4$ solution. In particular we give a general proof of pole-skipping in the retarded energy density Green's function of the dual quantum field theory whenever the spatial profile of energy fluctuations satisfies the shockwave equation governing the form of the OTOC. Furthermore, in the large black hole limit we are able to obtain a simple analytic expression for the OTOC for operator configurations on Hopf circles, and demonstrate that the associated Lyapunov exponent and butterfly velocity are robustly related to the locations of a family of pole-skipping points in the energy response. Finally, we note that in contrast to previous studies, our results are valid for any value of rotation and we are able to numerically demonstrate that the dispersion relations of sound modes in the energy response explicitly pass through our pole-skipping locations.
△ Less
Submitted 31 October, 2022;
originally announced November 2022.
-
Predicting electronic structures at any length scale with machine learning
Authors:
Lenz Fiedler,
Normand A. Modine,
Steve Schmerler,
Dayton J. Vogel,
Gabriel A. Popoola,
Aidan P. Thompson,
Sivasankaran Rajamanickam,
Attila Cangi
Abstract:
The properties of electrons in matter are of fundamental importance. They give rise to virtually all molecular and material properties and determine the physics at play in objects ranging from semiconductor devices to the interior of giant gas planets. Modeling and simulation of such diverse applications rely primarily on density functional theory (DFT), which has become the principal method for p…
▽ More
The properties of electrons in matter are of fundamental importance. They give rise to virtually all molecular and material properties and determine the physics at play in objects ranging from semiconductor devices to the interior of giant gas planets. Modeling and simulation of such diverse applications rely primarily on density functional theory (DFT), which has become the principal method for predicting the electronic structure of matter. While DFT calculations have proven to be very useful to the point of being recognized with a Nobel prize in 1998, their computational scaling limits them to small systems. We have developed a machine learning framework for predicting the electronic structure on any length scale. It shows up to three orders of magnitude speedup on systems where DFT is tractable and, more importantly, enables predictions on scales where DFT calculations are infeasible. Our work demonstrates how machine learning circumvents a long-standing computational bottleneck and advances science to frontiers intractable with any current solutions. This unprecedented modeling capability opens up an inexhaustible range of applications in astrophysics, novel materials discovery, and energy solutions for a sustainable future.
△ Less
Submitted 8 December, 2022; v1 submitted 20 October, 2022;
originally announced October 2022.
-
Permutation-adapted complete and independent basis for atomic cluster expansion descriptors
Authors:
James M. Goff,
Charles Sievers,
Mitchell A. Wood,
Aidan P. Thompson
Abstract:
Atomic cluster expansion (ACE) methods provide a systematic way to describe particle local environments of arbitrary body order. For practical applications it is often required that the basis of cluster functions be symmetrized with respect to rotations and permutations. Existing methodologies yield sets of symmetrized functions that are over-complete. These methodologies thus require an additiona…
▽ More
Atomic cluster expansion (ACE) methods provide a systematic way to describe particle local environments of arbitrary body order. For practical applications it is often required that the basis of cluster functions be symmetrized with respect to rotations and permutations. Existing methodologies yield sets of symmetrized functions that are over-complete. These methodologies thus require an additional numerical procedure, such as singular value decomposition (SVD), to eliminate redundant functions. In this work, it is shown that analytical linear relationships for subsets of cluster functions may be derived using recursion and permutation properties of generalized Wigner symbols. From these relationships, subsets (blocks) of cluster functions can be selected such that, within each block, functions are guaranteed to be linearly independent. It is conjectured that this block-wise independent set of permutation-adapted rotation and permutation invariant (PA-RPI) functions forms a complete, independent basis for ACE. Along with the first analytical proofs of block-wise linear dependence of ACE cluster functions and other theoretical arguments, numerical results are offered to demonstrate this. The utility of the method is demonstrated in the development of an ACE interatomic potential for tantalum. Using the new basis functions in combination with Bayesian compressive sensing sparse regression, some high degree descriptors are observed to persist and help achieve high-accuracy models.
△ Less
Submitted 26 February, 2024; v1 submitted 2 August, 2022;
originally announced August 2022.
-
Machine Learning Interatomic Potential for Simulations of Carbon at Extreme Conditions
Authors:
Jonathan T. Willman,
Kien Nguyen-Cong,
Ashley S. Williams,
Anatoly B. Belonoshko,
Stan G. Moore,
Aidan P. Thompson,
Mitchell A. Wood,
Ivan I. Oleynik
Abstract:
A Spectral Neighbor Analysis (SNAP) machine learning interatomic potential (MLIP) has been developed for simulations of carbon at extreme pressures (up to 5 TPa) and temperatures (up to 20,000 K). This was achieved using a large database of experimentally relevant quantum molecular dynamics (QMD) data, training the SNAP potential using a robust machine learning methodology, and performing extensiv…
▽ More
A Spectral Neighbor Analysis (SNAP) machine learning interatomic potential (MLIP) has been developed for simulations of carbon at extreme pressures (up to 5 TPa) and temperatures (up to 20,000 K). This was achieved using a large database of experimentally relevant quantum molecular dynamics (QMD) data, training the SNAP potential using a robust machine learning methodology, and performing extensive validation against QMD and experimental data. The resultant carbon MLIP demonstrates unprecedented accuracy and transferability in predicting the carbon phase diagram, melting curves of crystalline phases, and the shock Hugoniot, all within 3% of QMD. By achieving quantum accuracy and efficient implementation on leadership class high performance computing systems, SNAP advances frontiers of classical MD simulations by enabling atomic-scale insights at experimental time and length scales.
△ Less
Submitted 2 May, 2022;
originally announced May 2022.
-
Training Data Selection for Accuracy and Transferability of Interatomic Potentials
Authors:
David Montes de Oca Zapiain,
Mitchell A. Wood,
Nicholas Lubbers,
Carlos Z. Pereyra,
Aidan P. Thompson,
Danny Perez
Abstract:
Advances in machine learning (ML) techniques have enabled the development of interatomic potentials that promise both the accuracy of first principles methods and the low-cost, linear scaling, and parallel efficiency of empirical potentials. Despite rapid progress in the last few years, ML-based potentials often struggle to achieve transferability, that is, to provide consistent accuracy across co…
▽ More
Advances in machine learning (ML) techniques have enabled the development of interatomic potentials that promise both the accuracy of first principles methods and the low-cost, linear scaling, and parallel efficiency of empirical potentials. Despite rapid progress in the last few years, ML-based potentials often struggle to achieve transferability, that is, to provide consistent accuracy across configurations that significantly differ from those used to train the model. In order to truly realize the promise of ML-based interatomic potentials, it is therefore imperative to develop systematic and scalable approaches for the generation of diverse training sets that ensure broad coverage of the space of atomic environments. This work explores a diverse-by-construction approach that leverages the optimization of the entropy of atomic descriptors to create a very large ($>2\cdot10^{5}$ configurations, $>7\cdot10^{6}$ atomic environments) training set for tungsten in an automated manner, i.e., without any human intervention. This dataset is used to train polynomial as well as multiple neural network potentials with different architectures. For comparison, a corresponding family of potentials were also trained on an expert-curated dataset for tungsten. The models trained to entropy-optimized data exhibited vastly superior transferability compared to the expert-curated models. Furthermore, while the models trained with heavy user input (i.e., domain expertise) yield the lowest errors when tested on similar configurations, out-sample predictions are dramatically more robust when the models are trained on a deliberately diverse set of training data. Herein we demonstrate the development of both accurate and transferable ML potentials using automated and data-driven approaches for generating large and diverse training sets.
△ Less
Submitted 5 May, 2022; v1 submitted 24 January, 2022;
originally announced January 2022.
-
Long-term stellar activity variations and their effect on radial-velocity measurements
Authors:
Jean C. Costes,
Christopher A. Watson,
Ernst de Mooij,
Steven H. Saar,
Xavier Dumusque,
Andrew Collier Cameron,
David F. Phillips,
Maximilian N. Günther,
James S. Jenkins,
Annelies Mortier,
Andrew P. G. Thompson
Abstract:
Long-term stellar activity variations can affect the detectability of long-period and Earth-analogue extrasolar planets. We have, for 54 stars, analysed the long-term trend of five activity indicators: log$R'_\mathrm{HK}$, the cross-correlation function (CCF) bisector span, CCF full-width-at-half-maximum, CCF contrast, and the area of the Gaussian fit to the CCF; and studied their correlation with…
▽ More
Long-term stellar activity variations can affect the detectability of long-period and Earth-analogue extrasolar planets. We have, for 54 stars, analysed the long-term trend of five activity indicators: log$R'_\mathrm{HK}$, the cross-correlation function (CCF) bisector span, CCF full-width-at-half-maximum, CCF contrast, and the area of the Gaussian fit to the CCF; and studied their correlation with the RVs. The sign of the correlations appears to vary as a function of stellar spectral type, and the transition in sign signals a noteworthy change in the stellar activity properties where earlier type stars appear more plage dominated. These transitions become more clearly defined when considered as a function of the convective zone depth. Therefore, it is the convective zone depth (which can be altered by stellar metallicity) that appears to be the underlying fundamental parameter driving the observed activity correlations. In addition, for most of the stars, we find that the RVs become increasingly red-shifted as activity levels increase, which can be explained by the increase in the suppression of convective blue-shift. However, we also find a minority of stars where the RVs become increasingly blue-shifted as activity levels increase. Finally, using the correlation found between activity indicators and RVs, we removed RV signals generated by long-term changes in stellar activity. We find that performing simple cleaning of such long-term signals enables improved planet detection at longer orbital periods.
△ Less
Submitted 5 May, 2021;
originally announced May 2021.
-
Quantum-accurate magneto-elastic predictions with classical spin-lattice dynamics
Authors:
Svetoslav Nikolov,
Mitchell A. Wood,
Attila Cangi,
Jean-Bernard Maillet,
Mihai-Cosmin Marinica,
Aidan P. Thompson,
Michael P. Desjarlais,
Julien Tranchida
Abstract:
A data-driven framework is presented for building magneto-elastic machine-learning interatomic potentials (ML-IAPs) for large-scale spin-lattice dynamics simulations. The magneto-elastic ML-IAPs are constructed by coupling a collective atomic spin model with an ML-IAP. Together they represent a potential energy surface from which the mechanical forces on the atoms and the precession dynamics of th…
▽ More
A data-driven framework is presented for building magneto-elastic machine-learning interatomic potentials (ML-IAPs) for large-scale spin-lattice dynamics simulations. The magneto-elastic ML-IAPs are constructed by coupling a collective atomic spin model with an ML-IAP. Together they represent a potential energy surface from which the mechanical forces on the atoms and the precession dynamics of the atomic spins are computed. Both the atomic spin model and the ML-IAP are parametrized on data from first-principles calculations. We demonstrate the efficacy of our data-driven framework across magneto-structural phase transitions by generating a magneto-elastic ML-IAP for α-iron. The combined potential energy surface yields excellent agreement with first-principles magneto-elastic calculations and quantitative predictions of diverse materials properties including bulk modulus, magnetization, and specific heat across the ferromagnetic-paramagnetic phase transition.
△ Less
Submitted 2 February, 2021; v1 submitted 18 January, 2021;
originally announced January 2021.
-
Rapid Exploration of Optimization Strategies on Advanced Architectures using TestSNAP and LAMMPS
Authors:
Rahulkumar Gayatri,
Stan Moore,
Evan Weinberg,
Nicholas Lubbers,
Sarah Anderson,
Jack Deslippe,
Danny Perez,
Aidan P. Thompson
Abstract:
The exascale race is at an end with the announcement of the Aurora and Frontier machines. This next generation of supercomputers utilize diverse hardware architectures to achieve their compute performance, providing an added onus on the performance portability of applications. An expanding fragmentation of programming models would provide a compounding optimization challenge were it not for the ev…
▽ More
The exascale race is at an end with the announcement of the Aurora and Frontier machines. This next generation of supercomputers utilize diverse hardware architectures to achieve their compute performance, providing an added onus on the performance portability of applications. An expanding fragmentation of programming models would provide a compounding optimization challenge were it not for the evolution of performance-portable frameworks, providing unified models for mapping abstract hierarchies of parallelism to diverse architectures. A solution to this challenge is the evolution of performance-portable frameworks, providing unified models for mapping abstract hierarchies of parallelism to diverse architectures. Kokkos is one such performance portable programming model for C++ applications, providing back-end implementations for each major HPC platform. Even with a performance portable framework, restructuring algorithms to expose higher degrees of parallelism is non-trivial. The Spectral Neighbor Analysis Potential (SNAP) is a machine-learned inter-atomic potential utilized in cutting-edge molecular dynamics simulations. Previous implementations of the SNAP calculation showed a downward trend in their performance relative to peak on newer-generation CPUs and low performance on GPUs. In this paper we describe the restructuring and optimization of SNAP as implemented in the Kokkos CUDA backend of the LAMMPS molecular dynamics package, benchmarked on NVIDIA GPUs. We identify novel patterns of hierarchical parallelism, facilitating a minimization of memory access overheads and pushing the implementation into a compute-saturated regime. Our implementation via Kokkos enables recompile-and-run efficiency on upcoming architectures. We find a $\sim$22x time-to-solution improvement relative to an existing implementation as measured on an NVIDIA Tesla V100-16GB for an important benchmark.
△ Less
Submitted 25 November, 2020;
originally announced November 2020.
-
Accelerating Finite-temperature Kohn-Sham Density Functional Theory with Deep Neural Networks
Authors:
J. Austin Ellis,
Lenz Fiedler,
Gabriel A. Popoola,
Normand A. Modine,
J. Adam Stephens,
Aidan P. Thompson,
Attila Cangi,
Sivasankaran Rajamanickam
Abstract:
We present a numerical modeling workflow based on machine learning (ML) which reproduces the the total energies produced by Kohn-Sham density functional theory (DFT) at finite electronic temperature to within chemical accuracy at negligible computational cost. Based on deep neural networks, our workflow yields the local density of states (LDOS) for a given atomic configuration. From the LDOS, spat…
▽ More
We present a numerical modeling workflow based on machine learning (ML) which reproduces the the total energies produced by Kohn-Sham density functional theory (DFT) at finite electronic temperature to within chemical accuracy at negligible computational cost. Based on deep neural networks, our workflow yields the local density of states (LDOS) for a given atomic configuration. From the LDOS, spatially-resolved, energy-resolved, and integrated quantities can be calculated, including the DFT total free energy, which serves as the Born-Oppenheimer potential energy surface for the atoms. We demonstrate the efficacy of this approach for both solid and liquid metals and compare results between independent and unified machine-learning models for solid and liquid aluminum. Our machine-learning density functional theory framework opens up the path towards multiscale materials modeling for matter under ambient and extreme conditions at a computational scale and cost that is unattainable with current algorithms.
△ Less
Submitted 9 July, 2021; v1 submitted 10 October, 2020;
originally announced October 2020.
-
Simple and efficient algorithms for training machine learning potentials to force data
Authors:
Justin S. Smith,
Nicholas Lubbers,
Aidan P. Thompson,
Kipton Barros
Abstract:
Abstract Machine learning models, trained on data from ab initio quantum simulations, are yielding molecular dynamics potentials with unprecedented accuracy. One limiting factor is the quantity of available training data, which can be expensive to obtain. A quantum simulation often provides all atomic forces, in addition to the total energy of the system. These forces provide much more information…
▽ More
Abstract Machine learning models, trained on data from ab initio quantum simulations, are yielding molecular dynamics potentials with unprecedented accuracy. One limiting factor is the quantity of available training data, which can be expensive to obtain. A quantum simulation often provides all atomic forces, in addition to the total energy of the system. These forces provide much more information than the energy alone. It may appear that training a model to this large quantity of force data would introduce significant computational costs. Actually, training to all available force data should only be a few times more expensive than training to energies alone. Here, we present a new algorithm for efficient force training, and benchmark its accuracy by training to forces from real-world datasets for organic chemistry and bulk aluminum.
△ Less
Submitted 9 June, 2020;
originally announced June 2020.
-
Multi-fidelity machine-learning with uncertainty quantification and Bayesian optimization for materials design: Application to ternary random alloys
Authors:
Anh Tran,
Julien Tranchida,
Tim Wildey,
Aidan P. Thompson
Abstract:
We present a scale-bridging approach based on a multi-fidelity (MF) machine-learning (ML) framework leveraging Gaussian processes (GP) to fuse atomistic computational model predictions across multiple levels of fidelity. Through the posterior variance of the MFGP, our framework naturally enables uncertainty quantification, providing estimates of confidence in the predictions. We used Density Funct…
▽ More
We present a scale-bridging approach based on a multi-fidelity (MF) machine-learning (ML) framework leveraging Gaussian processes (GP) to fuse atomistic computational model predictions across multiple levels of fidelity. Through the posterior variance of the MFGP, our framework naturally enables uncertainty quantification, providing estimates of confidence in the predictions. We used Density Functional Theory as high-fidelity prediction, while a ML interatomic potential is used as the low-fidelity prediction. Practical materials design efficiency is demonstrated by reproducing the ternary composition dependence of a quantity of interest (bulk modulus) across the full aluminum-niobium-titanium ternary random alloy composition space. The MFGP is then coupled to a Bayesian optimization procedure and the computational efficiency of this approach is demonstrated by performing an on-the-fly search for the global optimum of bulk modulus in the ternary composition space. The framework presented in this manuscript is the first application of MFGP to atomistic materials simulations fusing predictions between Density Functional Theory and classical interatomic potential calculations.
△ Less
Submitted 5 August, 2020; v1 submitted 29 May, 2020;
originally announced June 2020.
-
The spectral impact of magnetic activity on disk-integrated HARPS-N solar observations: exploring new activity indicators
Authors:
A. P. G. Thompson,
C. A. Watson,
R. D. Haywood,
J. C. Costes,
E. de Mooij,
A. Collier Cameron,
X. Dumusque,
D. F. Phillips,
S. H. Saar,
A. Mortier,
T. W. Milbourne,
S. Aigrain,
H. M. Cegla,
D. Charbonneau,
R. Cosentino,
A. Ghedina,
D. W. Latham,
M. López-Morales,
G. Micela,
E. Molinari,
E. Poretti,
A. Sozzetti,
S. Thompson,
R. Walsworth
Abstract:
Stellar activity is the major roadblock on the path to finding true Earth-analogue planets with the Doppler technique. Thus, identifying new indicators that better trace magnetic activity (i.e. faculae and spots) is crucial to aid in disentangling these signals from that of a planet's Doppler wobble. In this work, we investigate activity related features as seen in disk-integrated spectra from the…
▽ More
Stellar activity is the major roadblock on the path to finding true Earth-analogue planets with the Doppler technique. Thus, identifying new indicators that better trace magnetic activity (i.e. faculae and spots) is crucial to aid in disentangling these signals from that of a planet's Doppler wobble. In this work, we investigate activity related features as seen in disk-integrated spectra from the HARPS-N solar telescope. We divide high-activity spectral echelle orders by low-activity master templates (as defined using both log R'HK and images from the Solar Dynamics Observatory, SDO), creating "relative spectra". With resolved images of the surface of the Sun (via SDO), the faculae and spot filling factors can be calculated, giving a measure of activity independent of, and in addition to, log R'HK. We find pseudo-emission (and pseudo-absorption) features in the relative spectra that are similar to those reported in our previous work on alpha Cen B. In alpha Cen B, the features are shown to correlate better to changes in faculae filling factor than spot filling factor. In this work we more confidently identify changes in faculae coverage of the visible hemisphere of the Sun as the source of features produced in the relative spectra. Finally, we produce trailed spectra to observe the RV component of the features, which show that the features move in a redward direction as one would expect when tracking active regions rotating on the surface of a star.
△ Less
Submitted 21 April, 2020;
originally announced April 2020.
-
Explicit Multi-element Extension of the Spectral Neighbor Analysis Potential for Chemically Complex Systems
Authors:
Mary Alice Cusentino,
Mitchell A. Wood,
Aidan P. Thompson
Abstract:
A natural extension of the descriptors used in the Spectral Neighbor Analysis Potential (SNAP) method is derived to treat atomic interactions in chemically complex systems. Atomic environment descriptors within SNAP are obtained from a basis function expansion of the weighted density of neighboring atoms. This new formulation instead partitions the neighbor density into partial densities for each…
▽ More
A natural extension of the descriptors used in the Spectral Neighbor Analysis Potential (SNAP) method is derived to treat atomic interactions in chemically complex systems. Atomic environment descriptors within SNAP are obtained from a basis function expansion of the weighted density of neighboring atoms. This new formulation instead partitions the neighbor density into partial densities for each chemical element, thus leading to explicit multi-element descriptors. For $N_{elem}$ chemical elements, the number of descriptors increases as $\mathcal{O}(N_{elem}^3)$, while the computational cost of the force calculation as implemented in LAMMPS is limited to $\mathcal{O}(N_{elem}^2)$ and the favorable linear scaling in the number of atoms is retained. We demonstrate these chemically aware descriptors by producing an interatomic potential for indium phosphide capable of capturing high-energy defects that result from radiation damage cascades. This new explicit multi-element SNAP method reproduces the relaxed defect formation energies with substantially greater accuracy than weighted-density SNAP, while retaining accurate representation of the bulk indium phosphide properties.
△ Less
Submitted 25 March, 2020;
originally announced March 2020.
-
NGTS-8b and NGTS-9b: two non-inflated hot-Jupiters
Authors:
Jean C. Costes,
Christopher A. Watson,
Claudia Belardi,
Ian P. Braker,
Matthew R. Burleigh,
Sarah L. Casewell,
Philipp Eigmüller,
Maximilian N. Günther,
James A. G. Jackman,
Louise D. Nielsen,
Maritza G. Soto,
Oliver Turner,
David R. Anderson,
Daniel Bayliss,
François Bouchy,
Joshua T. Briegal,
Edward M. Bryant,
Juan Cabrera,
Alexander Chaushev,
Szilard Csizmadia,
Anders Erikson,
Samuel Gill,
Edward Gillen,
Michael R. Goad,
Matthew J. Hooton
, et al. (13 additional authors not shown)
Abstract:
We report the discovery, by the Next Generation Transit Survey (NGTS), of two hot-Jupiters NGTS-8b and NGTS-9b. These orbit a V = 13.68 K0V star (Teff = 5241 +/- 50 K) with a period of 2.49970 days, and a V = 12.80 F8V star (Teff = 6330 +/- 130 K) in 4.43527 days, respectively. The transits were independently verified by follow-up photometric observations with the SAAO 1.0-m and Euler telescopes,…
▽ More
We report the discovery, by the Next Generation Transit Survey (NGTS), of two hot-Jupiters NGTS-8b and NGTS-9b. These orbit a V = 13.68 K0V star (Teff = 5241 +/- 50 K) with a period of 2.49970 days, and a V = 12.80 F8V star (Teff = 6330 +/- 130 K) in 4.43527 days, respectively. The transits were independently verified by follow-up photometric observations with the SAAO 1.0-m and Euler telescopes, and we report on the planetary parameters using HARPS, FEROS and CORALIE radial velocities. NGTS-8b has a mass, 0.93 +0.04 -0.03 MJ and a radius, 1.09 +/- 0.03 RJ similar to Jupiter, resulting in a density of 0.89 +0.08 -0.07 g cm-3. This is in contrast to NGTS-9b, which has a mass of 2.90 +/- 0.17 MJ and a radius of 1.07 +/- 0.06 RJ , resulting in a much greater density of 2.93 +0.53 -0.49 g cm-3. Statistically, the planetary parameters put both objects in the regime where they would be expected to exhibit larger than predicted radii. However, we find that their radii are in agreement with predictions by theoretical non-inflated models.
△ Less
Submitted 7 November, 2019;
originally announced November 2019.
-
A Performance and Cost Assessment of Machine Learning Interatomic Potentials
Authors:
Yunxing Zuo,
Chi Chen,
Xiangguo Li,
Zhi Deng,
Yiming Chen,
Jörg Behler,
Gábor Csányi,
Alexander V. Shapeev,
Aidan P. Thompson,
Mitchell A. Wood,
Shyue Ping Ong
Abstract:
Machine learning of the quantitative relationship between local environment descriptors and the potential energy surface of a system of atoms has emerged as a new frontier in the development of interatomic potentials (IAPs). Here, we present a comprehensive evaluation of ML-IAPs based on four local environment descriptors --- Behler-Parrinello symmetry functions, smooth overlap of atomic positions…
▽ More
Machine learning of the quantitative relationship between local environment descriptors and the potential energy surface of a system of atoms has emerged as a new frontier in the development of interatomic potentials (IAPs). Here, we present a comprehensive evaluation of ML-IAPs based on four local environment descriptors --- Behler-Parrinello symmetry functions, smooth overlap of atomic positions (SOAP), the Spectral Neighbor Analysis Potential (SNAP) bispectrum components, and moment tensors --- using a diverse data set generated using high-throughput density functional theory (DFT) calculations. The data set comprising bcc (Li, Mo) and fcc (Cu, Ni) metals and diamond group IV semiconductors (Si, Ge) is chosen to span a range of crystal structures and bonding. All descriptors studied show excellent performance in predicting energies and forces far surpassing that of classical IAPs, as well as predicting properties such as elastic constants and phonon dispersion curves. We observe a general trade-off between accuracy and the degrees of freedom of each model, and consequently computational cost. We will discuss these trade-offs in the context of model selection for molecular dynamics and other applications.
△ Less
Submitted 24 July, 2019; v1 submitted 20 June, 2019;
originally announced June 2019.
-
Data-driven Material Models for Atomistic Simulation
Authors:
Mitchell A. Wood,
Mary Alice Cusentino,
Brian D. Wirth,
Aidan P. Thompson
Abstract:
The central approximation made in classical molecular dynamics simulation of materials is the interatomic potential used to calculate the forces on the atoms. Great effort and ingenuity is required to construct viable functional forms and find accurate parameterizations for potentials using traditional approaches. Machine-learning has emerged as an effective alternative approach to develop accurat…
▽ More
The central approximation made in classical molecular dynamics simulation of materials is the interatomic potential used to calculate the forces on the atoms. Great effort and ingenuity is required to construct viable functional forms and find accurate parameterizations for potentials using traditional approaches. Machine-learning has emerged as an effective alternative approach to develop accurate and robust interatomic potentials. Starting with a very general model form, the potential is learned directly from a database of electronic structure calculations and therefore can be viewed as a multiscale link between quantum and classical atomistic simulations. Risk of inaccurate extrapolation exists outside the narrow range of time- and length-scales where the two methods can be directly compared. In this work, we use the Spectral Neighbor Analysis Potential (SNAP) and show how a fit can be produced with minimal interpolation errors which is also robust in extrapolating beyond training. To demonstrate the method, we have developed a new tungsten-beryllium potential suitable for the full range of binary compositions. Subsequently, large-scale molecular dynamics simulations were performed of high energy Be atom implantation onto the (001) surface of solid tungsten. The new machine learned W-Be potential generates a population of implantation structures consistent with quantum calculations of defect formation energies. A very shallow (<2nm) average Be implantation depth is predicted which may explain ITER diverter degradation in the presence of beryllium.
△ Less
Submitted 18 March, 2019; v1 submitted 25 February, 2019;
originally announced February 2019.
-
NGTS-2b: An inflated hot-Jupiter transiting a bright F-dwarf
Authors:
Liam Raynard,
Michael R. Goad,
Edward Gillen,
Louise D. Nielsen,
Christopher A. Watson,
Andrew P. G. Thompson,
James McCormac,
Daniel Bayliss,
Maritza Soto,
Szilard Csizmadia,
Alexander Chaushev,
Matthew R. Burleigh,
Richard Alexander,
David J. Armstrong,
François Bouchy,
Joshua T. Briegal,
Juan Cabrera,
Sarah L. Casewell,
Bruno Chazelas,
Benjamin F. Cooke,
Philipp Eigmüller,
Anders Erikson,
Boris T. Gänsicke,
Andrew Grange,
Maximilian N. Günther
, et al. (20 additional authors not shown)
Abstract:
We report the discovery of NGTS-2b, an inflated hot-Jupiter transiting a bright F5V star (2MASS J14202949-3112074; $T_{\rm eff}$=$6478^{+94}_{-89}$ K), discovered as part of the Next Generation Transit Survey (NGTS). The planet is in a P=4.51 day orbit with mass $0.74^{+0.13}_{-0.12}$ M$_{J}$, radius $1.595^{+0.047}_{-0.045}$ R$_{J}$ and density $0.226^{+0.040}_{-0.038}$ g cm$^{-3}$; therefore one…
▽ More
We report the discovery of NGTS-2b, an inflated hot-Jupiter transiting a bright F5V star (2MASS J14202949-3112074; $T_{\rm eff}$=$6478^{+94}_{-89}$ K), discovered as part of the Next Generation Transit Survey (NGTS). The planet is in a P=4.51 day orbit with mass $0.74^{+0.13}_{-0.12}$ M$_{J}$, radius $1.595^{+0.047}_{-0.045}$ R$_{J}$ and density $0.226^{+0.040}_{-0.038}$ g cm$^{-3}$; therefore one of the lowest density exoplanets currently known. With a relatively deep 1.0% transit around a bright V=10.96 host star, NGTS-2b is a prime target for probing giant planet composition via atmospheric transmission spectroscopy. The rapid rotation ($vsin$i=$15.2\pm0.8$ km s$^{-1}$) also makes this system an excellent candidate for Rossiter-McLaughlin follow-up observations, to measure the sky-projected stellar obliquity. NGTS-2b was confirmed without the need for follow-up photometry, due to the high precision of the NGTS photometry.
△ Less
Submitted 20 September, 2018; v1 submitted 26 May, 2018;
originally announced May 2018.
-
Ground-based detection of G star superflares with NGTS
Authors:
James A. G. Jackman,
Peter J. Wheatley,
Chloe E. Pugh,
Boris T. Gänsicke,
Edward Gillen,
Anne-Marie Broomhall,
David J. Armstrong,
Matthew R. Burleigh,
Alexander Chaushev,
Philipp Eigmüller,
Anders Erikson,
Michael R. Goad,
Andrew Grange,
Maximilian N. Günther,
James S. Jenkins,
James McCormac,
Liam Raynard,
Andrew P. G. Thompson,
Stéphane Udry,
Simon Walker,
Christopher A. Watson,
Richard G. West
Abstract:
We present high cadence detections of two superflares from a bright G8 star (V = 11.56) with the Next Generation Transit Survey (NGTS). We improve upon previous superflare detections by resolving the flare rise and peak, allowing us to fit a solar flare inspired model without the need for arbitrary break points between rise and decay. Our data also enables us to identify substructure in the flares…
▽ More
We present high cadence detections of two superflares from a bright G8 star (V = 11.56) with the Next Generation Transit Survey (NGTS). We improve upon previous superflare detections by resolving the flare rise and peak, allowing us to fit a solar flare inspired model without the need for arbitrary break points between rise and decay. Our data also enables us to identify substructure in the flares. From changing starspot modulation in the NGTS data we detect a stellar rotation period of 59 hours, along with evidence for differential rotation. We combine this rotation period with the observed \textit{ROSAT} X-ray flux to determine that the star's X-ray activity is saturated. We calculate the flare bolometric energies as $5.4^{+0.8}_{-0.7}\times10^{34}$ and $2.6^{+0.4}_{-0.3}\times10^{34}$ erg and compare our detections with G star superflares detected in the \textit{Kepler} survey. We find our main flare to be one of the largest amplitude superflares detected from a bright G star. With energies more than 100 times greater than the Carrington event, our flare detections demonstrate the role that ground-based instruments such as NGTS can have in assessing the habitability of Earth-like exoplanets, particularly in the era of \textit{PLATO}.
△ Less
Submitted 10 April, 2018;
originally announced April 2018.
-
Massively parallel symplectic algorithm for coupled magnetic spin dynamics and molecular dynamics
Authors:
J. Tranchida,
S. J. Plimpton,
P. Thibaudeau,
A. P. Thompson
Abstract:
A parallel implementation of coupled spin-lattice dynamics in the LAMMPS molecular dynamics package is presented. The equations of motion for both spin only and coupled spin-lattice dynamics are first reviewed, including a detailed account of how magneto-mechanical potentials can be used to perform a proper coupling between spin and lattice degrees of freedom. A symplectic numerical integration al…
▽ More
A parallel implementation of coupled spin-lattice dynamics in the LAMMPS molecular dynamics package is presented. The equations of motion for both spin only and coupled spin-lattice dynamics are first reviewed, including a detailed account of how magneto-mechanical potentials can be used to perform a proper coupling between spin and lattice degrees of freedom. A symplectic numerical integration algorithm is then presented which combines the Suzuki-Trotter decomposition for non-commuting variables and conserves the geometric properties of the equations of motion. The numerical accuracy of the serial implementation was assessed by verifying that it conserves the total energy and the norm of the total magnetization up to second order in the timestep size. Finally, a very general parallel algorithm is proposed that allows large spin-lattice systems to be efficiently simulated on large numbers of processors without degrading its mathematical accuracy. Its correctness as well as scaling efficiency were tested for realistic coupled spin-lattice systems, confirming that the new parallel algorithm is both accurate and efficient.
△ Less
Submitted 30 January, 2018;
originally announced January 2018.
-
Extending the Accuracy of the SNAP Interatomic Potential Form
Authors:
Mitchell A. Wood,
Aidan P. Thompson
Abstract:
The Spectral Neighbor Analysis Potential (SNAP) is a classical interatomic potential that expresses the energy of each atom as a linear function of selected bispectrum components of the neighbor atoms. An extension of the SNAP form is proposed that includes quadratic terms in the bispectrum components. The extension is shown to provide a large increase in accuracy relative to the linear form, whil…
▽ More
The Spectral Neighbor Analysis Potential (SNAP) is a classical interatomic potential that expresses the energy of each atom as a linear function of selected bispectrum components of the neighbor atoms. An extension of the SNAP form is proposed that includes quadratic terms in the bispectrum components. The extension is shown to provide a large increase in accuracy relative to the linear form, while incurring only a modest increase in computational cost. The mathematical structure of the quadratic SNAP form is similar to the embedded atom method (EAM), with the SNAP bispectrum components serving as counterparts to the two-body density functions in EAM. The effectiveness of the new form is demonstrated using an extensive set of training data for tantalum structures. Similarly to artificial neural network potentials, the quadratic SNAP form requires substantially more training data in order to prevent overfitting. The quality of this new potential form is measured through a robust cross-validation analysis.
△ Less
Submitted 20 February, 2018; v1 submitted 29 November, 2017;
originally announced November 2017.
-
Multiscale Modeling of Shock Wave Localization in Porous Energetic Material
Authors:
Mitchell A. Wood,
David E. Kittell,
Cole D. Yarrington,
Aidan P. Thompson
Abstract:
Shock wave interactions with defects, such as pores, are known to play a key role in the chemical initiation of energetic materials. The shock response of hexanitrostilbene is studied through a combination of large scale reactive molecular dynamics and mesoscale hydrodynamic simulations. In order to extend our simulation capability at the mesoscale to include weak shock conditions (< 6 GPa), atomi…
▽ More
Shock wave interactions with defects, such as pores, are known to play a key role in the chemical initiation of energetic materials. The shock response of hexanitrostilbene is studied through a combination of large scale reactive molecular dynamics and mesoscale hydrodynamic simulations. In order to extend our simulation capability at the mesoscale to include weak shock conditions (< 6 GPa), atomistic simulations of pore collapse are used to define a strain rate dependent strength model. Comparing these simulation methods allows us to impose physically-reasonable constraints on the mesoscale model parameters. In doing so, we have been able to study shock waves interacting with pores as a function of this viscoplastic material response. We find that the pore collapse behavior of weak shocks is characteristically different to that of strong shocks.
△ Less
Submitted 7 November, 2017;
originally announced November 2017.
-
NGTS-1b: A hot Jupiter transiting an M-dwarf
Authors:
Daniel Bayliss,
Edward Gillen,
Philipp Eigmuller,
James McCormac,
Richard D. Alexander,
David J. Armstrong,
Rachel S. Booth,
Francois Bouchy,
Matthew R. Burleigh,
Juan Cabrera,
Sarah L. Casewell,
Alexander Chaushev,
Bruno Chazelas,
Szilard Csizmadia,
Anders Erikson,
Francesca Faedi,
Emma Foxell,
Boris T. Gansicke,
Michael R. Goad,
Andrew Grange,
Maximilian N. Gunther,
Simon T. Hodgkin,
James Jackman,
James S. Jenkins,
Gregory Lambert
, et al. (18 additional authors not shown)
Abstract:
We present the discovery of NGTS-1b, a hot-Jupiter transiting an early M-dwarf host ($T_{eff}=3916^{+71}_{-63}~K$) in a P=2.674d orbit discovered as part of the Next Generation Transit Survey (NGTS). The planet has a mass of $0.812^{+0.066}_{-0.075}~M_{J}$, making it the most massive planet ever discovered transiting an M-dwarf. The radius of the planet is $1.33^{+0.61}_{-0.33}~R_{J}$. Since the t…
▽ More
We present the discovery of NGTS-1b, a hot-Jupiter transiting an early M-dwarf host ($T_{eff}=3916^{+71}_{-63}~K$) in a P=2.674d orbit discovered as part of the Next Generation Transit Survey (NGTS). The planet has a mass of $0.812^{+0.066}_{-0.075}~M_{J}$, making it the most massive planet ever discovered transiting an M-dwarf. The radius of the planet is $1.33^{+0.61}_{-0.33}~R_{J}$. Since the transit is grazing, we determine this radius by modelling the data and placing a prior on the density from the population of known gas giant planets. NGTS-1b is the third transiting giant planet found around an M-dwarf, reinforcing the notion that close-in gas giants can form and migrate similar to the known population of hot Jupiters around solar type stars. The host star shows no signs of activity, and the kinematics hint at the star being from the thick disk population. With a deep (2.5%) transit around a $K=11.9$ host, NGTS-1b will be a strong candidate to probe giant planet composition around M-dwarfs via JWST transmission spectroscopy.
△ Less
Submitted 30 October, 2017;
originally announced October 2017.
-
Centroid vetting of transiting planet candidates from the Next Generation Transit Survey
Authors:
Maximilian N. Günther,
Didier Queloz,
Edward Gillen,
James McCormac,
Daniel Bayliss,
Francois Bouchy,
Simon. R. Walker,
Richard G. West,
Philipp Eigmüller,
Alexis M. S. Smith,
David J. Armstrong,
Matthew Burleigh,
Sarah L. Casewell,
Alexander P. Chaushev,
Michael R. Goad,
Andrew Grange,
James Jackman,
James S. Jenkins,
Tom Louden,
Maximiliano Moyano,
Don Pollacco,
Katja Poppenhaeger,
Heike Rauer,
Liam Raynard,
Andrew P. G. Thompson
, et al. (3 additional authors not shown)
Abstract:
The Next Generation Transit Survey (NGTS), operating in Paranal since 2016, is a wide-field survey to detect Neptunes and super-Earths transiting bright stars, which are suitable for precise radial velocity follow-up and characterisation. Thereby, its sub-mmag photometric precision and ability to identify false positives are crucial. Particularly, variable background objects blended in the photome…
▽ More
The Next Generation Transit Survey (NGTS), operating in Paranal since 2016, is a wide-field survey to detect Neptunes and super-Earths transiting bright stars, which are suitable for precise radial velocity follow-up and characterisation. Thereby, its sub-mmag photometric precision and ability to identify false positives are crucial. Particularly, variable background objects blended in the photometric aperture frequently mimic Neptune-sized transits and are costly in follow-up time. These objects can best be identified with the centroiding technique: if the photometric flux is lost off-centre during an eclipse, the flux centroid shifts towards the centre of the target star. Although this method has successfully been employed by the Kepler mission, it has previously not been implemented from the ground. We present a fully-automated centroid vetting algorithm developed for NGTS, enabled by our high-precision auto-guiding. Our method allows detecting centroid shifts with an average precision of 0.75 milli-pixel, and down to 0.25 milli-pixel for specific targets, for a pixel size of 4.97 arcsec. The algorithm is now part of the NGTS candidate vetting pipeline and automatically employed for all detected signals. Further, we develop a joint Bayesian fitting model for all photometric and centroid data, allowing to disentangle which object (target or background) is causing the signal, and what its astrophysical parameters are. We demonstrate our method on two NGTS objects of interest. These achievements make NGTS the first ground-based wide-field transit survey ever to successfully apply the centroiding technique for automated candidate vetting, enabling the production of a robust candidate list before follow-up.
△ Less
Submitted 4 September, 2017; v1 submitted 25 July, 2017;
originally announced July 2017.
-
Quantum-Accurate Molecular Dynamics Potential for Tungsten
Authors:
Mitchell A. Wood,
Aidan P. Thompson
Abstract:
The purpose of this short contribution is to report on the development of a Spectral Neighbor Analysis Potential (SNAP) for tungsten. We have focused on the characterization of elastic and defect properties of the pure material in order to support molecular dynamics simulations of plasma-facing materials in fusion reactors. A parallel genetic algorithm approach was used to efficiently search for f…
▽ More
The purpose of this short contribution is to report on the development of a Spectral Neighbor Analysis Potential (SNAP) for tungsten. We have focused on the characterization of elastic and defect properties of the pure material in order to support molecular dynamics simulations of plasma-facing materials in fusion reactors. A parallel genetic algorithm approach was used to efficiently search for fitting parameters optimized against a large number of objective functions. In addition, we have shown that this many-body tungsten potential can be used in conjunction with a simple helium pair potential to produce accurate defect formation energies for the W-He binary system.
△ Less
Submitted 22 February, 2017;
originally announced February 2017.
-
The Changing Face of $α$ Centauri B: Probing plage and stellar activity in K-dwarfs
Authors:
A. P. G. Thompson,
C. A. Watson,
E. J. W. de Mooij,
D. B. Jess
Abstract:
A detailed knowledge of stellar activity is crucial for understanding stellar dynamos, as well as pushing exoplanet radial-velocity detection limits towards Earth analogue confirmation. We directly compare archival HARPS spectra taken at the minimum in $α$ Cen B's activity cycle to a high-activity state when clear rotational modulation of $\log{R'_{HK}}$ is visible. Relative to the inactive spectr…
▽ More
A detailed knowledge of stellar activity is crucial for understanding stellar dynamos, as well as pushing exoplanet radial-velocity detection limits towards Earth analogue confirmation. We directly compare archival HARPS spectra taken at the minimum in $α$ Cen B's activity cycle to a high-activity state when clear rotational modulation of $\log{R'_{HK}}$ is visible. Relative to the inactive spectra, we find a large number of narrow pseudo-emission features in the active spectra with strengths that are rotationally modulated. These features most likely originate from plage, spots, or a combination of both. They also display radial velocity variations of $\sim$300 m s$^{-1}$ - consistent with an active region rotating across the stellar surface. Furthermore, we see evidence that some of the lines originating from the `active immaculate' photosphere appear broader relative to the `inactive immaculate' case. This may be due to enhanced contributions of e.g. magnetic bright points to these lines, which then causes additional line broadening. More detailed analysis may enable measurements of plage and spot coverage using single spectra in the future.
△ Less
Submitted 6 February, 2017;
originally announced February 2017.
-
A Spectral Analysis Method for Automated Generation of Quantum-Accurate Interatomic Potentials
Authors:
Aidan P. Thompson,
Laura P. Swiler,
Christian R. Trott,
Stephen M. Foiles,
Garritt J. Tucker
Abstract:
We present a new interatomic potential for solids and liquids called Spectral Neighbor Analysis Potential (SNAP). The SNAP potential has a very general form and uses machine-learning techniques to reproduce the energies, forces, and stress tensors of a large set of small configurations of atoms, which are obtained using high-accuracy quantum electronic structure (QM) calculations. The local enviro…
▽ More
We present a new interatomic potential for solids and liquids called Spectral Neighbor Analysis Potential (SNAP). The SNAP potential has a very general form and uses machine-learning techniques to reproduce the energies, forces, and stress tensors of a large set of small configurations of atoms, which are obtained using high-accuracy quantum electronic structure (QM) calculations. The local environment of each atom is characterized by a set of bispectrum components of the local neighbor density projected on to a basis of hyperspherical harmonics in four dimensions. The bispectrum components are the same bond-orientational order parameters employed by the GAP potential [arXiv:0910.1019]. The SNAP potential, unlike GAP, assumes a linear relationship between atom energy and bispectrum components. The linear SNAP coefficients are determined using weighted least-squares linear regression against the full QM training set. This allows the SNAP potential to be fit in a robust, automated manner to large QM data sets using many bispectrum coefficients.
△ Less
Submitted 12 September, 2014;
originally announced September 2014.
-
Searching for Globally Optimal Functional Forms for Inter-Atomic Potentials Using Parallel Tempering and Genetic Programming
Authors:
A. Slepoy,
A. P. Thompson,
M. D. Peters
Abstract:
We develop a Genetic Programming-based methodology that enables discovery of novel functional forms for classical inter-atomic force-fields, used in molecular dynamics simulations. Unlike previous efforts in the field, that fit only the parameters to the fixed functional forms, we instead use a novel algorithm to search the space of many possible functional forms. While a follow-on practical pro…
▽ More
We develop a Genetic Programming-based methodology that enables discovery of novel functional forms for classical inter-atomic force-fields, used in molecular dynamics simulations. Unlike previous efforts in the field, that fit only the parameters to the fixed functional forms, we instead use a novel algorithm to search the space of many possible functional forms. While a follow-on practical procedure will use experimental and {\it ab inito} data to find an optimal functional form for a forcefield, we first validate the approach using a manufactured solution. This validation has the advantage of a well-defined metric of success. We manufactured a training set of atomic coordinate data with an associated set of global energies using the well-known Lennard-Jones inter-atomic potential. We performed an automatic functional form fitting procedure starting with a population of random functions, using a genetic programming functional formulation, and a parallel tempering Metropolis-based optimization algorithm. Our massively-parallel method independently discovered the Lennard-Jones function after searching for several hours on 100 processors and covering a miniscule portion of the configuration space. We find that the method is suitable for unsupervised discovery of functional forms for inter-atomic potentials/force-fields. We also find that our parallel tempering Metropolis-based approach significantly improves the optimization convergence time, and takes good advantage of the parallel cluster architecture.
△ Less
Submitted 18 August, 2006;
originally announced August 2006.