Ab Initio Calculations Research Articles

Influence of atomic substitution on the structural stability and half-metallicity of Fe2-xCrxCoSi (x = 0 to 1) alloys

Bulk Fe2-xCrxCoSi (x = 0, 0.25, 0.5, 0.75, and 1) Heusler alloys have been prepared by arc melting method. Alloys with 0.0 ≤ x ≤ 0.75 show a highly ordered phase pure XA structure. However, 17 % of an impurity (A15 phase) was detected in the alloy with x = 1. The saturation magnetization (Ms) and Curie temperature (TC) decreased linearly with an increase in Cr content for alloys with x ≤ 0.75. Ms measured at 5 K followed the Slater–Pauling rule for half-metals. Ab initio calculations with the GGA + U approach revealed that spin polarization (P) increased with an increase in Cr and eventually reached 100 % for the alloys with x = 0.5. Though Fe-Co disorder affects P, it is still high (>98 %) for alloys with x ≥ 0.5. Despite the smaller bandgap (0.61 eV) of Fe1.5Cr0.5CoSi than FeCrCoSi (0.91 eV), it has higher Ms, elevated TC, and 100 % P. These factors, coupled with the challenges associated with obtaining phase pure FeCrCoSi alloy, make Fe1.5Cr0.5CoSi a preferred candidate for spintronic device applications.

Electronic and Magnetic Properties of Cobalt Clusters on Pristine and Divacancy Graphene

Using ab initio calculations we investigate the adsorption of Co atoms, dimers and small cobalt clusters of 5 and 13 atoms on pristine graphene and graphene with a double vacancy. We report the atomic, electronic, magnetic and energetic properties of these systems. Stable adsorption configurations tend to maximize the number of cobalt-carbon bonds. On graphene, the adsorption energy of the clusters is only about 0.4 to 1 eV, and the clusters are relatively mobile on graphene. Interestingly, for different adsorbed Co13 isomers on graphene it is found that they converge to the same atomic structure. On graphene with a divacancy, the Co clusters bind in the divacancy site and isomerisation also occurs for the Co5 cluster system as well as for Co13. Co atoms and clusters can be effectively immobilized on the divacancy with corresponding adsorption energy being significantly enhanced by about 5 to 7 eV. All clusters act as electron donors in the interaction with the graphene/divacancy systems, and the amount of electron charge transfer increases with cluster size. Finite magnetic moments occur for all systems, where upon adsorption, the magnetic moment of the isolated Co atom (3μB) is significantly reduced due to electron transfer and bonding, resulting in values varying from ≈0.9-2.2 μB per Co atom. For the pristine graphene substrate, the total induced magnetic moments on the carbon atoms are negligible, while on the divacancy system, they are of the order of 0.1-0.3 μB. The attractive physical properties of these hybrid systems could find applications in catalysis and materials science.

Vacuum-thermal alteration of lunar soil: Evidence from iron whiskers on troilite in Chang’e-5 samples

The formation of a unique microstructure of minerals on the surface of airless bodies is attributed to space weathering. However, it is difficult to distinguish the contributions of meteorite impacts and solar wind to the modification of lunar soil, resulting in limited research on the space weathering mechanism of airless bodies. The thermochemical reactivity of troilite can be used to distinguish the contributions of impact events and solar wind to the modification of lunar soil and provide evidence for space weathering of lunar soil. We examined the structure of a single particle of troilite in the Chang’e-5 lunar soil and determined whether an impact caused the thermal reaction. Microanalysis showed that troilite underwent substantial mass loss during thermal desulfurization, forming a crystallographically aligned porous structure with iron whiskers, an oxygen-rich layer, and other crystallographic and thermochemical evidence. We used an ab initio deep neural network model and thermodynamic calculations to conduct experiments and determine the anisotropy and crystal growth of troilite. The surface microstructure of troilite was transformed by the thermal reaction in the vacuum on the lunar surface. Similar structures have been found in near-Earth objects (NEOs), indicating that small bodies underwent the same impact-induced thermal events. Thus, thermal reactions in a vacuum are likely ubiquitous in the solar system and critical for space weathering alterations of the soil of airless bodies.

Modelling of (N-vinylcarbazole)/Fullerene nanoheterojunction for organic solar cells and photovoltaics applications – A DFT Approach

The current manuscript depicts the construction of a novel hybrid carbon nanomaterials with fascinating electronic and chemical properties and tunable energy bandgap of a nano heterojunction comprising (9-vinyl carbazole), VK and fullerene, C60. A comprehensive investigation of VK/C60 adsorption surface is conducted using ab-initio density functional theory calculations with van der Waals corrections (GGA+vdW) implemented. The novel nanoheterostructure interfaces synthesized by DFT methods have low-bandgap characteristics (1.112–1.223 eV), broaden absorption wavelengths and photovoltaic characteristics. This novel VK/C60 nanoheterostructure composites have finite bandgaps and good stability and, therefore, have promising applications for photovoltaic devices. Optoelectronic features of all three developed systems are recommended for future solar cell applications. The studies reveal that the VK/C60 nanoheterojunctions can behave as efficient electron donor-acceptor components in organic solar cells (OSCs) that make them potential candidates for the development of low-cost optoelectronic devices. The band gap values calculated in this work are low enough to make the nanoheterostructures exceedingly promising for photovoltaic applications.

Phase equilibrium and separation process of byproduct vinyl acrylate in the ethylene vapor phase method producing vinyl acetate

Vinyl acrylate is one of the byproducts in the ethylene vapor phase method for the production of vinyl acetate (VAc), which not only has a toxic effect on the catalyst but also affects the purity of the product. Currently, phase equilibrium data for vinyl acrylate with water are missing, making the destination of vinyl acrylate in the distillation process uncertain. This has led to a lack of effective removal means, which restricts the optimization of the ethylene vapor phase method. This paper uses a method combining quantum chemistry and molecular dynamics simulation (MD) to investigate the existential form of the water molecule. And based on the ab initio calculation of the first principle, the force field parameters of the special water molecular structure are calculated, and the TIP4P-P force field that can describe the special water structure is established. Using the Gibbs Ensemble Monte Carlo (GEMC) method, the reliability of the improved water molecule force field for modeling the phase equilibrium of a binary system is first verified. Further, vapor–liquid phase equilibrium data for the vinyl acrylate-water binary system at 110 kPa and 150 kPa are calculated and the corresponding x-y and T-x-y phase diagrams are plotted. After performing the thermodynamic consistency test, the phase equilibrium data are fitted to obtain the binary interaction parameters of vinyl acrylate with water. The fitted parameters are applied to the process simulation of the VAc refining process to optimize the conditions and obtain an effective separation scheme for vinyl acrylate. The concentration distribution of vinyl acrylate in the distillation column is investigated, thus clarifying the destination of vinyl acrylate in the distillation process and facilitating the optimization of the ethylene vapor phase method.

Optoelectronic and Transport Properties of Layer-dependent Two-Dimensional Perovskite Cs3Bi2I9

The emergence of two-dimensional (2D) perovskite provides an ideal platform for designing and fabricating microelectronic and optoelectronic devices. Here, we present detailed ab initio calculations to comprehensively explore the layer-dependent optoelectronic and transport properties of the newly synthesized 2D perovskite Cs3Bi2I9. The calculations reveal that the indirect band gap of monolayer, bilayer, and trilayer Cs3Bi2I9 decreases from 2.41 to 2.35 eV. The charge carriers of few-layered perovskite Cs3Bi2I9 are dominated by electrons, and the highest electron and hole carrier mobilities are 140 and 62 cm²V⁻¹s⁻¹ along the b axis for bilayer Cs3Bi2I9. Exciton binding energies decrease from 2.48 to 1.19 eV with an increment of layer, and the calculated exciton level drops down into the valence band to generate potential exciton insulator. 2D Cs3Bi2I9 exhibits potential in the field of ultraviolet detection and photoluminescent devices due to large exciton energy and ultraviolet absorption.

Electron-Spin Relaxation in Boron-Doped Graphene Nanoribbons.

Boron-doped graphene nanoribbons are promising platforms for developing organic materials with magnetic properties. Boron dopants can be used to create localized magnetic states in nanoribbons with tunable interactions. Controlling the coherence times of these magnetic states is the very first step in designing materials for quantum computation or information storage. In this work, we address the connection between the relaxation time and the position of the dopants for a series of boron-doped graphene nanofragments. We combine Redfield theory and ab initio calculations of magnetic properties to unveil the mechanism that governs spin relaxation in solution. We demonstrate that relaxation times can be in the order of 1 ms for the selected graphene nanofragments. A detailed analysis of the relaxation mechanism reveals that the spin decoherence is fundamentally driven by fluctuations of the spin-orbit coupling, and the hyperfine interaction facilitated by the thermal motion of the graphene nanofragments. The close connection between relaxation time, hyperfine interaction and the spin-orbit coupling offers the perspective of designing attractive materials with long-lived spin states.

Harnessing Free Sulfate Groups in Glycosylation Reactions.

General strategies for synthesizing sulfated oligosaccharides employ the protection of sulfate groups for glycosylation or post-sulfation after the synthesis of oligosaccharides. However, the chemical behavior of free sulfate groups in glycosylation reactions has not been thoroughly studied. We examined several glycosyl donors with free sulfate groups, but neither glycosyl imidates nor thioglycosides achieved products. Conversely, activating 6-sulfated GalNAc donors with either a 2-NTroc group or 2-O-acetyl group using an ortho-hexynylbenzoate group at the anomeric position yielded β-glycosides in good yield. These results indicate that glycosylations with a free sulfate group can be performed and that the neighboring group participation of 2-NHTroc and 2-O-acetyl groups works even with an unprotected sulfate group at 6-position. Ab initio calculations supported the formation of acyloxionium-cation via 2-O-acyl participation. Additionally, glycosylation reactions under various counter cations of the sulfate group, such as Na+, Li+, K+, Ba2+ and pyridinium were examined. Results showed that sodium and lithium salt donors yielded the products in good yield. These results were also supported by ab initio calculations. Practical glycosylation reactions between disaccharide donor-acceptor pairs with free sulfate groups successfully yielded the sulfated tetrasaccharides. This study also discusses how a sodium salt acts as a protecting group during glycosylation.

Synergistic Catalytic Sites in High-Entropy Metal Hydroxide Organic Framework for Oxygen Evolution Reaction.

The integration of multiple elements in a high-entropy state is crucial in the design of high-performance, durable electrocatalysts. High-entropy metal hydroxide organic frameworks (HE-MHOFs) are synthesized under mild solvothermal conditions. This novel crystalline metal-organic framework (MOF) features a random, homogeneous distribution of cations within high-entropy hydroxide layers. HE-MHOF exhibits excellent electrocatalytic performance for the oxygen evolution reaction (OER), reaching a current density of 100mA cm-2 at ≈1.64 VRHE, and demonstrates remarkable durability, maintaining a current density of 10mA cm-2 for over 100 h. Notably, HE-MHOF outperforms precious metal-based electrocatalysts despite containing only ≈60% OER active metals. Ab initio calculations and operando X-ray absorption spectroscopy (XAS) demonstrate that the high-entropy catalyst contains active sites that facilitate a multifaceted OER mechanism. This study highlights the benefits of high-entropy MOFs in developing noble metal-free electrocatalysts, reducing reliance on precious metals, lowering metal loading (especially for Ni, Co, and Mn), and ultimately reducing costs for sustainable water electrolysis technologies.

Efficient Selective CO2 Composite Sorbent from Amino Acid Ionic Liquids and Silica Gel: 2H NMR Spectroscopy Provides Insight on the CO2-Binding Mechanism and in-Pore Microscopic Viscosity.

A promising CO2 sorbent based on the [Emim][Gly] (1-ethyl-3-methylimidazolium glycinate)/silica gel composite has been studied. CO2 sorption experiments have shown that the optimum loading of [Emim][Gly] ionic liquid in silica gel is 40 wt.%. The 2-step process CO2 binding mechanism in [Emim][Gly] has been proposed based on the results of sorption experiments, 2H NMR spectroscopy and ab-initio calculations. The impact of CO2 on the microscopic viscosity and the dynamical melting of ionic liquid has been thoroughly investigated. 2H NMR spectroscopy has revealed that CO2 strongly binds cation and anion in [Emim][Gly], forcing them to move in a correlated fashion.

Electronic structure of 2(5H)-thiophenone studied by UPS and soft X-ray spectroscopy

Abstract The electronic structure of 2(5H)-thiophenone in the gas phase was investigated by ultraviolet photoelectron spectroscopy (UPS) and x-ray photoemission spectroscopy (XPS) and near edge x-ray absorption fine structure (NEXAFS) spectroscopy at the C 1s, O 1s, and S 2p edges. All assignments of the experimental results are supported by both ab-initio electron propagator OVGF calculations for the valence photoemission bands and density functional theory (DFT) and relativistic time dependent DFT calculations for the core levels XPS and NEXAFS spectra. Overall good agreement between experiment and theory is observed; this is especially true for core electron excitations which has permitted an unambiguous assignment of the observed spectral features in terms of single-particle excitations to virtual molecular orbitals. The assignment of the valence band spectra based on OVGF calculations, although satisfactory, points to the importance of electron correlations effects that partially break the single particle picture of ionization.

Understanding the Stability of Finite Double Walled Phenine and its Hybrid Structures using DFT Investigation

Abstract We have performed ab initio calculations based on density functional theory to investigate the structural and electronic properties
of finite double wall phenine nanotubes and their hybrid structures with armchair carbon nanotubes. The result shows that double wall
phenine nanotubes and the hybrid structures are quite stable. The stability of structure is strongly depends on inter-layer separation be-
tween the two tubes. The double walled phenine nanotubes are more stable than their corresponding pristine armchair tubes and exhibit
a semi-conducting nature with a HOMO-LUMO gap of 2.78eV and 2.80eV. The hybrid structures have a small electronic gap between
0.13eV-0.50eV similar to double wall armchair nanotubes. The results highlight the possibility of the synthesis of novel semi-conducting
double wall phenine structures.

Cluster Perturbation Theory for Core Excited States and Core Ionization Potentials Using Core-Valence Separation.

The development of accurate and fast computational procedures for the ab initio calculation of X-ray spectroscopies is paramount to facilitate theoretical analysis of modern X-ray experiments on molecules. Herein, we present the extension of Cluster Perturbation theory to comprehend the calculation of core excited states and core ionization potentials using the core-valence separation approximation, which has seen widespread success for various quantum chemistry methods. We derive the theoretical framework for introducing core-valence separation into Cluster Perturbation series for excitation energies and display the performance of the methodology in S(D) orbital excitation spaces. The obtained core excitation energies on a test set of medium sized organic molecules show that carbon, nitrogen, and oxygen K-edge excitation energies can be determined with errors below 2 eV relative to the CCSD reference results using the developed CPS(D) excitation energy models which can be used for systems way beyond the reach of conventional CCSD.

Pressure-induced structural stability, metallization and optoelectronic properties of transition metal chalcogenide BaTiS3

The pressure effects on the structural stability, metallization and optoelectronic properties of transition metal chalcogenide BaTiS3 are investigated by performing ab-initio calculations. Our calculations show that BaTiS3 has a hexagonal structure at ambient pressure, which matches well with the available experimental values. As pressure increases, a transition from hexagonal to orthorhombic phase occurs at around 10.8 GPa. The band structure and density of states of hexagonal BaTiS3 have been calculated and analysed. Moreover, the effect of pressure on the optical properties has been studied in terms of dielectric function as well as absorption coefficient, reflectivity and extinction coefficient. The application of hydrostatic pressure has shifted the maximum absorption toward the visible range, revealing that hexagonal BaTiS3 can be used for high pressure optoelectronic applications.

Trigonal Planar Heteroleptic Lanthanide(III) Bis(silyl)amide Complexes Containing Aminoxyl Radicals and Anions.

Modulation of the crystal field (CF) in lanthanide (Ln) complexes can enhance optical and magnetic properties, and large CF splitting can be achieved with low coordination numbers in specific geometries. We previously reported that the homoleptic near-linear Sm2+ complex [SmII{N(SiiPr3)2}2] (1-Sm) is oxidized by the 2,2,6,6-tetramethylpiperidinyl-1-oxy (TEMPO•) radical to give the heteroleptic, approximately trigonal planar Sm3+ complex, [SmIII{N(SiiPr3)2}2(TEMPO-)] (2-Sm). Here, we report the synthesis of homologous [LnIII{N(SiiPr3)2}2(TEMPO-)] (2-Ln; Ln = Tm, Yb) complexes by the oxidation of the parent [Ln{N(SiiPr3)2}2] (1-Ln; Ln = Tm, Yb) with TEMPO•; complexes 2-Ln all contain TEMPO- anions. The homoleptic bent Ln3+ complexes [LnIII{N(SiiPr3)2}2][B(C6F5)4] (3-Ln; Ln = Sm, Tm, Yb) were also treated with TEMPO• to yield the heteroleptic, approximately trigonal planar Ln3+ complexes [LnIII{N(SiiPr3)2}2(TEMPO•)][B(C6F5)4] (4-Ln; Ln = Sm, Tm, Yb); the cations of 4-Ln all contain TEMPO• radicals. We have compared the electronic structures of the two geometrically similar families of Ln3+ complexes with the TEMPO- anion (2-Ln) or TEMPO• radical (4-Ln) using a combination of UV-vis-NIR and EPR spectroscopy, magnetic measurements, and ab initio calculations. These studies revealed no single-molecule magnet behavior for 2-Yb despite evidence for sizable CF splitting and a high degree of purity of the ground stabilized mJ = |±7/2⟩ state, while the radical TEMPO• in 4-Yb did not significantly improve performance.

Generalized Path Integral Energy and Heat Capacity Estimators of Quantum Oscillators and Crystals Using Harmonic Mapping.

Imaginary-time path integral (PI) is a rigorous tool to treat nuclear quantum effects in static properties. However, with its high computational demand, it is crucial to devise precise estimators. We introduce generalized PI estimators for the energy and heat capacity that utilize coordinate mapping. While it can reduce to the standard thermodynamic and centroid virial (CVir) estimators, the formulation can also take advantage of harmonic character of quantum oscillators and crystals to construct a coordinate mapping. The method is not applicable to fluids or systems with imaginary modes such as double-well potentials. This yields harmonically mapped averaging (HMA) estimators, with mappings that decouple (HMAc) or couple (HMAq) the centroid and internal modes. The HMAq is constructed with normal mode coordinates (HMAq-NM) with quadratic scaling of cost or harmonic oscillator staging (HMAq-SG) coordinates with linear scaling. The estimator performance is examined for a 1D anharmonic oscillator and a 3D Lennard-Jones crystal using path integral molecular dynamics (PIMD) simulation. The HMA estimators consistently provide more precise estimates compared to CVir, with the best performance obtained by HMAq-NM, followed by HMAq-SG, and then HMAc. We also examine the effect of anharmonicity (for AO), intrinsic quantumness, and Trotter number. The HMA formulation introduced assumes the availability of forces and Hessian matrix; however, an equally efficient finite difference alternative is possible when these derivatives are inaccessible. The remarkable improvement in precision offered by HMAq estimators provides a framework for efficient PI simulation of more challenging systems, such as those based on ab initio calculations.

Electron scattering study on acetic acid and methyl formate

Abstract In this work we report the results of a theoretical calculation of the elastic, differential scattering, and excitation cross-sections on electron interactions with the C2H4O2 isomers (methyl formate and acetic acid) using the ab initio R-matrix method in the energy range of 0.1-20 eV. The computations were performed using static exchange (SE), static exchange plus polarization (SEP) and Close-Coupling (CC) models with electronic structure calculation for these molecules performed using GAMESS. In the electron scattering cross section we have identified π* type resonance in both the isomers. Ionization cross-sections for both the molecules from ionization threshold to 500 eV using BEB method are also presented here. We endeavoured to explore the isomeric effect on various cross sections among these two isomers and included the third isomer, namely glycolaldehyde, as reported in our previous publication.

Terahertz time-domain investigation of atacamite: Spectral analysis and theoretical insights for cultural heritage applications

Spectroscopic methods based on terahertz (THz) radiation can represent an innovative approach in the Cultural Heritage field. Specifically, contemporary and historical pigments can exhibit unique fingerprints. However, despite the high selectiveness of terahertz technique, the origin of spectral characteristics is not often known. In this work, the characteristic terahertz-time-domain spectroscopy spectrum of atacamite is reported in the spectral range 0.1–3.7 THz along with ab initio calculation of lattice dynamics to assess the nature of experimental features. The obtained results enable the acquisition of reliable new data in Heritage sciences, which is essential for developing accurate reference pigment databases.

Catching the π-Stacks: Prediction of Aggregate Structures of Porphyrin.

π-π interactions decisively shape the supramolecular structure and functionality of π-conjugated molecular semiconductor materials. Despite the customizable molecular building blocks, predicting their supramolecular structure remains a challenge. Traditionally, force field methods have been used due to the complexity of these structures, but advances in computational power have enabled ab initio approaches such as density functional theory (DFT). DFT is particularly suitable for finding energetically favorable structures of dye aggregates, which are determined by a large number of different interactions, but a systematic aggregate search can still be very challenging due to the large number of possible geometries. In this work, we show ways to overcome this challenge. We investigate how finely translational and rotational lattices must be structured to identify all energetic minima of π-stack structures, focusing on porphyrins as a prototype challenge. Our approach involves single-point DFT calculations of systematically varied dimer geometries, identification of local energy minima, hierarchical grouping of geometrically similar structures, and optimization of the energetically favorable representatives of each geometric family. This ab initio method provides a general framework for the systematic prediction of aggregate structures and reveals geometrically diverse and energetically favorable dimers.

Scalable bottom-up synthesis of Co-Ni-doped graphene.

Introducing heteroatoms into graphene is a powerful strategy to modulate its catalytic, electronic, and magnetic properties. At variance with the cases of nitrogen (N)- and boron (B)-doped graphene, a scalable method for incorporating transition metal atoms in the carbon (C) mesh is currently lacking, limiting the applicative interest of model system studies. This work presents a during-growth synthesis enabling the incorporation of cobalt (Co) alongside nickel (Ni) atoms in graphene on a Ni(111) substrate. Single atoms are covalently stabilized within graphene double vacancies, with a Co load ranging from 0.07 to 0.22% relative to C atoms, controllable by synthesis parameters. Structural characterization involves variable-temperature scanning tunneling microscopy and ab initio calculations. The Co- and Ni-codoped layer is transferred onto a transmission electron microscopy grid, confirming stability through scanning transmission electron microscopy and electron energy loss spectroscopy. This method holds promise for applications in spintronics, gas sensing, electrochemistry and catalysis, and potential extension to graphene incorporation of similar metals.