Search

Hybrid organic–inorganic formate perovskites, AB(HCOO)3, are a large family of compounds that exhibit a variety of phase transitions and diverse properties, such as (anti)ferroelectricity, ferroelasticity, (anti)ferromagnetism, and multiferroism. While many properties of these materials have already been characterized, we are not aware of any study that focuses on the comprehensive property assessment of a large number of formate perovskites. A comparison of the properties of materials within the family is challenging due to systematic errors attributed to different techniques or the lack of data. For example, complete piezoelectric, dielectric, and elastic tensors are not available. In this work, we utilize first-principles density functional theory based simulations to overcome these challenges and to report structural, mechanical, dielectric, piezoelectric, and ferroelectric properties of 29 formate perovskites. We find that these materials exhibit elastic stiffness in the range 0.5–127.0 GPa; highly anisotropic linear compressibility, including zero and even negative values; dielectric constants in the range 0.1–102.1; highly anisotropic piezoelectric response with the longitudinal values in the range 1.18–21.12 pC/N; and spontaneous polarizations in the range 0.2–7.8 μC/cm2. Furthermore, we propose and computationally characterize a few formate perovskites that have not been reported yet. [ABSTRACT FROM AUTHOR]

Kohn–Sham density functional theory (DFT) is nowadays widely used for electronic structure theory simulations, and the accuracy and efficiency of DFT rely on approximations of the exchange-correlation functional. By including the kinetic energy density τ, the meta-generalized-gradient approximation (meta-GGA) family of functionals achieves better accuracy and flexibility while retaining the efficiency of semi-local functionals. For example, the strongly constrained and appropriately normed (SCAN) meta-GGA functional has been proven to yield accurate results for solid and molecular systems. We implement meta-GGA functionals with both numerical atomic orbitals and plane wave bases in the ABACUS package. Apart from the exchange-correlation potential, we also discuss the evaluation of force and stress. To validate our implementation, we perform finite-difference tests and convergence tests with the SCAN, rSCAN, and r2SCAN meta-GGA functionals. We further test water hexamers, weakly interacting molecules from the S22 dataset, as well as 13 semiconductors using the three functionals. The results show satisfactory agreement with previous calculations and available experimental values. [ABSTRACT FROM AUTHOR]

A polarization consistent framework, where dielectric screening is affected consistently in polarizable continuum model (PCM) calculations, is employed for the study of solvation energies. The computational framework combines a screened range-separated-hybrid functional (SRSH) with PCM calculations, SRSH-PCM, where dielectric screening is imposed in both PCM self-consistent reaction field (SCRF) iterations and the electronic structure Hamiltonian. We begin by demonstrating the impact of modifying the Hamiltonian to include such dielectric screening in SCRF iterations by considering the solutions of electrostatically embedded Hartree–Fock (HF) exact exchange equations. Long-range screened HF-PCM calculations are shown to capture properly the linear dependence of gap energy of frontier orbitals on the inverse of the dielectric constant, whereas unscreened HF-PCM orbital energies are fallaciously semi-constant with respect to the dielectric constant and, therefore, inconsistent with the ionization energy gaps. Similar trends affect density functional theory (DFT) calculations that aim to achieve predictive quality. Importantly, the dielectric screened calculations are shown to significantly affect DFT- and HF PCM-based solvation energies, where screened solvation energies are smaller compared to the unscreened values. Importantly, SRSH-PCM, therefore, appears to reduce the tendency of DFT-PCM to overestimate solvation energies, where we find the effect to increase with the dielectric constant and the polarity of the molecular solute, trends that enhance the quality of DFT-PCM calculations of solvation energy. Understanding the relationship of dielectric screening in the Hamiltonian and DFT-PCM calculations can ultimately benefit on-going efforts for the design of predictive and parameter free descriptions of solvation energies. [ABSTRACT FROM AUTHOR]

A lattice-based version of the aggregation-volume-bias Monte Carlo method that was introduced recently has allowed for the extension of the calculation of the nucleation free energies from liquid clusters to solid clusters. Here, it was used to calculate the nucleation free energies of both bcc and fcc clusters formed by Lennard–Jones particles. Under the simulation conditions considered in this study, a cross-over of the thermodynamic stability from the bcc to the fcc structure was observed directly from the free energy results. In addition, the free energies obtained for both types of clusters were used to extrapolate bulk phase information, including chemical potential and surface tension, which revealed that bcc clusters are favored due to the lower surface tension. These results corroborate a recent classical density functional theory study. This work also demonstrates that this approach can be used to predict the entire thermodynamic landscape (i.e., free energies for clusters of different structures and sizes, including an infinitely large cluster, which is the bulk phase), which is important to answer fundamental questions related to crystallization such as the origin of polymorphism. [ABSTRACT FROM AUTHOR]

The knowledge of the frontier orbital, highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO), energies is vital for studying chemical and electrochemical stability of compounds, their corrosion inhibition potential, reactivity, etc. Density functional theory (DFT) calculations provide a direct route to estimate these energies either in the gas-phase or condensed phase. However, the application of DFT methods becomes computationally intensive when hundreds of thousands of compounds are to be screened. Such is the case when all the isomers for the 1-alkyl-3-alkylimidazolium cation [CnCmim]+ (n = 1–10, m = 1–10) are considered. Enumerating the isomer space of [CnCmim]+ yields close to 386 000 cation structures. Calculating frontier orbital energies for each would be computationally very expensive and time-consuming using DFT. In this article, we develop a machine learning model based on the extreme gradient boosting method using a small subset of the isomer space and predict the HOMO and LUMO energies. Using the model, the HOMO energies are predicted with a mean absolute error (MAE) of 0.4 eV and the LUMO energies are predicted with a MAE of 0.2 eV. Inferences are also drawn on the type of the descriptors deemed important for the HOMO and LUMO energy estimates. Application of the machine learning model results in a drastic reduction in computational time required for such calculations. [ABSTRACT FROM AUTHOR]

Despite the wide applications, the ab initio modeling of the ceria based catalyst is challenging. The partial occupation in the 4f orbitals creates a fundamental challenge for commonly used density functional theory (DFT) methods, including semilocal functionals with Hubbard U correction to force localization and hybrid functionals. In this work, we benchmark the random phase approximation (RPA) for ceria surface properties, including surface energy and hydrogenation energy, compared to the results utilizing the DFT + U approach or hybrid functionals. We show that, for the latter approaches, different surface properties require opposite directions of parameter tuning. This forms a dilemma for the parameter based DFT methods, as the improvement of a certain property by tuning parameters will inevitably lead to the worsening of other properties. Our results suggest that the parameter-free many-body perturbation theory methods exemplified by RPA are a promising strategy to escape the dilemma and provide highly accurate descriptions, which will allow us to better understand the catalytic reactions in ceria related systems. [ABSTRACT FROM AUTHOR]

We have witnessed considerable research interest in the recent literature about the development and applications of quantities from the information-theoretic approach (ITA) in density functional theory. These ITA quantities are explicit density functionals, whose local distributions in real space are continuous and well-behaved. In this work, we further develop ITA by systematically analyzing the topological behavior of its four representative quantities, Shannon entropy, two forms of Fisher information, and relative Shannon entropy (also called information gain or Kullback−Leibler divergence). Our results from their topological analyses for 103 molecular systems provide new insights into bonding interactions and physiochemical properties, such as electrophilicity, nucleophilicity, acidity, and aromaticity. We also compare our results with those from the electron density, electron localization function, localized orbital locator, and Laplacian functions. Our results offer a new methodological approach and practical tool for applications that are especially promising for elucidating chemical bonding and reactivity propensity. [ABSTRACT FROM AUTHOR]

Ryabinkin–Kohut–Staroverov (RKS) theory builds a bridge between wave function theory and density functional theory by using quantities from the former to produce accurate exchange-correlation potentials needed by the latter. In this work, the RKS method is developed and tested alongside Slater atomic orbital basis functions for the first time. To evaluate this approach, full configuration interaction computations in the Slater orbital basis are employed to give quality input to RKS, allowing full correlation to be present along with correct nuclei cusps and asymptotic decay of the wavefunction. SlaterRKS is shown to be an efficient algorithm to arrive at exchange-correlation potentials without unphysical artifacts in moderately-sized basis sets. Furthermore, enforcement of the nuclear cusp conditions will be shown to be vital for the success of the Slater-basis RKS method. Examples of weakly and strongly correlated molecular systems will demonstrate the main features of SlaterRKS. [ABSTRACT FROM AUTHOR]

Coarse-grained molecular dynamics (CGMD) simulations address lengthscales and timescales that are critical to many chemical and material applications. Nevertheless, contemporary CGMD modeling is relatively bespoke and there are no black-box CGMD methodologies available that could play a comparable role in discovery applications that density functional theory plays for electronic structure. This gap might be filled by machine learning (ML)-based CGMD potentials that simplify model development, but these methods are still in their early stages and have yet to demonstrate a significant advantage over existing physics-based CGMD methods. Here, we explore the potential of Δ-learning models to leverage the advantages of these two approaches. This is implemented by using ML-based potentials to learn the difference between the target CGMD variable and the predictions of physics-based potentials. The Δ-models are benchmarked against the baseline models in reproducing on-target and off-target atomistic properties as a function of CG resolution, mapping operator, and system topology. The Δ-models outperform the reference ML-only CGMD models in nearly all scenarios. In several cases, the ML-only models manage to minimize training errors while still producing qualitatively incorrect dynamics, which is corrected by the Δ-models. Given their negligible added cost, Δ-models provide essentially free gains over their ML-only counterparts. Nevertheless, an unexpected finding is that neither the Δ-learning models nor the ML-only models significantly outperform the elementary pairwise models in reproducing atomistic properties. This fundamental failure is attributed to the relatively large irreducible force errors associated with coarse-graining that produces little benefit from using more complex potentials. [ABSTRACT FROM AUTHOR]

Recent work on atomic-precision dopant incorporation technologies has led to the creation of both boron and aluminum δ -doped layers in silicon with densities above the solid solubility limit. We use density functional theory to predict the band structure and effective mass values of such δ layers, first modeling them as ordered supercells. Structural relaxation is found to have a significant impact on the impurity band energies and effective masses of the boron layers, but not the aluminum layers. However, disorder in the δ layers is found to lead to a significant flattening of the bands in both cases. We calculate the local density of states and doping potential for these δ -doped layers, demonstrating that their influence is highly localized with spatial extents at most 4 nm. We conclude that acceptor δ -doped layers exhibit different electronic structure features dependent on both the dopant atom and spatial ordering. This suggests prospects for controlling the electronic properties of these layers if the local details of the incorporation chemistry can be fine-tuned. [ABSTRACT FROM AUTHOR]

Ferroelectric oxides have gained research attention in the field of ferroelectric photovoltaics (PV) after the discovery of power conversion efficiency exceeding the Shockley–Queisser limit in BaTiO3 (BTO) crystals. However, advancement in this field is hindered by the wide bandgap (>3 eV) nature of ferroelectric oxides. In this work, a novel lead-free ferroelectric (1 − x)BTO − xBi(Ni2/3Nb1/3)O3 system was proposed and demonstrated to show bandgap reduction without compromising the polarization. Notably, the system displayed a bandgap reduction from 3.1 to 2.4 eV upon varying the composition from x = 0.0 to 0.05. Particularly, the optimal composition x = 0.02 showed enhancement in polarization (Pmax = 16 μC/cm2) and anomalous PV response with an open-circuit voltage of 6 V at 300 K. The origin of the bandgap reduction and polarization retention is explored experimentally by Raman spectroscopic measurements and analyzed theoretically by density functional theory. Our results revealed that the oxygen octahedral distortions and Ni2+ doping favor bandgap lowering, and Bi3+ ions stabilize the ferroelectric polarization. This study provides insight into the origin of bandgap tuning and paves the route for exploring new low-bandgap ferroelectric material with room temperature polarization. [ABSTRACT FROM AUTHOR]

Dirac cones are difficult to achieve in a square lattice with full symmetry. Here, we have theoretically investigated a bipartite tetragonal lattice composed of tetragons and octagons using both Tight-Binding (TB) model and density functional theory (DFT) calculations. The TB model predicts that the system exhibits nodal line semi-metallic properties when the on-site energies of all atoms are identical. When the on-site energies differ, the formation of an elliptical Dirac cone is predicted. Its physical properties (anisotropy, tilting, merging, and emerging) can be regulated by the hopping energies. An exact analytical formula is derived to determine the position of the Dirac point by the TB parameters, and a criterion for the existence of Dirac cones is obtained. The "divide-and-coupling" method is applied to understand the origin of the Dirac cone, which involves dividing the bands into several groups and examining the couplings among inter-groups and intra-groups. Various practical systems computed by DFT methods, e.g., t-BN, t-Si, 4,12,2-graphyne, and t-SiC, are also examined, and they all possess nodal lines or Dirac cones as predicted by the TB model. The results provide theoretical foundation for designing novel Dirac materials with tetragonal symmetry. [ABSTRACT FROM AUTHOR]

Machine learning potentials (MLPs) are poised to combine the accuracy of ab initio predictions with the computational efficiency of classical molecular dynamics (MD) simulation. While great progress has been made over the last two decades in developing MLPs, there is still much to be done to evaluate their model transferability and facilitate their development. In this work, we construct two deep potential (DP) models for liquid water near graphene surfaces, Model S and Model F, with the latter having more training data. A concurrent learning algorithm (DP-GEN) is adopted to explore the configurational space beyond the scope of conventional ab initio MD simulation. By examining the performance of Model S, we find that an accurate prediction of atomic force does not imply an accurate prediction of system energy. The deviation from the relative atomic force alone is insufficient to assess the accuracy of the DP models. Based on the performance of Model F, we propose that the relative magnitude of the model deviation and the corresponding root-mean-square error of the original test dataset, including energy and atomic force, can serve as an indicator for evaluating the accuracy of the model prediction for a given structure, which is particularly applicable for large systems where density functional theory calculations are infeasible. In addition to the prediction accuracy of the model described above, we also briefly discuss simulation stability and its relationship to the former. Both are important aspects in assessing the transferability of the MLP model. [ABSTRACT FROM AUTHOR]

Theoretical and experimental vibrational signatures of H2O and OH− (dissociated water) adsorbed on stoichiometric ceria{111} surfaces are compared. The experimental ones were collected from low-coverage experiments in the literature, and the theoretical anharmonic frequencies were generated using density functional theory calculations employing the optPBE-vdW functional for coverages from 0.5 to a few monolayers. It is found that (i) the experiments and our calculations overall agree well, lending credibility to both; (ii) the calculations manage to resolve the large class of H-bonded motifs into frequency classes that can guide experimental assignments; (iii) it is possible to find a geometrical H-bond definition that also captures the OH vibrational frequency downshifts well: R(H⋯O) ≤ 2.5 Å and the O–H⋯O angle θ ≥ 100°; and (iv) the frequency vs electric field relations for water and hydroxides (i.e., dissociated water) follow different and well-separated curves. [ABSTRACT FROM AUTHOR]

We study the interaction between heterogeneously charged surfaces in an electrolyte solution by employing classical Density Functional Theory (cDFT) and Monte Carlo simulations. We observe a consistent behavior between cDFT and Monte Carlo simulations regarding force curves and two-dimensional density profiles. Armed with the validated cDFT, we explore the system's behavior under parameters that are challenging to simulate directly. Our findings include the impacts of domain size, domain charge, domain charge configuration, and bulk electrolyte concentration on the osmotic pressure. Remarkably, the force curve is more sensitive to the domain size for an asymmetric configuration than a symmetry configuration; the bulk concentration weakly influences the force curve independent of the system configurations. [ABSTRACT FROM AUTHOR]

We have designed a [Fe(SH)4H]− model with the fifth proton binding either to Fe or S. We show that the energy difference between these two isomers (∆E) is hard to estimate with quantum-mechanical (QM) methods. For example, different density functional theory (DFT) methods give ∆E estimates that vary by almost 140 kJ/mol, mainly depending on the amount of exact Hartree–Fock included (0%–54%). The model is so small that it can be treated by many high-level QM methods, including coupled-cluster (CC) and multiconfigurational perturbation theory approaches. With extrapolated CC series (up to fully connected coupled-cluster calculations with singles, doubles, and triples) and semistochastic heat-bath configuration interaction methods, we obtain results that seem to be converged to full configuration interaction results within 5 kJ/mol. Our best result for ∆E is 101 kJ/mol. With this reference, we show that M06 and B3LYP-D3 give the best results among 35 DFT methods tested for this system. Brueckner doubles coupled cluster with perturbaitve triples seems to be the most accurate coupled-cluster approach with approximate triples. CCSD(T) with Kohn–Sham orbitals gives results within 4–11 kJ/mol of the extrapolated CC results, depending on the DFT method. Single-reference CC calculations seem to be reasonably accurate (giving an error of ∼5 kJ/mol compared to multireference methods), even if the D1 diagnostic is quite high (0.25) for one of the two isomers. [ABSTRACT FROM AUTHOR]

Rare earth sulfides are environmentally friendly and have excellent physical optoelectronic properties and coloring applications. In particular, red γ-Ce2S3 is an important and popular inorganic optical functional material. However, it has become crucial to improve the production efficiency of γ-Ce2S3 and its color performance and stability. By adding ammonium thiocyanate to the precursor, eco-friendly γ-Ce2S3 with excellent color performance was produced in this study. In addition, the effect of different amounts of Na+ doping on the final product's purity and the effect of chroma value were systematically investigated. Results show that an increase in Na+ can effectively lower the synthesis temperature of γ-Ce2S3, improve production efficiency, and optimize color performance prior to saturation. Finally, using density functional theory, the influence of the vacancy structure of γ-Ce2S3 and the change of energy band after the addition of Na+ were calculated, and a detailed experimental characterization in conjunction with density functional theory calculation thoroughly explained the significance of ion doping for reaction synthesis. [ABSTRACT FROM AUTHOR]

Borophene, the lightest member of mono-elemental 2D materials family, has attracted much attention due to its intriguing polymorphism. Among many polymorphs, digitally discovered 8-Pmmn stands out owing to its unique tilted-Dirac fermions. However, the property of interfaces between 8-Pmmn and metal substrates has so far remained unexplored, which has critical importance of its application in any electronic devices. Here, with the help of density functional theory, we show that the unique tilted-Dirac property is completely lost when 8-Pmmn borophene is interfaced with common electrode materials such as Au, Ag, and Ti. This is attributed to the high chemical reactivity of borophene as observed from crystal orbital Hamilton population and electron localization function analysis. In an effort to restore the Dirac property, we insert a graphene/hexagonal-boron-nitride (hBN) layer between 8-Pmmn and metal, a technique used in recent experiments for other 2D materials. We show that while the insertion of graphene successfully restores the Dirac nature for all three metals, hBN fails to do so while interfacing with Ti. The quantum chemical insights presented in this work may aid in to access the Dirac properties of 8-Pmmn in experiments. [ABSTRACT FROM AUTHOR]

The facet-dependent adsorption of CO on oxidized and reduced CeO2 single crystal surfaces is reviewed, with emphasis on the effect of CO coverage and the ability of state-of-the-art quantum-mechanical methods to provide reliable energies and an accurate description of the IR vibrational frequency of CO. Comparison with detailed, high-resolution experimental infrared reflection absorption spectroscopy data obtained for single crystal samples allows the assignment of the different CO vibrational bands observed on all three low-index ceria surfaces. Good agreement is achieved with the hybrid density functional theory approach with the HSE06 functional and with saturation coverage. It is shown that CO is very sensitive to the structure of cerium oxide surfaces and to the presence of oxygen vacancies. The combined theoretical-experimental approach offers new opportunities for a better characterization of ceria nanoparticles and for unraveling changes occurring during reactions involving CO at higher pressures. [ABSTRACT FROM AUTHOR]

The UV photochemistry of small heteroaromatic molecules serves as a testbed for understanding fundamental photo-induced chemical transformations in moderately complex compounds, including isomerization, ring-opening, and molecular dissociation. Here, a combined experimental-theoretical study of 268 nm UV light-induced dynamics in 2-iodothiophene (C4H3IS) is performed. The dynamics are experimentally monitored with a femtosecond extreme ultraviolet (XUV) probe that measures iodine N-edge 4d core-to-valence transitions. Experiments are complemented by density functional theory calculations of both the pump-pulse induced valence excitations and the XUV probe-induced core-to-valence transitions. Possible intramolecular relaxation dynamics are investigated by ab initio molecular dynamics simulations. Gradual absorption changes up to ∼0.5 to 1 ps after excitation are observed for both the parent molecular species and emerging iodine fragments, with the latter appearing with a characteristic rise time of 160 ± 30 fs. Comparison of spectral intensities and energies with the calculations identifies an iodine dissociation pathway initiated by a predominant π → π* excitation. In contrast, initial excitation to a nearby n⟂ → σ* state appears unlikely based on a significantly smaller oscillator strength and the absence of any corresponding XUV absorption signatures. Excitation to the π → π* state is followed by contraction of the C–I bond, enabling a nonadiabatic transition to a dissociative π → σ C − I * state. For the subsequent fragmentation, a relatively narrow bond-length region along the C–I stretch coordinate between 230 and 280 pm is identified, where the transition between the parent molecule and the thienyl radical + iodine atom products becomes prominent in the XUV spectrum due to rapid localization of two singly occupied molecular orbitals on the two fragments. [ABSTRACT FROM AUTHOR]

Understanding the structure and chemical bonding in water dimers is central to the study of many (photo-)electrochemical oxidation reactions. Two structures of the water dimer radical cation, namely, proton-transfer and hemi-bonded structures, have been suggested using density functional theory (DFT) and coupled cluster singles, doubles, and perturbative triples [CCSD(T)]. Both structures are identified by us as local minima, and their relative stability strongly depends on the level of theory. The exact exchange correlates linearly to the energy difference between both local minima. DFT functionals with less than 20 percent exact exchange predict the hemi-bonded structure to be more stable, while more than 20 percent of the exact exchange stabilizes the proton-transfer structure. The latter structure is also confirmed by CCSD(T) benchmark computations. These computations, furthermore, indicate that the oxidized water dimer consists of a hydronium cation (H3O+) and an HO· radical. These results are reproduced by DFT functionals with more than 50% of exact exchange (BHandH, M06-2X, and M06-HF). The transition barrier for the interconversion from the proton-transfer to the hemi-bonded structure is 0.6 eV, while the reverse reaction has a barrier of 0.1 eV. [ABSTRACT FROM AUTHOR]

Contemporary graph learning algorithms are not well-suited for large molecules since they do not consider the hierarchical interactions among the atoms, which are essential to determining the molecular properties of macromolecules. In this work, we propose Multiresolution Graph Transformers (MGT), the first graph transformer architecture that can learn to represent large molecules at multiple scales. MGT can learn to produce representations for the atoms and group them into meaningful functional groups or repeating units. We also introduce Wavelet Positional Encoding (WavePE), a new positional encoding method that can guarantee localization in both spectral and spatial domains. Our proposed model achieves competitive results on three macromolecule datasets consisting of polymers, peptides, and protein-ligand complexes, along with one drug-like molecule dataset. Significantly, our model outperforms other state-of-the-art methods and achieves chemical accuracy in estimating molecular properties (e.g., highest occupied molecular orbital, lowest unoccupied molecular orbital, and their gap) calculated by Density Functional Theory in the polymers dataset. Furthermore, the visualizations, including clustering results on macromolecules and low-dimensional spaces of their representations, demonstrate the capability of our methodology in learning to represent long-range and hierarchical structures. Our PyTorch implementation is publicly available at https://github.com/HySonLab/Multires-Graph-Transformer. [ABSTRACT FROM AUTHOR]

Quantum mechanical calculations for material modeling using Kohn–Sham density functional theory (DFT) involve the solution of a nonlinear eigenvalue problem for N smallest eigenvector-eigenvalue pairs, with N proportional to the number of electrons in the material system. These calculations are computationally demanding and have asymptotic cubic scaling complexity with the number of electrons. Large-scale matrix eigenvalue problems arising from the discretization of the Kohn–Sham DFT equations employing a systematically convergent basis traditionally rely on iterative orthogonal projection methods, which are shown to be computationally efficient and scalable on massively parallel computing architectures. However, as the size of the material system increases, these methods are known to incur dominant computational costs through the Rayleigh–Ritz projection step of the discretized Kohn–Sham Hamiltonian matrix and the subsequent subspace diagonalization of the projected matrix. This work explores the potential of polynomial expansion approaches based on recursive Fermi-operator expansion as an alternative to the subspace diagonalization of the projected Hamiltonian matrix to reduce the computational cost. Subsequently, we perform a detailed comparison of various recursive polynomial expansion approaches to the traditional approach of explicit diagonalization on both multi-node central processing unit and graphics processing unit architectures and assess their relative performance in terms of accuracy, computational efficiency, scaling behavior, and energy efficiency. [ABSTRACT FROM AUTHOR]

The effects of Mn and As doping in β-FeSi2 have been studied by theoretical simulations and electrical characterizations by analyzing Hall parameters within the temperature range of 20–300 K using mobility and the dual band model. The Hall resistivity ρ of doped samples increases linearly from a negative to a positive magnetic field (B), demonstrating the normal Hall effect at room temperature. High temperature Hall concentration increases significantly with the gradual increase in both Mn and As doping due to more and more ionization of the deep donor level. High temperature activation energies of Mn doped β-FeSi2 are considerably greater than that of low temperature energies, which demonstrates clear evidence of the dual band model. From density functional theory calculations, the origin of the dual band model has been validated from the electronic structure of β-FeSi2. Both density of states and charge transfer to the system upon doping have been investigated through the density functional theory, which demonstrates the Mn and As doped systems to be p-type and n-type, respectively. Both Mn and As doped β-FeSi2 exhibit p-type and n-type conductivities for spin down and spin up channels, respectively, in the presence of an external magnetic field, which will encourage its applications in novel spintronic devices. In addition, a β-FeSi2 based homo-junction diode fabricated from the Mn and As doped β-FeSi2 exhibits a cut-in voltage of 0.82 V, a reverse breakdown voltage of −10 V, and an ideality factor of 3.87. Thus, doped β-FeSi2 will be very much useful for fabricating an efficient and cost-effective solar cell if fabricated physically. [ABSTRACT FROM AUTHOR]

Amorphous titanium dioxide (a-TiO2) is widely used as a coating material in applications such as electrochemistry and self-cleaning surfaces where its interface with water has a central role. However, little is known about the structures of the a-TiO2 surface and aqueous interface, particularly at the microscopic level. In this work, we construct a model of the a-TiO2 surface via a cut-melt-and-quench procedure based on molecular dynamics simulations with deep neural network potentials (DPs) trained on density functional theory data. After interfacing the a-TiO2 surface with water, we investigate the structure and dynamics of the resulting system using a combination of DP-based molecular dynamics (DPMD) and ab initio molecular dynamics (AIMD) simulations. Both AIMD and DPMD simulations reveal that the distribution of water on the a-TiO2 surface lacks distinct layers normally found at the aqueous interface of crystalline TiO2, leading to an ∼10 times faster diffusion of water at the interface. Bridging hydroxyls (Ti2–ObH) resulting from water dissociation decay several times more slowly than terminal hydroxyls (Ti–OwH) due to fast Ti–OwH2 → Ti–OwH proton exchange events. These results provide a basis for a detailed understanding of the properties of a-TiO2 in electrochemical environments. Moreover, the procedure of generating the a-TiO2-interface employed here is generally applicable to studying the aqueous interfaces of amorphous metal oxides. [ABSTRACT FROM AUTHOR]

The equilibrium structure of selenophenol has been investigated using rotational spectroscopy and high-level quantum mechanical calculations, offering electronic and structural insight into the scarcely studied selenium compounds. The jet-cooled broadband microwave spectrum was measured in the 2–8 GHz cm-wave region using broadband (chirped-pulse) fast-passage techniques. Additional measurements up to 18 GHz used narrow-band impulse excitation. Spectral signatures were obtained for six isotopic species of selenium (80Se, 78Se, 76Se, 82Se, 77Se, and 74Se), together with different monosubstituted 13C species. The (unsplit) rotational transitions associated with the non-inverting μa-dipole selection rules could be partially reproduced with a semirigid rotor model. However, the internal rotation barrier of the selenol group splits the vibrational ground state into two subtorsional levels, doubling the dipole-inverting μb transitions. The simulation of the double-minimum internal rotation gives a very low barrier height (B3PW91: 42 cm−1), much smaller than for thiophenol (277 cm−1). A monodimensional Hamiltonian then predicts a huge vibrational separation of 72.2 GHz, justifying the non-observation of μb transitions in our frequency range. The experimental rotational parameters were compared with different MP2 and density functional theory calculations. The equilibrium structure was determined using several high-level ab initio calculations. A final Born–Oppenheimer ( r e BO ) structure was obtained at the coupled-cluster CCSD(T)_ae/cc-wCVTZ level of theory, including small corrections for the wCVTZ → wCVQZ basis set enlargement calculated at the MP2 level. The mass-dependent method with predicates was used to produce an alternative r m (2) structure. The comparison between the two methods confirms the high accuracy of the r e BO structure and offers information on other chalcogen-containing molecules. [ABSTRACT FROM AUTHOR]

Spin crossover (SCO) complexes, which exhibit changes in spin state in response to external stimuli, have applications in molecular electronics and are challenging materials for computational design. We curate a dataset of 95 Fe(II) SCO complexes (SCO-95) from the Cambridge Structural Database that have available low- and high-temperature crystal structures and, in most cases, confirmed experimental spin transition temperatures (T1/2). We study these complexes using density functional theory (DFT) with 30 functionals spanning across multiple rungs of "Jacob's ladder" to understand the effect of exchange–correlation functional on electronic and Gibbs free energies associated with spin crossover. We specifically assess the effect of varying the Hartree–Fock exchange fraction (aHF) in structures and properties within the B3LYP family of functionals. We identify three best-performing functionals, a modified version of B3LYP (aHF = 0.10), M06-L, and TPSSh, that accurately predict SCO behavior for the majority of the complexes. While M06-L performs well, MN15-L, a more recently developed Minnesota functional, fails to predict SCO behavior for all complexes, which could be the result of differences in datasets used for parametrization of M06-L and MN15-L and also the increased number of parameters for MN15-L. Contrary to observations from prior studies, double-hybrids with higher aHF values are found to strongly stabilize high-spin states and therefore exhibit poor performance in predicting SCO behavior. Computationally predicted T1/2 values are consistent among the three functionals but show limited correlation to experimentally reported T1/2 values. These failures are attributed to the lack of crystal packing effects and counter-anions in the DFT calculations that would be needed to account for phenomena such as hysteresis and two-step SCO behavior. The SCO-95 set thus presents opportunities for method development, both in terms of increasing model complexity and method fidelity. [ABSTRACT FROM AUTHOR]

The expensive cost of computing exact exchange in periodic systems limits the application range of density functional theory with hybrid functionals. To reduce the computational cost of exact change, we present a range-separated algorithm to compute electron repulsion integrals for Gaussian-type crystal basis. The algorithm splits the full-range Coulomb interactions into short-range and long-range parts, which are, respectively, computed in real and reciprocal space. This approach significantly reduces the overall computational cost, as integrals can be efficiently computed in both regions. The algorithm can efficiently handle large numbers of k points with limited central processing unit (CPU) and memory resources. As a demonstration, we performed an all-electron k-point Hartree–Fock calculation for LiH crystal with one million Gaussian basis functions, which was completed on a desktop computer in 1400 CPU hours. [ABSTRACT FROM AUTHOR]

Tailoring the specific properties to practical applications by structural modification is of vital importance for the envisioned development of two-dimensional ferroelectric materials. Herein, a comprehensive investigation on the effects of single doping on the ferroelectric properties and electronic transport in a monolayer of α-In2Se3 was carried out via the combination of first-principles density functional theory calculations and electron–phonon coupling simulations. Our results show that single-doping in In2Se3 can reduce effective mass of carriers and thereby enhance the high carrier mobility potential of the material. Moreover, the ferroelectric phonon mode in single-doped In2Se2X features a lower scattering rate, associating with the single-doping atom, and indicates reduced hindrance to carrier transport during ferroelectric switching. Compared to pristine In2Se3, the obtained smaller ferroelectric barriers (<1 eV) of single-doped ones promote more promising ferroelectricity from the analysis of the ferroelectric soft mode. Interestingly, the observed variations in ferroelectric behaviors resulting from doping of different elements highlight the significance of single-doping in modifying the ferroelectric properties of monolayers. Furthermore, strain engineering results reveal that single doping obviously affects the dependence of gap on strains: linear relationship for doping ones and nonlinearity for pristine one. Our study provides valuable insights into achieving higher carrier mobility in these critical materials. [ABSTRACT FROM AUTHOR]

We present a combined experimental and theoretical study of surface-related magnetic states in TiO 2. Our experiments on nano-sized thin films of pure TiO 2 have suggested that the observed room-temperature magnetism originates from defects, in particular, from the surface of thin films as well as from point defects, such as oxygen vacancies located mainly at the surface. Clarifying this phenomenon is very important for harnessing magnetic properties of pristine TiO 2 films in future spintronic applications but a detailed experimental investigation is very demanding. Therefore, quantum-mechanical density functional theory calculations were performed for (i) bulk anatase TiO 2 , (ii) bulk-like TiO 2 -terminated vacancy-free (001) surfaces, (iii) vacancy-containing TiO-terminated (001) surfaces, (iv) TiO 0.75 -terminated (001) surfaces with additional 25% surface oxygen vacancies, as well as (v) oxygen-terminated (001)-surfaces. Our fixed-spin-moment calculations identified both the bulk and the bulk-like terminated vacancy-free TiO 2 -terminated (001) surfaces as non-magnetic. In contrast, oxygen vacancies in the case of TiO-terminated and TiO 0.75 -terminated (001) surfaces lead to ferromagnetic and rather complex ferrimagnetic states, respectively. The spin-polarized atoms are the Ti atoms (due to the d -states) located in the surface and sub-surface atomic planes. Last, the O-terminated surfaces are also magnetic due to the surface and sub-surface oxygen atoms and sub-surface Ti atoms (but their surface energy is high). [ABSTRACT FROM AUTHOR]

We report a first-principles study on the electronic and optical properties of delafossite CuMO2 (M = Al, Ga and In) using the recently developed hybrid functional pseudopotentials. We obtain trends of the fundamental and optical gaps with increasing M-atomic number, in agreement with experiment. In particular, we reproduce the experimental fundamental gap, optical gap, and Cu 3d energy of CuAlO2 almost perfectly, in contrast to the various calculations that have traditionally focused on valence electrons, which are unable to reproduce these key properties simultaneously. Since all that distinguishes our calculations is simply the use of a different Cu pseudopotential with a partially exact exchange interaction, this suggests that an inappropriate description of the electron–ion interaction may play a role in the density functional theory bandgap problem for CuAlO2. Applying Cu hybrid pseudopotentials to CuGaO2 and CuInO2 is also effective, yielding optical gaps that are very close to experiment. However, due to the limited experimental data for these two oxides, a comprehensive comparison as that for CuAlO2 is not possible. Furthermore, our calculations yield large exciton binding energies for delafossite CuMO2, all around 1 eV. [ABSTRACT FROM AUTHOR]

Materials are needed to preferentially remove heavy metal pollutants from waste streams. Here, atomic scale simulation, based on density functional theory, is used to predict encapsulation energies and electronic structures of the one-dimensional yttrium silicide electride (Y5Si3:e–) hosting the heavy metals such as zinc (Zn), cadmium (Cd), mercury (Hg), chromium (Cr), nickel (Ni), copper (Cu), or lead (Pb). Assuming that the metals are present as dimers in gaseous waste, encapsulation is highly exoergic for those metals with incomplete p or d electrons (Cr, Ni, Cu, and Pb). Of the heavy metals with complete d shells (Zn, Cd, and Hg), only Zn exhibits strong encapsulation. Differences in encapsulation behavior are analyzed in terms of electronic structures, electron affinities, atomic radii, and Bader charges. [ABSTRACT FROM AUTHOR]

The research presented in the paper concerns the elastic properties of struvite. The article combines theoretical and experimental research. Experimental studies were carried out on struvite single crystals grown in sodium metasilicate gel by single diffusion. This unique method leads to obtaining crystals of sufficiently large size to conduct, for the first time, experimental measurement of elastic properties of monocrystalline struvite. Using the nanoindentation method, the Ez = 29.1 ± 0.7 GPa value of the component of Young's modulus was determined for a struvite single crystal. In addition, the elastic constants C11, C22, and C33 were determined using micro-Brillouin spectroscopy. Theoretical calculations of the abovementioned properties have been carried out by employing density functional theory methods. Scaling of the theoretical elastic constants leads to obtaining good agreement with the experimental values. Values of the Ex and Ey components of the Young's modulus, not available from the experimental nanoindentation technique, have been determined theoretically as 23 GPa and 27 GPa, respectively. Differences in the values of elastic components and Young's modulus components are related to the layered crystal structure of struvite and directional character of the hydrogen-bonding pattern. [ABSTRACT FROM AUTHOR]

Silver sulfide in monoclinic phase (α-Ag2S) has attracted significant attention owing to its metal-like ductility and promising thermoelectric properties near room temperature. However, first-principles studies on this material by density functional theory calculations have been challenging as both the symmetry and atomic structure of α-Ag2S predicted from such calculations are inconsistent with experimental findings. Here, we propose that a dynamical approach is imperative for correctly describing the structure of α-Ag2S. The approach is based on a combination of ab initio molecular dynamics simulation and deliberately chosen density functional considering both proper treatment of the van der Waals interaction and on-site Coulomb interaction. The obtained lattice parameters and atomic site occupations of α-Ag2S are in good agreement with experimental data. A stable phonon spectrum at room temperature can be obtained from this structure, which also yields a bandgap in accord with experimental measurements. The dynamical approach thus paves the way for studying this important ductile semiconductor in not only thermoelectric but also optoelectronic applications. [ABSTRACT FROM AUTHOR]

We report a Kohn–Sham density functional theory calculation of a system with more than 200 000 atoms and 800 000 electrons using a real-space high-order finite-difference method to investigate the electronic structure of large spherical silicon nanoclusters. Our system of choice was a 20 nm large spherical nanocluster with 202 617 silicon atoms and 13 836 hydrogen atoms used to passivate the dangling surface bonds. To speed up the convergence of the eigenspace, we utilized Chebyshev-filtered subspace iteration, and for sparse matrix–vector multiplications, we used blockwise Hilbert space-filling curves, implemented in the PARSEC code. For this calculation, we also replaced our orthonormalization + Rayleigh–Ritz step with a generalized eigenvalue problem step. We utilized all of the 8192 nodes (458 752 processors) on the Frontera machine at the Texas Advanced Computing Center. We achieved two Chebyshev-filtered subspace iterations, yielding a good approximation of the electronic density of states. Our work pushes the limits on the capabilities of the current electronic structure solvers to nearly 106 electrons and demonstrates the potential of the real-space approach to efficiently parallelize large calculations on modern high-performance computing platforms. [ABSTRACT FROM AUTHOR]

Superadiabatic dynamical density functional theory (superadiabatic-DDFT), a first-principles approach based on inhomogeneous two-body correlation functions, is employed to investigate the response of interacting Brownian particles to time-dependent external driving. Predictions for the superadiabatic dynamics of the one-body density are made directly from the underlying interparticle interactions without the need for either adjustable fit parameters or simulation input. The external potentials we investigate have been chosen to probe distinct aspects of structural relaxation in dense, strongly interacting liquid states. Nonequilibrium density profiles predicted by the superadiabatic theory are compared with those obtained from both adiabatic DDFT and event-driven Brownian dynamics simulation. Our findings show that superadiabatic-DDFT accurately predicts the time-evolution of the one-body density. [ABSTRACT FROM AUTHOR]

We report APW+lo (augmented plane wave plus local orbital) density functional theory (DFT) calculations of large molecular systems using the domain specific SIRIUS multi-functional DFT package. The APW and FLAPW (full potential linearized APW) task and data parallelism options and the advanced eigen-system solver provided by SIRIUS can be exploited for performance gains in ground state Kohn–Sham calculations on large systems. This approach is distinct from our prior use of SIRIUS as a library backend to another APW+lo or FLAPW code. We benchmark the code and demonstrate performance on several magnetic molecule and metal organic framework systems. We show that the SIRIUS package in itself is capable of handling systems as large as a several hundred atoms in the unit cell without having to make technical choices that result in the loss of accuracy with respect to that needed for the study of magnetic systems. [ABSTRACT FROM AUTHOR]

With the growing reliance of modern supercomputers on accelerator-based architecture such a graphics processing units (GPUs), the development and optimization of electronic structure methods to exploit these massively parallel resources has become a recent priority. While significant strides have been made in the development GPU accelerated, distributed memory algorithms for many modern electronic structure methods, the primary focus of GPU development for Gaussian basis atomic orbital methods has been for shared memory systems with only a handful of examples pursing massive parallelism. In the present work, we present a set of distributed memory algorithms for the evaluation of the Coulomb and exact exchange matrices for hybrid Kohn–Sham DFT with Gaussian basis sets via direct density-fitted (DF-J-Engine) and seminumerical (sn-K) methods, respectively. The absolute performance and strong scalability of the developed methods are demonstrated on systems ranging from a few hundred to over one thousand atoms using up to 128 NVIDIA A100 GPUs on the Perlmutter supercomputer. [ABSTRACT FROM AUTHOR]