Multi-Objective Optimization for Rapid Identification of Novel Compound Metals for Interconnect Applications.

Interconnect materials play the critical role of routing energy and information in integrated circuits. However, established bulk conductors, such as copper, perform poorly when scaled down beyond 10 nm, limiting the scalability of logic devices. Here, a multi-objective search is developed, combined with first-principles calculations, to rapidly screen over 15,000 materials and discover new interconnect candidates. This approach simultaneously optimizes the bulk electronic conductivity, surface scattering time, and chemical stability using physically motivated surrogate properties accessible from materials databases. Promising local interconnects are identified that have the potential to outperform ruthenium, the current state-of-the-art post-Cu material, and also semi-global interconnects with potentially large skin depths at the GHz operation frequency. The approach is validated on one of the identified candidates, CoPt, using both ab initio and experimental transport studies, showcasing its potential to supplant Ru and Cu for future local interconnects.

Structural phase transition in monolayer MoTe2 driven by electrostatic doping.

Monolayers of transition-metal dichalcogenides (TMDs) exhibit numerous crystal phases with distinct structures, symmetries and physical properties. Exploring the physics of transitions between these different structural phases in two dimensions may provide a means of switching material properties, with implications for potential applications. Structural phase transitions in TMDs have so far been induced by thermal or chemical means; purely electrostatic control over crystal phases through electrostatic doping was recently proposed as a theoretical possibility, but has not yet been realized. Here we report the experimental demonstration of an electrostatic-doping-driven phase transition between the hexagonal and monoclinic phases of monolayer molybdenum ditelluride (MoTe2). We find that the phase transition shows a hysteretic loop in Raman spectra, and can be reversed by increasing or decreasing the gate voltage. We also combine second-harmonic generation spectroscopy with polarization-resolved Raman spectroscopy to show that the induced monoclinic phase preserves the crystal orientation of the original hexagonal phase. Moreover, this structural phase transition occurs simultaneously across the whole sample. This electrostatic-doping control of structural phase transition opens up new possibilities for developing phase-change devices based on atomically thin membranes.

Predicting molecule size distribution in hydrocarbon pyrolysis using random graph theory.

Hydrocarbon pyrolysis is a complex process involving large numbers of chemical species and types of chemical reactions. Its quantitative description is important for planetary sciences, in particular, for understanding the processes occurring in the interior of icy planets, such as Uranus and Neptune, where small hydrocarbons are subjected to high temperature and pressure. We propose a computationally cheap methodology based on an originally developed ten-reaction model and the configurational model from random graph theory. This methodology generates accurate predictions for molecule size distributions for a variety of initial chemical compositions and temperatures ranging from 3200 to 5000 K. Specifically, we show that the size distribution of small molecules is particularly well predicted, and the size of the largest molecule can be accurately predicted provided that this molecule is not too large.

Pit viper thermography: the pit organ used by crotaline snakes to detect thermal contrast has poor spatial resolution.

Pit vipers detect infrared radiation by means of temperature contrasts created on their pit organ membranes. Signals from pit organs integrate with visual signals in the optic tectum, leading to the conjecture that the facial pits operate as an extension of the visual system. Because similar mechanisms underlie thermal imaging technology, imagery from thermal cameras is often used to infer how pit vipers perceive their environment. However, pit organs lack a focusing mechanism, and biophysical models predict that pit organs should have poor spatial resolution compared with thermal imaging cameras. Nevertheless, behavioral studies occasionally suggest pits may have better resolution than predicted by biophysical models, indicating that processing in the central nervous system may improve imaging. To estimate the spatial resolution of the neural image informing behavior, we recorded snake responses evoked by targets moving across backgrounds composed of two contrasting temperatures with an average temperature equal to the target temperature. An unresolved background would appear uniform; thus, the target would be detectable only if the background pattern were resolved. Western rattlesnakes (Crotalus oreganus) displayed no statistically significant responses to targets presented in front of patterned backgrounds, regardless of the temperature contrasts or spatial frequencies within the background, but responded strongly to targets presented in front of homogeneous backgrounds. We found no evidence that the pit organ system can resolve spatial details subtending an angle of 9 deg or less. We discuss the implications of these results for understanding pit organ function in ecologically relevant habitats with thermal heterogeneity.

Atomic-Level Features for Kinetic Monte Carlo Models of Complex Chemistry from Molecular Dynamics Simulations.

The high computational cost of evaluating atomic interactions recently motivated the development of computationally inexpensive kinetic models, which can be parameterized from molecular dynamics (MD) simulations of the complex chemistry of thousands of species or other processes and accelerate the prediction of the chemical evolution by up to four orders of magnitude. Such models go beyond the commonly employed potential energy surface fitting methods in that they are aimed purely at describing kinetic effects. So far, such kinetic models utilize molecular descriptions of reactions and have been constrained to only reproduce molecules previously observed in MD simulations. Therefore, these descriptions fail to predict the reactivity of unobserved molecules, for example, in the case of large molecules or solids. Here, we propose a new approach for the extraction of reaction mechanisms and reaction rates from MD simulations, namely, the use of atomic-level features. Using the complex chemical network of hydrocarbon pyrolysis as an example, it is demonstrated that kinetic models built using atomic features are able to explore chemical reaction pathways never observed in the MD simulations used to parameterize them, a critical feature to describe rare events. Atomic-level features are shown to construct reaction mechanisms and estimate reaction rates of unknown molecular species from elementary atomic events. Through comparisons of the model ability to extrapolate to longer simulation time scales and different chemical compositions than the ones used for parameterization, it is demonstrated that kinetic models employing atomic features retain the same level of accuracy and transferability as the use of features based on molecular species, while being more compact and parameterized with less data. We also find that atomic features can better describe the formation of large molecules enabling the simultaneous description of small molecules and condensed phases.

New Assembly-Free Bulk Layered Inorganic Vertical Heterostructures with Infrared and Optical Bandgaps.

In principle, a nearly endless number of unique van der Waals heterostructures can be created through the vertical stacking of two-dimensional (2D) materials, resulting in unprecedented potential for material design. However, this widely employed synthetic method for generating van der Waals heterostructures is slow, imprecise, and prone to introducing interlayer contaminants when compared with synthesis methods that are scalable to industrially relevant scales. Herein, we study the properties of a new class of layered bulk inorganic materials that has recently been reported that we call assembly-free bulk layered inorganic heterostructures, wherein the individual layers are of dissimilar chemical composition, distinguishing them from commonly studied layered materials. We find that these bulk materials exhibit properties similar to vertical heterostructures but without the complex and unscalable stacking process. Using state-of-the-art computational approaches, we study the electronic properties of livingstonite (HgSb4S8), a naturally occurring mineral that is a bulk lattice-commensurate heterostructure. We find that isolated bilayers of livingstonite have an intralayer HSE-06 band gap of 2.08 eV. This is the first report of a naturally occurring van der Waals heterostructure with a calculated band gap in the visible spectrum. We also studied the electronic properties of tetragonal Ti3Bi4O12, Sm2Ti3Bi2O12, orthorhombic Ti3Bi4O12, Nb3Bi5O15, LaTiNbBi2O9, and AgPbBrO and found some of them are potentially well-suited for photovoltaic applications. We also provide characterization of the electronic structure of the isolated bilayer and monolayer subcomponents of the bulk heterostructures. The report of the properties of these materials significantly enhances the library of known van der Waals heterostructures.

Low Resistivity and High Breakdown Current Density of 10 nm Diameter van der Waals TaSe3 Nanowires by Chemical Vapor Deposition.

Micron-scale single-crystal nanowires of metallic TaSe3, a material that forms -Ta-Se3-Ta-Se3- stacks separated from one another by a tubular van der Waals (vdW) gap, have been synthesized using chemical vapor deposition (CVD) on a SiO2/Si substrate, in a process compatible with semiconductor industry requirements. Their electrical resistivity was found unaffected by downscaling from the bulk to as little as 7 nm in nanowire width and height, in striking contrast to the resistivity of copper for the same dimensions. While the bulk resistivity of TaSe3 is substantially higher than that of bulk copper, at the nanometer scale the TaSe3 wires become competitive to similar-sized copper ones. Moreover, we find that the vdW TaSe3 nanowires sustain current densities in excess of 108 A/cm2 and feature an electromigration energy barrier twice that of copper. The results highlight the promise of quasi-one-dimensional transition metal trichalcogenides for electronic interconnect applications and the potential of van der Waals materials for downscaled electronics.

Transferable Kinetic Monte Carlo Models with Thousands of Reactions Learned from Molecular Dynamics Simulations.

Molecular dynamics (MD) simulation of complex chemistry typically involves thousands of atoms propagating over millions of time steps, generating a wealth of data. Traditionally these data are used to calculate some aggregate properties of the system and then discarded, but we propose that these data can be reused to study related chemical systems. Using approximate chemical kinetic models and methods from statistical learning, we study hydrocarbon chemistries under extreme thermodynamic conditions. We discover that a single MD simulation can contain sufficient information about reactions and rates to predict the dynamics of related yet different chemical systems using kinetic Monte Carlo (KMC) simulation. Our learned KMC models identify thousands of reactions and run 4 orders of magnitude faster than MD. The transferability of these models suggests that we can viably reuse data from existing MD simulations to accelerate future simulation studies and reduce the number of new MD simulations required.

Screening billions of candidates for solid lithium-ion conductors: A transfer learning approach for small data.

Machine learning (ML) methods have the potential to revolutionize materials design, due to their ability to screen materials efficiently. Unlike other popular applications such as image recognition or language processing, large volumes of data are not available for materials design applications. Here, we first show that a standard learning approach using generic descriptors does not work for small data, unless it is guided by insights from physical equations. We then propose a novel method for transferring such physical insights onto more generic descriptors, allowing us to screen billions of unknown compositions for Li-ion conductivity, a scale which was previously unfeasible. This is accomplished by using the accurate model trained with physical insights to create a large database, on which we train a new ML model using the generic descriptors. Unlike previous applications of ML, this approach allows us to screen materials which have not necessarily been tested before (i.e., not on ICSD or Materials Project). Our method can be applied to any materials design application where a small amount of data is available, combined with high details of physical understanding.

ODE integration schemes for plane-wave real-time time-dependent density functional theory.

Integration schemes are implemented with a plane-wave basis in the context of real-time time-dependent density functional theory. Crank-Nicolson methods and three classes of explicit integration schemes are explored and assessed in terms of their accuracy and stability properties. Within the framework of plane-wave density functional theory, a graphene monolayer system is used to investigate the error, stability, and serial computational cost of these methods. The results indicate that Adams-Bashforth and Adams-Bashforth-Moulton methods of orders 4 and 5 outperform commonly used methods, including Crank-Nicolson and Runge-Kutta methods, in simulations where a relatively low error is desired. Parallel runtime scaling of the most competitive serial methods is presented, further demonstrating that the Adams-Bashforth and Adams-Bashforth-Moulton methods are efficient methods for propagating the time-dependent Kohn-Sham equations. Our integration schemes are implemented as an extension to the Quantum ESPRESSO code.

Data Mining for New Two- and One-Dimensional Weakly Bonded Solids and Lattice-Commensurate Heterostructures.

Layered materials held together by weak interactions including van der Waals forces, such as graphite, have attracted interest for both technological applications and fundamental physics in their layered form and as an isolated single-layer. Only a few dozen single-layer van der Waals solids have been subject to considerable research focus, although there are likely to be many more that could have superior properties. To identify a broad spectrum of layered materials, we present a novel data mining algorithm that determines the dimensionality of weakly bonded subcomponents based on the atomic positions of bulk, three-dimensional crystal structures. By applying this algorithm to the Materials Project database of over 50,000 inorganic crystals, we identify 1173 two-dimensional layered materials and 487 materials that consist of weakly bonded one-dimensional molecular chains. This is an order of magnitude increase in the number of identified materials with most materials not known as two- or one-dimensional materials. Moreover, we discover 98 weakly bonded heterostructures of two-dimensional and one-dimensional subcomponents that are found within bulk materials, opening new possibilities for much-studied assembly of van der Waals heterostructures. Chemical families of materials, band gaps, and point groups for the materials identified in this work are presented. Point group and piezoelectricity in layered materials are also evaluated in single-layer forms. Three hundred and twenty-five of these materials are expected to have piezoelectric monolayers with a variety of forms of the piezoelectric tensor. This work significantly extends the scope of potential low-dimensional weakly bonded solids to be investigated.

Nanosecond homogeneous nucleation and crystal growth in shock-compressed SiO2.

Understanding the kinetics of shock-compressed SiO2 is of great importance for mitigating optical damage for high-intensity lasers and for understanding meteoroid impacts. Experimental work has placed some thermodynamic bounds on the formation of high-pressure phases of this material, but the formation kinetics and underlying microscopic mechanisms are yet to be elucidated. Here, by employing multiscale molecular dynamics studies of shock-compressed fused silica and quartz, we find that silica transforms into a poor glass former that subsequently exhibits ultrafast crystallization within a few nanoseconds. We also find that, as a result of the formation of such an intermediate disordered phase, the transition between silica polymorphs obeys a homogeneous reconstructive nucleation and grain growth model. Moreover, we construct a quantitative model of nucleation and grain growth, and compare its predictions with stishovite grain sizes observed in laser-induced damage and meteoroid impact events.

Dynamic Structural Response and Deformations of Monolayer MoS2 Visualized by Femtosecond Electron Diffraction.

Two-dimensional materials are subject to intrinsic and dynamic rippling that modulates their optoelectronic and electromechanical properties. Here, we directly visualize the dynamics of these processes within monolayer transition metal dichalcogenide MoS2 using femtosecond electron scattering techniques as a real-time probe with atomic-scale resolution. We show that optical excitation induces large-amplitude in-plane displacements and ultrafast wrinkling of the monolayer on nanometer length-scales, developing on picosecond time-scales. These deformations are associated with several percent peak strains that are fully reversible over tens of millions of cycles. Direct measurements of electron-phonon coupling times and the subsequent interfacial thermal heat flow between the monolayer and substrate are also obtained. These measurements, coupled with first-principles modeling, provide a new understanding of the dynamic structural processes that underlie the functionality of two-dimensional materials and open up new opportunities for ultrafast strain engineering using all-optical methods.

Strain engineering in monolayer materials using patterned adatom adsorption.

We utilize reactive empirical bond order (REBO)-based interatomic potentials to explore the potential for the engineering of strain in monolayer materials using lithographically or otherwise patterned adatom adsorption. In the context of graphene, we discover that the monolayer strain results from a competition between the in-plane elasticity and out-of-plane relaxation deformations. For hydrogen adatoms on graphene, the strain outside the adsorption region vanishes due to out-of-plane relaxation deformations. Under some circumstances, an annular adsorption pattern generates homogeneous tensile strains of approximately 2% in graphene inside the adsorption region, approximately 30% of the strain in the adsorbed region. We find that an elliptical adsorption pattern produces strains of as large as 5%, close to the strain in the adsorbed region. Also, nonzero maximum shear strain (∼ 4%) can be introduced by the elliptical adsorption pattern. We find that an elastic plane stress model provides qualitative guidance for strain magnitudes and conditions under which strain-diminishing buckling can be avoided. We identify geometric conditions under which this effect has potential to be scaled to larger areas. Our results elucidate a method for strain engineering at the nanoscale in monolayer devices.

Flexural electromechanical coupling: a nanoscale emergent property of boron nitride bilayers.

The symmetry properties of atomically thin boron nitride (BN) monolayers endow them with piezoelectric properties, whereas the bulk parent crystal of stacked BN layers is not piezoelectric. This suggests potential for unusual electromechanical properties in the few layer regime. In this work, we explore this regime and discover that a bilayer consisting of two BN monolayers exhibits a strong mechanical coupling between curvature and electric fields. Using a mechanical model with parameters obtained from density functional theory, we find that these bilayers amplify in-plane piezoelectric displacements by exceedingly large factors on the order of 10(3)-10(4). We find that this type of electromechanical coupling is an emergent nanoscale property that occurs only for the case of two stacked BN monolayers.

Optical characterization of chemistry in shocked nitromethane with time-dependent density functional theory.

We compute the optical properties of the liquid-phase energetic material nitromethane (CH3NO2) for the first 100 ps behind the front of a simulated shock at 6.5 km/s, close to the experimentally observed detonation shock speed of the material. We utilize molecular dynamics trajectories computed using the multiscale shock technique (MSST) for time-resolved optical spectrum calculations based on both linear response time-dependent DFT (TDDFT) and the Kubo-Greenwood formula with Kohn-Sham DFT wave functions. We find that the TDDFT method predicts an optical conductivity 25-35% lower than the Kubo-Greenwood calculation and provides better agreement with the experimentally measured index of refraction of unreacted nitromethane. We investigate the influence of electronic temperature on the Kubo-Greenwood spectra and find no significant effect at optical wavelengths. In both Kubo-Greenwood and TDDFT, the spectra evolve nonmonotonically in time as shock-induced chemistry takes place. We attribute the time-resolved absorption at optical wavelengths to time-dependent populations of molecular decomposition products, including NO, CNO, CNOH, H2O, and larger molecules. These calculations offer direction for guiding and interpreting ultrafast optical measurements on reactive materials.

Comparison of ReaxFF, DFTB, and DFT for phenolic pyrolysis. 2. Elementary reaction paths.

Reaction paths for the loss of CO, H2, and H2O from atomistic models of phenolic resin are determined using the hybrid B3LYP approach. B3LYP energetics are confirmed using CCSD(T). The energetics along the B3LYP paths are also evaluated using the PW91 generalized gradient approximation (GGA), the more approximate self-consistent charge density functional tight binding (SCC-DFTB), and the reactive force field (ReaxFF). Compared with the CCSD(T)/cc-pVTZ level for bond and reaction energies and barrier heights, the B3LYP, PW91, DFTB(mio), DFTB(pbc), and ReaxFF have average absolute errors of 3.8, 5.1, 17.4, 13.2, and 19.6 kcal/mol, respectively. The PW91 is only slightly less accurate than the B3LYP approach, while the more approximate approaches yield somewhat larger errors. The SCC-DFTB paths are in better agreement with B3LYP than are those obtained with ReaxFF.

Comparison of ReaxFF, DFTB, and DFT for phenolic pyrolysis. 1. Molecular dynamics simulations.

A systematic comparison of atomistic modeling methods including density functional theory (DFT), the self-consistent charge density-functional tight-binding (SCC-DFTB), and ReaxFF is presented for simulating the initial stages of phenolic polymer pyrolysis. A phenolic polymer system is simulated for several hundred picoseconds within a temperature range of 2500 to 3500 K. The time evolution of major pyrolysis products including small-molecule species and char is examined. Two temperature zones are observed which demark cross-linking versus fragmentation. The dominant chemical products for all methods are similar, but the yields for each product differ. At 3500 K, DFTB overestimates CO production (300-400%) and underestimates free H (~30%) and small C(m)H(n)O molecules (~70%) compared with DFT. At 3500 K, ReaxFF underestimates free H (~60%) and fused carbon rings (~70%) relative to DFT. Heterocyclic oxygen-containing five- and six-membered carbon rings are observed at 2500 K. Formation mechanisms for H2O, CO, and char are discussed. Additional calculations using a semiclassical method for incorporating quantum nuclear energies of molecules were also performed. These results suggest that chemical equilibrium can be affected by quantum nuclear effects at temperatures of 2500 K and below. Pyrolysis reaction mechanisms and energetics are examined in detail in a companion manuscript.