Ab Initio Calculations Research Articles

Fast prediction of anharmonic vibrational spectra for complex organic molecules

Interpreting Raman and IR vibrational spectra in complex organic molecules lacking symmetries poses a formidable challenge. In this study, we propose an innovative approach for simulating vibrational spectra and attributing observed peaks to molecular motions, even when highly anharmonic, without the need for computationally expensive ab initio calculations. Our approach stems from the time-dependent stochastic self-consistent harmonic approximation to capture quantum nuclear fluctuations in atom dynamics while describing interatomic interaction through state-of-the-art reactive machine-learning force fields. Finally, we employ an isotropic charge model and a bond capacitor model trained on ab initio data to predict the intensity of IR and Raman signals.

Evaluating Sulfur as a P‐Type Dopant in Cu3N Using Ab Initio Methods

Copper nitride (Cu3N) is an environmentally friendly semiconducting material with bipolar doping capability and is of interest to various applications. As deposited Cu3N films have inherent n‐type conductivity, further controllable n‐type doping is possible by introducing metal impurities. First‐principles methods based on density functional theory and beyond have been employed to study the p‐type doping behavior of sulfur atoms in Cu3N. The structural, electronic, optical, and thermal properties of pure Cu3N and sulfur‐doped Cu3N are computed for single and 3 × 3 × 3 supercells. Sulfur doping causes a shift from intrinsic n‐type to p‐type behavior. This study confirms that sulfur atoms in sulfur‐doped copper nitride preferentially occupy interstitial positions over nitrogen substitution, face‐centered, or copper substitution sites. Due to this change and an increased lattice constant, Cu3N becomes a softer material with a larger bandgap in the single‐cell alloy. Doped Cu3N supercell results show significant changes in optical properties appropriate for solar and other photoelectric applications. Cu3N:S exhibits remarkable enhancements in power factor and thermal and electrical conductivity, indicating potentially better performance in thermoelectric applications. The dielectric constant and absorption coefficient also significantly change with the incorporation of sulfur into Cu3N.

Second-Order Mass-Weighting Scheme for Atom-Centered Density Matrix Propagation Molecular Dynamics.

The atom-centered density matrix propagation (ADMP) method is an extended Lagrangian approach to ab initio molecular dynamics, which includes the density matrix in an orthonormalized atom-centered Gaussian basis as additional, fictitious, electronic degrees of freedom, classically propagated along with the nuclear ones. A high adiabaticity between the nuclear and electronic subsystems is mandatory in order to keep the trajectory close to the Born-Oppenheimer (BO) surface. In this regard, the fictitious electronic mass μ, being a symmetric, nondiagonal matrix in its most general form, represents a free parameter, exploitable to optimize the propagation of the electronic density. Although mass-weighting schemes in ADMP exist, a systematic procedure to define an optimal value of the fictitious masses is not available yet. In this work, in order to rationally evaluate the electronic mass, fictitious electronic normal modes are defined through the diagonalization of the Hessian of the electronic density matrix. If the same frequency is imposed on all such modes (compatible with the chosen integration time step), then the corresponding μ matrix can be calculated and then employed for the following propagation. Analysis of several ADMP test simulations reveals that such Hessian-based mass-weighting approach is able to ensure, together with a 0.1/0.2 fs time steps, a high separation between the (real) nuclear and the (fictitious) electronic frequencies, which determines a high adiabaticity. This high, unprecedented, accuracy in the propagation leads, in turn, to low errors in the estimated nuclear vibrational frequencies, making the ADMP method totally comparable to a fully converged BO molecular dynamics simulation but more computationally efficient. This work, therefore, contributes to a further development of the ADMP ab initio molecular dynamics method, aimed at improving its accuracy through a more rational evaluation of the fictitious electronic mass parameter.

Theoretical investigations of two-dimensional intrinsic magnets derived from transition-metal borides M3B4 (M = Cr, Mn, and Fe)

ABSTRACT Two-dimensional (2D) magnetic materials with high critical temperatures (T C ) and robust magnetic anisotropy energies (MAE) hold significant potential for spintronic applications. However, most of 2D magnetic materials are derived from the van der Waals (vdW) layered bulks, which greatly limits the synthesis of 2D magnetic materials. Here, 2D M3B4 (M = Cr, Mn, and Fe; B = Boron), derived from hexagonal and orthorhombic M3AlB4 phases by selectively etching Al layers, was studied for its structural stability, electronic structure, and magnetic properties. By utilizing ab initio calculations and Monte Carlo simulations, we found that the orthorhombic Cr3B4 shows ferromagnetic (FM) metal and possesses an in-plane magnetic easy axis, while the remaining hexagonal and orthorhombic M3B4 structures exhibit antiferromagnetic (AFM) metals with a magnetic easy axis which is perpendicular to the two-dimensional plane. The critical temperatures of these 2D M3B4 structures are found to be above the 130 K. Notably, the ort-Mn3B4 possesses highest T C (~600 K) and strongest MAE (~220 µeV/atom) among these borides-based 2D magnetic materials. Our findings reveal that the 2D M3B4 compounds exhibit much better resistance to deformation compared to M2B2 MBenes and other 2D magnetic materials. The combination of high critical temperature, robust MAE, and excellent mechanical properties makes 2D Mn3B4 monolayer exhibits a favorable potential for spintronic applications. Our research also sheds light on the magnetic coupling mechanism of 2D M3B4, providing valuable insights into its fundamental characteristics.

Atomic and electronic properties of 2D Chevrel phases: A case study of the superatomic two-dimensional semiconductor Re6Se8Cl2

The design of two-dimensional superatomic materials, which form their atomic structures through covalently bonded clusters with variable chemical compositions, will enable the development of new materials with promised electronic properties that are beneficial for modern nanoelectronics. This paper presents ab initio calculations of the atomic and electronic structures of both bulk and 2D Re6Se8Cl2. The calculations were carried out using density functional theory, incorporating noncollinear spin density and the pseudopotential method. The results include data on the atomic structure, band gap value, formation energy of the Re6Se8Cl2 2D layer, and the redistribution of atomic charges within the structures. The differences in effective masses for electrons and holes in the two-dimensional and bulk Re6Se8Cl2 materials are demonstrated, along with an explanation of how these differences impact their transport properties. The findings are expected to be of great significance for the design, synthesis, and implementation of new two-dimensional superatomic materials with controlled properties in modern nanoelectronics.

Giant Valley Zeeman Splitting in Vanadium-Doped WSe2 Monolayers.

2D dilute magnetic semiconductors (DMS) based on transition metal dichalcogenides (TMD) offer an innovative pathway for advancing spintronic technologies, including the potential to exploit phenomena such as the valley Zeeman effect. However, the impact of magnetic ordering on the valley degeneracy breaking and on the enhancement of the optical transitions g-factors of these materials remains an open question. Here, a giant effective g-factors ranging between ≈-27 and -69 for the bound exciton at 4K in vanadium-doped WSe2 monolayers, obtained through magneto-photoluminescence (PL) experiments is reported. This giant g-factor disappears at room temperature, suggesting that this response is associated with a magnetic ordering of the vanadium impurity states at low temperatures. Ab initio calculations for the vanadium-doped WSe2 monolayer confirm the existence of magnetic ordering of the vanadium states, which leads to degeneracy breaking of the valence bands at K and K'. A phenomenological analysis is employed to correlate this splitting with the measured enhanced effective g-factor. The findings shed light on the potential of defect engineering of 2D materials for spintronicapplications.

Stimuli-Induced Fluorescence Switching in Azine-Containing Fluorophores Displaying Resonance-Stabilized ESIPT Emission.

This article reports the synthesis, along with structural and photophysical characterization of 2-(2'-hydroxyphenyl)benzazole derivatives functionalized with various azaheterocycles (pyridine, pyrimidine, terpyridine). These compounds show dual-state emission properties, that is intense fluorescence both in solution and in the solid-state with a range of fluorescent color going from blue to orange. Moreover, the nature of their excited state can be tuned by the presence of external stimuli such as protons or metal cations. In the absence of stimuli, these dyes show emission stemming from anionic species obtained after deprotonation (D* transition), whereas upon protonation or metal chelation, ESIPT process occurs leading to a stabilized and highly emissive K* transition. With the help of extensive ab initio calculations, we confirm that external stimuli can switch the nature of the transitions, making this series of dyes attractive candidates for the development of stimuli-responsive fluorescent ratiometric probes.

Simulation of lithium hydroxide decomposition using deep potential molecular dynamics.

Chemical reactions and vapor-liquid equilibria for molten lithium hydroxide (LiOH) were studied using molecular dynamics simulations and a deep potential (DP) model. The neural network for the model was trained on quantum density functional theory data for a range of conditions. The DP model allows simulations over timescales of hundreds of ns, which provide equilibrium compositions for the systems of interest. Single-phase NPT simulations of the liquid show the decomposition of LiOH into lithium oxide (Li2O) and dissolved water (H2O). These DP results were validated by direct abinitio molecular dynamics simulations that confirmed the accuracy of the model with respect to reaction kinetics and equilibrium properties of the melt. The reactive vapor-liquid behavior of this system was subsequently studied using direct coexistence interfacial DP simulations. Partial pressures of H2O in the vapor are found to be in close agreement with available experimental measurements. By fitting temperature-dependent expressions for the reaction equilibrium and Henry's law constants, the equilibrium composition for any given initial composition and temperature can be quantitatively modeled. For high initial concentrations of Li2O or H2O, mixtures of LiOH + Li2O/H2O are found to undergo phase separation. The present study illustrates how DP-based molecular dynamics simulations can be used for quantitative modeling of multiphase reactive behavior with the accuracy of the underlying abinitio quantum chemical methods.

PyCI: A Python-scriptable library for arbitrary determinant CI.

PyCI is a free and open-source Python library for setting up and running arbitrary determinant-driven configuration interaction (CI) computations, as well as their generalizations to cases where the coefficients of the determinant are nonlinear functions of optimizable parameters. PyCI also includes functionality for computing the residual correlation energy, along with the ability to compute spin-polarized one- and two-electron (transition) reduced density matrices. PyCI was originally intended to replace the abinitio quantum chemistry functionality in the HORTON library but emerged as a standalone research tool, primarily intended to aid in method development, while maintaining high performance so that it is suitable for practical calculations. To this end, PyCI is written in Python, adopting principles of modern software development, including comprehensive documentation, extensive testing, continuous integration/delivery protocols, and package management. Computationally intensive steps, notably operations related to generating Slater determinants and computing their expectation values, are delegated to low-level C++ code. This article marks the official release of the PyCI library, showcasing its functionality and scope.

Engaging the Concepts of Bimetallicity and Mechanical Strain for N2 Activation: A Computational Exploration.

N2 activation is a vital step in the process toward NH3 production. NH3 synthesis has been considered a crucial process for the production of value-added chemicals and/or hydrogen carriers over recent years. In this work, density functional theory (ab initio) calculations are implemented for a thorough screening of bimetallic alloy surfaces using Fe, Ru, and Mo as the matrix (host) metals and Ag, Au, Co, Cu, Fe, Mo, Ni, Pd, Pt, Rh, and Rh as heterometals toward exploring the N2 catalytic activation (electronic and chemical characteristics); the monometallic surfaces are used for critical comparison in terms of their N2 activation behavior. In particular, adsorption geometries/energetics, density of states (DOS), and charge transfer are discussed. From the N2 activation on the surfaces, we could precisely capture the transition state of the N2 dissociation reaction/step. The effect of the metal alloying (geometrical and electronic factors) as well as the effect of applied mechanical strain, as a tuning factor of alloying, are both studied and thoroughly discussed. DOS studies revealed that the d-band center moved toward the negative direction for all late-TM-based alloys, thereby allowing the nitrogen molecule to adsorb weakly as compared to the early-TM surface alloys. In terms of the mechanical strain, for most of the alloy surfaces studied, apart from the Mo/Fe(110) one, the N2 binding energy varies as a linear function of the applied strain. The mechanical effect trend is in agreement with the charge transfer descending order followed: Fe/Mo(110) > Rh/Mo(110) > Au/Mo(110) > Pt/Mo(110) > Ni/Mo(110) > Ru/Mo(110) > Cu/Mo(110) > Ag/Mo(110) > Pd/Mo(110) > Au/Mo(110), pointing out that Fe-functionalized Mo(110) surface presents the highest charge transfer of -2.14 |e| to the N2 molecule. This study aspires to provide navigation criteria through the abundant design criteria of N2 activation catalysts.

Tailoring electronic structure and thermodynamic stability of (Al, In)-substituted GaAs: Ab-initio insights into bulk and (001) surfaces

Ab-initio calculations are conducted to investigate the structural properties, electronic structure, and thermodynamic stabilities of bulk and (001) surface systems of the GaAs semiconductor in its zinc-blende crystal phase. For the bulk case, we focus on the impact of substituting gallium atoms with M atoms (M = Al, In) at various concentrations (5.55 %, 11 %, and 22 %). We present the resulting structural parameters, along with the formation and substitutional energies associated with these substitutions. For surface systems, we explore indium (In) or aluminum (Al) substitutions at coverages of 0.25, 0.50, and 1.0 monolayers. This includes investigating the effects of substituting M atoms within the first atomic layer of the slabs. Ab-initio thermodynamic calculations reveal that high Al substitutions, in both bulk and surface, result in more favorable formation and surface energies. In contrast, In substitutions exhibit the opposite behavior, with increasing In content leading to less favorable energetics. Moreover, density of states analysis reveals that bulk systems maintained a semiconducting character, while all studied surface configurations exhibited a metallic electronic behavior. This study provides valuable insights into the influence of Al and In substitutions on the structural, electronic, and thermodynamic properties of GaAs, both in bulk and on the (001) surface.

The suprachelate effect: a pathway to supramolecules with Cu(II) and ethylenediamine

ABSTRACT This paper portrays the solution-phase dynamics as copper(II) and ethylenediamine explore a multitude of different complexes. The data from five spectrophotometric titrations were globally analysed, evidencing four predominant species ([Cu]2+, [Cu(en-N,N’)]2+, [Cu(en-N,N’)2]2+, [Cu(en-N)4]2+) along with their molar absorptivity curves and associative binding constants. The data also seem to support a fifth species, [Cu2(µ-en-N,N’)]4+, in which ethylenediamine bridges two Cu(II) centres. The thermodynamic stability of all five species is corroborated by ab initio computational calculations. The potential existence of [Cu(en-N)4]2+ highlights the suprachelate effect – going beyond the chelate effect – where multidenticity is overtaken by monodenticity. Such dangling multidentate ligands are available to bind to additional metal centres and thus build towards self-assembling supramolecules.

Signature of pressure-induced topological phase transition in ZrTe5

The layered van der Waals material ZrTe5 is known as a candidate topological insulator (TI), however its topological phase and the relation with other properties such as an apparent Dirac semimetallic state is still a subject of debate. We employ a semiclassical multicarrier transport (MCT) model to analyze the magnetotransport of ZrTe5 nanodevices at hydrostatic pressures up to 2 GPa. The temperature dependence of the MCT results between 10 and 300 K is assessed in the context of thermal activation, and we obtain the positions of conduction and valence band edges in the vicinity of the chemical potential. We find evidence of the closing and re-opening of the band gap with increasing pressure, which is consistent with a phase transition from weak to strong TI. This matches expectations from ab initio band structure calculations, as well as previous observations that CVT-grown ZrTe5 is a weak TI in ambient conditions.

The Degradation of the 2-Hydroxyethylhydrazinium Nitrate Ionic Liquid System.

The mechanism by which the 2-hydroxyethylhydrazinium nitrate (HEHN) ionic liquid (IL) dissociates via the formation and loss of nitric acid is investigated utilizing the effective fragment potential (EFP) method and various ab initio methods. Additionally, the interactions that dictate the stability of the HEHN IL and contribute to the formation of nitric acid are further assessed using the quasi-atomic orbital (QUAO) bonding analysis. From this work, it is suggested that the formation of nitric acid is dictated by the presence of hydrogen-bonding interactions which appear to become increasingly ionic in nature as the system approaches the bulk, limiting the formation of nitric acid.

Exploring the Mechanisms of LiNiO2 Cathode Degradation by the Electrolyte Interfacial Deprotonation Reaction.

Nickel-rich layered oxides stand as ideal cathode candidates for high specific capacity and energy density next-generation lithium-ion batteries. However, increasing the Ni content significantly exacerbates structural degradation under high operating voltage, which greatly restricts large-scale commercialization. While strategies are being developed to improve cathode material stability, little is known about the effects of electrolyte-electrode interaction on the structural changes of cathode materials. Here, using LiNiO2 in contact with electrolytes with different proton-generating levels as model systems, we present a holistic picture of proton-induced structural degradation of LiNiO2. Through ab initio molecular dynamics calculations based on density functional theory, we investigated the mechanisms of electrolyte deprotonation, protonation-induced Ni dissolution, and cathode degradation and the impacts of dissolved Ni on the Li metal anode surfaces. We show that the proton-transfer reaction from electrolytes to cathode surfaces leads to dissolution of Ni cations in the form of NiOOHx, which stimulates cation mixing and oxygen loss in the lattice accelerating its layered-spinel-rock-salt phase transition. Migration of dissolved Ni2+ ions to the anode side causes their reduction into the metallic state and surface deposition. This work reveals that interactions between the electrolyte and cathode that result in protonation can be a dominant factor for the structural stability of Ni-rich cathodes. Considering this factor in electrolyte design should be of benefit for the development of future batteries.

Study of the cross sections for the resonant nuclear reactions induced by 6Li + 2H collisions based on ab initio computations of characteristics of highly excited states of 8Be spectrum

Study of the cross sections for the resonant nuclear reactions induced by ⁶Li + ²H collisions based on ab initio computations of characteristics of highly excited states of ⁸Be spectrum

Cooperativity and halonium transfer in the ternary NCI···CH3I···-CN halogen-bonded complex: An ab initio gas phase study.

The strength and nature of the two halogen bonds in the NCI···CH3I···-CN halogen-bonded ternary complex are studied in the gas phase via ab initio calculations. Different indicators of halogen bond strength were employed to examine the interactions including geometries, complexation energies, Natural Bond Order (NBO) Wiberg bond indices, and Atoms in Molecules (AIM)-based charge density topological properties. The results show that the halogen bond is strong and partly covalent in nature when CH3I donates the halogen bond, but weak and noncovalent in nature when CH3I accepts the halogen bond. Significant halogen bond cooperativity emerges in the ternary complex relative to the corresponding heterodimer complexes, NCI···CH3I and CH3I···-CN, respectively. For example, the CCSD(T) complexation energy of the ternary complex (-18.27kcal/mol) is about twice the sum of the complexation energies of the component dimers (-9.54kcal/mol). The halonium transfer reaction that converts the ternary complex into an equivalent one was also investigated. The electronic barrier for the halonium transfer was calculated to be 6.70kcal/mol at the CCSD(T) level. Although the MP2 level underestimates and the MP3 overestimates the barrier, their calculated MP2.5 average barrier (6.44kcal/mol) is close to that of the more robust CCSD(T) level. Insights on the halonium ion transfer reaction was obtained by examining the reaction energy and force profiles along the intrinsic reaction coordinate, IRC. The corresponding evolution of other properties such as bond lengths, Wiberg bond indices, and Mulliken charges provides specific insight on the extent of structural rearrangements and electronic redistribution throughout the entire IRC space. The MP2 method was used for geometry optimizations. Energy calculations were performed using the CCSD(T) method. The aug-cc-pVTZ basis set was employed for all atoms other than iodine for which the aug-cc-pVTZ-PP basis set was used instead.

Generation of femtosecond spin-polarized current pulses at Fe/MgO interface by quasi-static voltage

The generation of short spin-current pulses is essential for fast spintronic devices. So far, spin current pulses are generated by femtosecond laser pulses which drive spins from a ferromagnetic metal layer. However, the need for miniaturization, simplicity and energy efficiency favour electric-field control of spintronic devices over optic or thermal control. Here, we combine ab initio calculations of electronic density of states at MgO/Fe interface with continuous model for charge transport to investigate the dynamics of the spin-dependent potential. We demonstrate that the voltage-driven instability of the electronic band structure due to the electronic resonant states at the Fe/MgO interface results in the generation of the femtosecond spin-polarized current pulse with the spin polarization up to P=7 00 % that propagates from the interface to the bulk. The dynamics of the current pulses driven by the Stoner instability depends neither on the dielectric relaxation time nor on the details of how the instability is achieved by changing the voltage, i.e. as long as the voltage changes are slow (quasi-static) with respect to the time determined by the spin diffusion constant, being of the order of fs. The presence of the instability can be detected by THz time-domain spectroscopy or pump-probe techniques.

Structural, optoelectronic, and magnetic properties of Q‑carbon studied by hybrid density functional theory ab initio calculations and experiment

This work presents a theoretical study of the structural, optoelectronic, and magnetic properties of Q‑carbon, utilizing hybrid functionals within density functional theory (DFT) ab initio calculations. Moreover, experiments were conducted to measure structural parameters, dielectric functions, and magnetic properties, using various techniques. From the DFT simulations, structural, electronic, and optical properties have been simulated. The band gap energy was calculated using a hybrid approach, which combine HSE functional with different values of exact Hartree-Fock (HF) exchange (α). We propose a theoretical investigation of cubic and quasi-cubic forms of Q‑carbon. Our findings indicate that a non-cubic Q‑carbon cell is energetically more stable and possesses higher cell mass density and bulk modulus compared to a cubic cell. Calculations indicated that Q‑carbon exhibits semiconductor behavior; an indirect band gap of 3.10 eV and a direct band gap of about 3.4 eV was obtained for α = 0.44, in good agreement with experimental values. The density of states analysis allowed us to identify the p orbital of sp2 hybridized atoms as being prevalent on the top of the valence band and in the lowest energy conduction band states. Magnetic moment values of 0.37 μB and 0.41μB was found in two sp2 bonded C atoms, in accordance with the experimental observation. Several optical properties were calculated. The theoretical findings are compatible with the experimental results shown here and available in the literature, with excellent agreement found between the calculated optical band gap and that obtained from absorption measurements.

Construction of Al–Si interatomic potential based on Bayesian active learning

Nanoscale simulations for optimizing the performance and processing of Al–Si alloys are currently facing two major obstacles: the scarcity of high-quality semi-empirical potentials tailored to complex alloy systems, and the prohibitively high computational cost associated with ab initio molecular dynamics simulations. In order to enhance simulation efficiency and accuracy of the Al–Si alloys’ microstructural evolution, this study employs a dynamic active learning technique, FLARE, to develop a non-parametric machine learning potential that combines the high accuracy of density functional theory (DFT) with the efficiency of classical molecular dynamics (MD). Without relying on extensive initial ab initio molecular dynamics data or existing databases, collection of the necessary data is progressively made during the active learning process, thereby constructing a potential capable of accurately simulating the structure and dynamics of high-temperature Al–Si alloys. By comparing with experimental measurements and ab initio molecular dynamics calculations, the high accuracy and computational efficiency of this potential is demonstrated in predicting energy, force, structure, and dynamic properties. The results provide novel theoretical insights for optimizing Al–Si alloy processing and underscore the usefulness of active learning methods in constructing high-accuracy potentials.