Search

Hybrid organic-inorganic nanomaterials can exhibit transitional behavior that deviates from models developed for all-organic or all-inorganic materials systems. Here, we reveal the complexity of exciton-phonon interactions in a recently discovered 2D layered hybrid organic-inorganic semiconductor, silver phenylselenolate (AgSePh). Using femtosecond resonant impulsive vibrational spectroscopy and non-resonant Raman scattering, we measure multiple hybrid organic-inorganic vibrational modes and identify a subset of these modes that strongly couple to the electronic excited state. Calculations by density functional perturbation theory show that these strongly coupled modes exhibit large out-of-plane silver atomic motions and silver-selenium spacing displacements. Moreover, analysis of photoluminescence fine-structure splitting and temperature-dependent peak-shifting/linewidth-broadening suggests that light emission in AgSePh is most strongly affected by a compound 100 cm-1 mode involving the wagging motion of phenylselenolate ligands and accompanying metal-chalcogen stretching. Finally, red-shifting of vibrational modes with increasing temperature reveals a high degree of anharmonicity arising from non-covalent interactions between phenyl rings. These findings reveal the unique effects of hybrid vibrational modes in organic-inorganic semiconductors and motivate future work aimed at specifically engineering such interactions through chemical and structural modification.

The widespread application of multimodal machine learning models like GPT-4 has revolutionized various research fields including computer vision and natural language processing. However, its implementation in materials informatics remains underexplored, despite the presence of materials data across diverse modalities, such as composition and structure. The effectiveness of machine learning models trained on large calculated datasets depends on the accuracy of calculations, while experimental datasets often have limited data availability and incomplete information. This paper introduces a novel approach to multimodal machine learning in materials science via composition-structure bimodal learning. The proposed COmposition-Structure Bimodal Network (COSNet) is designed to enhance learning and predictions of experimentally measured materials properties that have incomplete structure information. Bimodal learning significantly reduces prediction errors across distinct materials properties including Li conductivity in solid electrolyte, band gap, refractive index, dielectric constant, energy, and magnetic moment, surpassing composition-only learning methods. Furthermore, we identified that data augmentation based on modal availability plays a pivotal role in the success of bimodal learning.

Historically, materials informatics has relied on human-designed descriptors of materials structures. In recent years, graph neural networks (GNNs) have been proposed for learning representations of crystal structures from data end-to-end producing vectorial embeddings that are optimized for downstream prediction tasks. However, a systematic scheme is lacking to analyze and understand the limits of GNNs for capturing crystal structures. In this work, we propose to use human-designed descriptors as a bank of human knowledge to test whether black-box GNNs can capture the knowledge of crystal structures. We find that current state-of-the-art GNNs cannot capture the periodicity of crystal structures well, and we analyze the limitations of the GNN models that result in this failure from three aspects: local expressive power, long-range information, and readout function. We propose an initial solution, hybridizing descriptors with GNNs, to improve the prediction of GNNs for materials properties, especially phonon internal energy and heat capacity with 90% lower errors, and we analyze the mechanisms for the improved prediction. All the analysis can be extended easily to other deep representation learning models, human-designed descriptors, and systems such as molecules and amorphous materials.

Open material databases storing hundreds of thousands of material structures and their corresponding properties have become the cornerstone of modern computational materials science. Yet, the raw outputs of the simulations, such as the trajectories from molecular dynamics simulations and charge densities from density functional theory calculations, are generally not shared due to their huge size. In this work, we describe a cloud-based platform to facilitate the sharing of raw data and enable the fast post-processing in the cloud to extract new properties defined by the user. As an initial demonstration, our database currently includes 6286 molecular dynamics trajectories for amorphous polymer electrolytes and 5.7 terabytes of data. We create a public analysis library at https://github.com/TRI-AMDD/htp_md to extract multiple properties from the raw data, using both expert designed functions and machine learning models. The analysis is run automatically with computation in the cloud, and results then populate a database that can be accessed publicly. Our platform encourages users to contribute both new trajectory data and analysis functions via public interfaces. Newly analyzed properties will be incorporated into the database. Finally, we create a front-end user interface at https://www.htpmd.matr.io for browsing and visualization of our data. We envision the platform to be a new way of sharing raw data and new insights for the computational materials science community., Comment: 21 pages, 7 figures

Machine learning materials properties measured by experiments is valuable yet difficult due to the limited amount of experimental data. In this work, we use a multi-fidelity random forest model to learn the experimental formation enthalpy of materials with prediction accuracy higher than the empirically corrected PBE functional (PBEfe) and meta-GGA functional (SCAN), and it outperforms the hotly studied deep neural-network based representation learning and transfer learning. We then use the model to calibrate the DFT formation enthalpy in the Materials Project database, and discover materials with underestimated stability. The multi-fidelity model is also used as a data-mining approach to find how DFT deviates from experiments by the explaining the model output.

Polymer electrolytes are promising candidates for the next generation lithium-ion battery technology. Large scale screening of polymer electrolytes is hindered by the significant cost of molecular dynamics (MD) simulation in amorphous systems: the amorphous structure of polymers requires multiple, repeated sampling to reduce noise and the slow relaxation requires long simulation time for convergence. Here, we accelerate the screening with a multi-task graph neural network that learns from a large amount of noisy, unconverged, short MD data and a small number of converged, long MD data. We achieve accurate predictions of 4 different converged properties and screen a space of 6247 polymers that is orders of magnitude larger than previous computational studies. Further, we extract several design principles for polymer electrolytes and provide an open dataset for the community. Our approach could be applicable to a broad class of material discovery problems that involve the simulation of complex, amorphous materials., Comment: 29 pages, 6 figures + supplementary information

Two-dimensional (2D) materials have received considerable attention as possible electrodes in Li-ion batteries (LIBs), although a deeper understanding of the Li adsorption behavior as well as broad screening of the materials space is still needed. In this work, we build a high-throughput screening scheme that incorporates a learned interaction. First, density functional theory and graph convolution networks are utilized to calculate minimum Li adsorption energies for a small set of 2D metallic materials. The data is then used to find a dependence of the minimum Li adsorption energies on the sum of ionization potential, work function of the 2D metal, and coupling energy between Li+ and substrate. Our results show that variances of elemental properties and density are the most correlated features with coupling. To illustrate the applicability of this approach, the model is employed to show that some fluorides and chromium oxides are potential high-voltage materials with adsorption energies < -7 eV, and the found physics is used as the design principle to enhance the Li adsorption ability of graphene. This physics-driven approach shows higher accuracy and transferability compared with purely data-driven models.

Self-affine morphology of random interfaces governs their functionalities across tribological, geological, (opto-)electrical and biological applications. However, the knowledge of how energy carriers or generally classical/quantum waves interact with structural irregularity is still incomplete. In this work, we study vibrational energy transport through random interfaces exhibiting different correlation functions on the two-dimensional hexagonal lattice. We show that random interfaces at the atomic scale are Cantor composites populated on geometrical fractals, thus multifractals, and calculate their quantized conductance using atomistic approaches. We obtain a universal scaling law, which contains self-similarity for mass perturbation, and exponential scaling of structural irregularity quantified by fractal dimension. The multifractal nature and Cantor-composite picture may also be extendable to charge and photon transport across random interfaces.

Thermal conductivity is a fundamental material property but challenging to predict, with less than 5% out of about $10^5$ synthesized inorganic materials being documented. In this work, we extract the structural chemistry that governs lattice thermal conductivity, by combining graph neural networks and random forest approaches. We show that both mean and variation of unit-cell configurational properties, such as atomic volume and bond length, are the most important features, followed by mass and elemental electronegativity. We chart the structural chemistry of lattice thermal conductivity into extended van-Arkel triangles, and predict the thermal conductivity of all known inorganic materials in the Inorganic Crystal Structure Database. For the latter, we develop a transfer learning framework extendable for other applications.

Electronic structure of layered LiNiO2 has been controversial despite numerous theoretical and experimental reports regarding its nature. We investigate the charge densities, lithium intercalation potentials and Li diffusion barrier energies of LixNiO2 (0.0 < x < 1.0) system using a truly ab-initio method, diffusion quantum Monte Carlo (DMC). We compare the charge densities from DMC and density functional theory (DFT) and show that local and semi-local DFT functionals yield spin polarization densities with incorrect sign on the oxygen atoms. SCAN functional and Hubbard-U correction improves the polarization density around Ni and O atoms, resulting in smaller deviations from the DMC densities. DMC accurately captures the p-d hybridization between the Ni-O atoms, yielding accurate lithium intercalation voltages, polarization densities and reaction barriers.

Electron charge density distribution of materials is one of the key quantities in computational materials science as theoretically it determines the ground state energy and practically it is used in many materials analyses. However, the scaling of density functional theory calculations with number of atoms limits the usage of charge-density-based calculations and analyses. Here we introduce a machine learning scheme with local-environment-based graphs and graph convolutional neural networks to predict charge density on grid-points from crystal structure. We show the accuracy of this scheme through a comparison of predicted charge densities as well as properties derived from the charge density, and the scaling is O(N). More importantly, the transferability is shown to be high with respect to different compositions and structures, which results from the explicit encoding of geometry.

Understanding the dynamical processes that govern the performance of functional materials is essential for the design of next generation materials to tackle global energy and environmental challenges. Many of these processes involve the dynamics of individual atoms or small molecules in condensed phases, e.g. lithium ions in electrolytes, water molecules in membranes, molten atoms at interfaces, etc., which are difficult to understand due to the complexity of local environments. In this work, we develop graph dynamical networks, an unsupervised learning approach for understanding atomic scale dynamics in arbitrary phases and environments from molecular dynamics simulations. We show that important dynamical information can be learned for various multi-component amorphous material systems, which is difficult to obtain otherwise. With the large amounts of molecular dynamics data generated everyday in nearly every aspect of materials design, this approach provides a broadly useful, automated tool to understand atomic scale dynamics in material systems., Comment: 25 + 7 pages, 5 + 3 figures

We present a new approximation to ionic conductivity well suited to dynamical atomic-scale simulations, based on the Nernst-Einstein equation. In our approximation, ionic aggregates constitute the elementary charge carriers, and are considered as non-interacting species. This approach conveniently captures the dominant effect of ion-ion correlations on conductivity, short range interactions in the form of clustering. In addition to providing better estimates to the conductivity at a lower computational cost than exact approaches, this new method allows to understand the physical mechanisms driving ion conduction in concentrated electrolytes. As an example, we consider Li$^+$ conduction in poly(ethylene oxide), a standard solid-state polymer electrolyte. Using our newly developed approach, we are able to reproduce recent experimental results reporting negative cation transference numbers at high salt concentrations, and to confirm that this effect can be caused by a large population of negatively charged clusters involving cations.

Solar thermal fuels (STFs) are an unconventional paradigm for solar energy conversion and storage which is attracting renewed attention. In this concept, a material absorbs sunlight and stores the energy chemically via an induced structural change, which can later be reversed to release the energy as heat. An example is the azobenzene molecule which has a cis-trans photoisomerization with these properties, and can be tuned by chemical substitution and attachment to templates such as carbon nanotubes, small molecules, or polymers. By analogy to the Shockley-Queisser limit for photovoltaics, we analyze the maximum attainable efficiency for STFs from fundamental thermodynamic considerations. Microscopic reversibility provides a bound on the quantum yield of photoisomerization due to fluorescence, regardless of details of photochemistry. We emphasize the importance of analyzing the free energy, not just enthalpy, of the metastable molecules, and find an efficiency limit for conversion to stored chemical energy equal to the Shockley-Queisser limit. STF candidates from a recent high-throughput search are analyzed in light of the efficiency limit., Comment: 16 pages, 4 figures

The combination of high throughput computation and machine learning has led to a new paradigm in materials design by allowing for the direct screening of vast portions of structural, chemical, and property space. The use of these powerful techniques leads to the generation of enormous amounts of data, which in turn calls for new techniques to efficiently explore and visualize the materials space to help identify underlying patterns. In this work, we develop a unified framework to hierarchically visualize the compositional and structural similarities between materials in an arbitrary material space with representations learned from different layers of graph convolutional neural networks. We demonstrate the potential for such a visualization approach by showing that patterns emerge automatically that reflect similarities at different scales in three representative classes of materials: perovskites, elemental boron, and general inorganic crystals, covering material spaces of different compositions, structures, and both. For perovskites, elemental similarities are learned that reflects multiple aspects of atom properties. For elemental boron, structural motifs emerge automatically showing characteristic boron local environments. For inorganic crystals, the similarity and stability of local coordination environments are shown combining different center and neighbor atoms. The method could help transition to a data-centered exploration of materials space in automated materials design., Comment: 22 + 7 pages, 6 + 5 figures

Next generation batteries based on lithium (Li) metal anodes have been plagued by the dendritic electrodeposition of Li metal on the anode during cycling, resulting in short circuit and capacity loss. Suppression of dendritic growth through the use of solid electrolytes has emerged as one of the most promising strategies for enabling the use of Li metal anodes. We perform a computational screening of over 12,000 inorganic solids based on their ability to suppress dendrite initiation in contact with Li metal anode. Properties for mechanically isotropic and anisotropic interfaces that can be used in stability criteria for determining the propensity of dendrite initiation are usually obtained from computationally expensive first-principles methods. In order to obtain a large dataset for screening, we use machine learning models to predict the mechanical properties of several new solid electrolytes. We train a convolutional neural network on the shear and bulk moduli purely on structural features of the material. We use AdaBoost, Lasso and Bayesian ridge regression to train the elastic constants, where the choice of the model depended on the size of the training data and the noise that it can handle. Our models give us direct interpretability by revealing the dominant structural features affecting the elastic constants. The stiffness is found to increase with a decrease in volume per atom, increase in minimum anion-anion separation, and increase in sublattice (all but Li) packing fraction. Cross-validation/test performance suggests our models generalize well. We predict over 20 mechanically anisotropic interfaces between Li metal and 6 solid electrolytes which can be used to suppress dendrite growth. Our screened candidates are generally soft and highly anisotropic, and present opportunities for simultaneously obtaining dendrite suppression and high ionic conductivity in solid electrolytes., Comment: 34 pages, 4 Figures, 3 Table, 7 pages of Supporting Information

Despite their potential to exceed the theoretical Shockley-Queisser limit, ferroelectric photovoltaics (FPVs) have performed inefficiently due to their extremely low photocurrents. Incorporating Bi2FeCrO6 (BFCO) as the light absorber in FPVs has recently led to impressively high and record photocurrents [Nechache et al. Nature Photon. 2015, 9, 61], reviving the FPV field. However, our understanding of this remarkable phenomenon is far from satisfactory. Here, we use first-principles calculations to determine that such excellent performance mainly lies in the efficient separation of electron-hole (e-h) pairs. We show that photoexcited electrons and holes in BFCO are spatially separated on the Fe and Cr sites, respectively. This separation is much more pronounced in disordered BFCO phases, which show exceptional PV responses. We further set out to design a strategy for next-generation FPVs, not limited to BFCO, by exploring 44 additional Bi-based double-perovskite oxides. We suggest 9 novel active-layer materials that can offer strong e-h separations and a desired band gap energy for application in FPVs. Our work indicates that charge separation is the most important issue to be addressed for FPVs to compete with conventional devices., Comment: 33 pages, 5 main figures, 5 supporting figures

The use of machine learning methods for accelerating the design of crystalline materials usually requires manually constructed feature vectors or complex transformation of atom coordinates to input the crystal structure, which either constrains the model to certain crystal types or makes it difficult to provide chemical insights. Here, we develop a crystal graph convolutional neural networks framework to directly learn material properties from the connection of atoms in the crystal, providing a universal and interpretable representation of crystalline materials. Our method provides a highly accurate prediction of density functional theory calculated properties for eight different properties of crystals with various structure types and compositions after being trained with $10^4$ data points. Further, our framework is interpretable because one can extract the contributions from local chemical environments to global properties. Using an example of perovskites, we show how this information can be utilized to discover empirical rules for materials design., Comment: 6+9 pages, 3+6 figures

Nanoporous materials are a promising platform for thermoelectrics in that they offer high thermal conductivity tunability while preserving good electrical properties, a crucial requirement for high- effciency thermal energy conversion. Understanding the impact of the pore arrangement on thermal transport is pivotal to engineering realistic materials, where pore disorder is unavoidable. Although there has been considerable progress in modeling thermal size effects in nanostructures, it has remained a challenge to screen such materials over a large phase space due to the slow simulation time required for accurate results. We use density functional theory in connection with the Boltzmann transport equation, to perform calculations of thermal conductivity in disordered porous materials. By leveraging graph theory and regressive analysis, we identify the set of pores representing the phonon bottleneck and obtain a descriptor for thermal transport, based on the sum of the pore-pore distances between such pores. This approach provides a simple tool to estimate phonon suppression in realistic porous materials for thermoelectric applications and enhance our understanding of heat transport in disordered materials.

Rapid, non-destructive characterization of molecular level chemistry for organic matter (OM) is experimentally challenging. Raman spectroscopy is one of the most widely used techniques for non-destructive chemical characterization, although it currently does not provide detailed identification of molecular components in OM, due to the combination of diffraction-limited spatial resolution and poor applicability of peak-fitting algorithms. Here, we develop a genome-inspired collective molecular structure fingerprinting approach, which utilizes ab initio calculations and data mining techniques to extract molecular level chemistry from the Raman spectra of OM. We illustrate the power of such an approach by identifying representative molecular fingerprints in OM, for which the molecular chemistry is to date inaccessible using non-destructive characterization techniques. Chemical properties such as aromatic cluster size distribution and H/C ratio can now be quantified directly using the identified molecular fingerprints. Our approach will enable non-destructive identification of chemical signatures with their correlation to the preservation of biosignatures in OM, accurate detection and quantification of environmental contamination, as well as objective assessment of OM with respect to their chemical contents.

Pressure-induced changes in the Fermi surface of Cd up to 40 GPa are studied using highly accurate density-functional theory calculations. The topology of Fermi surface changes at pressures of 2, 8, 12 , 18 and 28 GPa indicating electronic topological transitions (ETTs). Structural parameters, compressibility data and elastic constant reveal anomalies across these ETTs. The computed equation of state at 300 K is in excellent agreement with the experimental data. In view of the highly controversial nature of these ETTs the present studies are aimed at motivating the experimentalists for their direct detection at high pressures., Comment: 5 figures

Strain in a material induces shifts in vibrational frequencies, which is a probe of the nature of the vibrations and interatomic potentials, and can be used to map local stress/strain distributions via Raman microscopy. This method is standard for crystalline silicon devices, but due to lack of calibration relations, it has not been applied to amorphous materials such as hydrogenated amorphous silicon (a-Si:H), a widely studied material for thin-film photovoltaic and electronic devices. We calculated the Raman spectrum of a-Si:H \ab initio under different strains $\epsilon$ and found peak shifts $\Delta \omega = \left( -460 \pm 10\ \mathrm{cm}^{-1} \right) {\rm Tr}\ \epsilon$. This proportionality to the trace of the strain is the general form for isotropic amorphous vibrational modes, as we show by symmetry analysis and explicit computation. We also performed Raman measurements under strain and found a consistent coefficient of $-510 \pm 120\ \mathrm{cm}^{-1}$. These results demonstrate that a reliable calibration for the Raman/strain relation can be achieved even for the broad peaks of an amorphous material, with similar accuracy and precision as for crystalline materials., Comment: 12 pages, 3 figures + supplementary 8 pages, 4 figures

In this book chapter we provide the definition of "Simulating Nanoscale Heat Transport", broadly regarded as modeling heat conduction beyond Fourier's law. We primarely focus on incoherent transport, which is dominated by scattering between phonons and the boundaries of a given material. After introducing the Boltzmann Transport Equation (BTE) for phonons, we elaborate on different numerical techniques for solving nanoscale heat transport includind deterministic and Monte Carlo based method. A discussion on the need for multiscale modeling concludes the chapter.

The evolution of the electronic absorption edge of type I, II and III kerogen is studied by diffuse reflectance UV-Visible absorption spectroscopy. The functional form of the electronic absorption edge for all kerogens measured is in excellent agreement with the "Urbach tail" phenomenology. The Urbach decay width extracted from the exponential fit within the visible range is strongly correlated with the aliphatic/aromatic ratio in isolated kerogen, regardless of the kerogen type. No correlation is found between the decay width and the average size of aromatic clusters, which is explained in terms of a non-linear increase in optical absorption with increasing size of the aromatic clusters determined by 13C NMR. Further, absorption spectra calculated with density functional theory calculations on proxy ensemble models of kerogen are in excellent agreement with the experimental results. The correlation of the decay width with conventional maturity indicators such as vitrinite reflectance is found to be good within a particular kerogen type, but not consistent across different kerogen types, reflecting systematic variations in bulk composition for different type kerogen types with the same vitrinite reflectance. Thus, diffuse reflectance visible absorption spectroscopy is presented as a rapid, calibrated and non-destructive method to monitor both the maturity and the chemical composition of kerogen. The chemical insight of kerogen in relation to its optical absorption provided by this methodology may serve for rapid screening of kerogen for electronics and optical devices in place of functionalized produced carbon., Comment: 26 pages, 10 figures. To appear in "Carbon"

Tuning thermal transport in nanostructured materials is a powerful approach to develop high-efficiency thermoelectric materials. Using a recently developed approach based on the phonon mean free path dependent Boltzmann transport equation, we compute the effective thermal conductivity of nanoporous materials with pores of various shapes and arrangements. We assess the importance of pore-pore distance in suppressing thermal transport, and identify the pore arrangement that minimizes the thermal conductivity, composed of a periodic arrangement of two misaligned rows of triangular pores. Such a configuration yields a reduction in the thermal conductivity of more than $60 \%$ with respect the simple circular aligned case with the same porosity., Comment: 4 pages, 4 figures, 1 table

By sculpting the magnetic field applied to magneto-acoustic materials, phonons can be used for information processing. Using a combination of analytic and numerical techniques, we demonstrate designs for diodes (isolators) and transistors that are independent of their conventional, electronic formulation. We analyze the experimental feasibility of these systems, including the sensitivity of the circuits to likely systematic and random errors., Comment: 5 pages, 4 figures

We develop a computational framework, based on the Boltzmann transport equation, with the ability to compute the thermal transport in nanostructured materials of any geometry using as the only input the bulk thermal conductivity accumulation function. The main advantage of our method is twofold. First, while the scattering times and dispersion curves are unknown for most materials, the phonon mean free path distribution can be directly obtained by experiments. As a consequence, a wider range of materials can be simulated than with a frequency-dependent approach. Second, when phonon dispersions are available from first principles calculations, our approach allows one to include easily the whole Brillouen zone in the calculations without discretizing the phonon frequencies for all polarizations, reducing considerably the computational effort. Furthermore, after deriving the ballistic and diffusive limits of our model, we develop a multi-scale scheme that couples phonon transport across different scales, enabling efficient simulations of materials with wide phonon mean free path distributions. After we validate our model against the frequency-dependent approach, we apply the method to porous silicon membranes and find good agreement with experiments on mesoscale pores. By enabling the investigation of thermal transport in unexplored nanostructured materials, our method has potential to advance high-efficiency thermoelectric devices.

Phononic computing -- the use of (typically thermal) vibrations for information processing -- is a nascent technology; its capabilities are still being discovered. We analyze an alternative form of phononic computing inspired by optical, rather than electronic, computing. Using the acoustic Faraday effect, we design a phonon gyrator and thereby a means of performing computation through the manipulation of polarization in transverse phonon currents. Moreover, we establish that our gyrators act as generalized transistors and can construct digital logic gates. Exploiting the wave nature of phonons and the similarity of our logic gates, we demonstrate parallel computation within a single circuit, an effect presently unique to phonons. Finally, a generic method of designing these parallel circuits is introduced and used to analyze the feasibility of magneto-acoustic materials in realizing these circuits.

We present an approach to calculation of point defect optical and thermal ionization energies based on the highly accurate quantum Monte Carlo methods. The use of an inherently many-body theory that directly treats electron correlation offers many improvements over the typically-employed density functional theory Kohn-Sham description. In particular, the use of quantum Monte Carlo methods can help overcome the band gap problem and obviate the need for ad-hoc corrections. We demonstrate our approach to the calculation of the optical and thermal ionization energies of the F-center defect in magnesium oxide, and obtain excellent agreement with experimental and/or other high-accuracy computational results.