Search

Utilizing the sparsity of the electronic structure problem, fragmentation methods have been researched for decades with great success, pushing the limits of ab initio quantum chemistry ever further. Recently, this set of methods has been expanded to include a fundamentally different approach called excitonic renormalization, providing promising initial results. It builds a supersystem Hamiltonian in a second-quantized-like representation from transition-density tensors of isolated fragments, contracted with biorthogonalized molecular integrals. This makes the method fully modular in terms of the quantum chemical methods applied to each fragment and enables massive truncation of the state-space required. Proof-of-principle tests have previously shown that an excitonically renormalized Hamiltonian can efficiently scale to hundreds of fragments, but the ad hoc approach to building the Hamiltonian was not scalable to larger fragments. On the other hand, initial tests of the originally proposed modular Hamiltonian build, presented here, show the accuracy to be poor on account of its non-Hermitian character. In this study, we bridge the gap between these with an operator expansion that is shown to converge rapidly, tending toward a Hermitian Hamiltonian while retaining the modularity, yielding an accurate, scalable method. The accuracy is tested here for a beryllium dimer. At distances near equilibrium and longer, the zeroth-order method is comparable to coupled-cluster singles, doubles, and perturbative triples and the first-order method is comparable to full configuration interaction (FCI). The second-order method agrees with FCI for distances well up the inner repulsive wall of the potential. Deviations occurring at shorter bond distances are discussed along with approaches to scaling to larger fragments. [ABSTRACT FROM AUTHOR]

Since the long-range structure of the Gd-doped system (Ce1−xGdxO2−x/2) remains in the fluorite structure in the x-range [0–0.3], it is of interest to study this system in the local range. Molecular dynamics simulations seem to be particularly suitable for this purpose. Among all the possible parameters used to describe the local order/disorder, we focus on the behavior of the bond lengths. Doping cations and oxygen vacancies are introduced inside the fluorite structure following two arrangements, the random (R) and the neighbor (N). We investigate the evolution of the bond lengths with the temperature and doping level and identify different trends depending on the bond type (cation–anion, cation–cation). The bond length thermal expansion can be drastically different from the long-range thermal expansion mainly for the bonds implying Gd3+ cation. The presence of C-type Gd2O3 clusters is discussed through the study of the double component of the Gd–Gd Radial Distribution Function. Comments are provided for evolution of distances beyond the first and second shells of the cations' environment. This study also demonstrates the impact of point defect distribution (R or N) on the evolution of the local structure, while maintaining the long-range fluorite structure. [ABSTRACT FROM AUTHOR]

Photodissociation of the CH2Cl radical is investigated by using high-level multireference configuration interaction ab initio methods, including the spin–orbit coupling. All possible fragmentation pathways, namely, CH2Cl + hν → CH2 + Cl, HCCl + H, and CCl + H2, have been analyzed. The potential-energy curves of the ground and several excited electronic states along the corresponding dissociating bond distance of each pathway have been calculated. Inclusion of the spin–orbit couplings is found to be crucial because it strongly determines the shape of the curves of the different excited states and, therefore, their photodissociation dynamics behavior. Analysis of the potential curves indicates that the pathways producing CH2 + Cl and HCCl + H can occur through a fast direct dissociation mechanism, while the pathway leading to CCl + H2 involves much slower dissociation mechanisms such as internal conversion between electronic states, predissociation, or tunneling through exit barriers. The main implications are that the two faster channels are predicted to be dominant, while the slower pathway is expected to be very unlikely and rather irrelevant. Appreciable actinic fluxes of solar irradiation are available at stratospheric altitudes where ozone is abundant, in the wavelength range where absorption of the first low-lying excited states of CH2Cl has been observed experimentally. Our results show that in this excitation energy range, the above-mentioned two dominant dissociation pathways are open and then could contribute to stratospheric ozone depletion. [ABSTRACT FROM AUTHOR]

Glass and crystal are inherently different material states in terms of their structural and physical features; consequently, the direct quantitative connection between crystal and glass is lacking. Herein, we first show that the liquid fragility m, which is featured by the negative departure degree of viscosity with the temperature at the glass transition temperature (Tg), has a direct exponential correlation with the ratio of the average longest metal–oxygen and the average phosphorus, silicon, or boron–oxygen bond lengths of the crystal in various oxides including phosphates, silicates, and borates. Such a result can be rationalized by the fact that the fragility m in these glass-formers is associated with the total network rigidity determined by the weakest bond due to the "bucket effect" and the bond pair inheritance of glass from that of the crystal. Our work connects direct features between glass and crystal with identical composition, providing a new viewpoint bridging glass and crystal. [ABSTRACT FROM AUTHOR]

A linear chain of N particles connected by Morse bonds with periodic boundaries in a Brownian bath is considered as a model for polymer fracture. The potential of mean force (pomf) with respect to the length of a bond and its energy and entropic components are computed via Langevin Dynamics simulations at various chain elongations λ > 1. A narrow range of λ values is identified over which equilibrium is established between an intact (i) and a fractured (f) state over a pomf barrier. While the barrier and (f) regions of the pomf are well described by a well-known (N − 1)-bond adiabatic approximation, the shape of the (i) region departs from it, exhibiting a bimodal character. A new, 1-bond adiabatic approximation is proposed to explain this. The lower envelope of the two adiabatic approximations provides an excellent description of the pomf. The relationship between (f) states for individual bonds and for the entire chain is elucidated. A new approach, based on analysis of equilibrium tension autocorrelation functions, is developed to extract rate constants for the fracture and reformation of the chain in dependence of λ. Results from it agree with tension relaxation during nonequilibrium simulations initiated at equispaced configurations of the particles. Rate constants for chain fracture conform to a Boltzmann–Arrhenius–Zhurkov dependence on the tension averaged over the intact state of the chain, the activation length being higher than estimated from pomf extrema. These findings constitute a step toward a predictive multiscale simulation scheme for fracture in polymeric materials. [ABSTRACT FROM AUTHOR]

The coupled structural and electronic parameters of intramolecular hydrogen bonding play an important role in ultrafast chemical reactions, such as proton transfer processes. We perform one- and two-dimensional vibrational–electronic (1D and 2D VE) spectroscopy experiments to understand the couplings between vibrational and electronic coordinates in 10-Hydroxybenzo[h]quinoline, an ultrafast proton transfer system. The experiments reveal that the OH stretch (νOH) is strongly coupled to the electronic excitation, and Fourier analysis of the 1D data shows coherent oscillations from the low frequency backbone vibrational modes coupled to the νOH mode, resulting in an electronically detected vibronic signal. In-plane low-frequency vibrations at 242 and 386 cm−1 change the hydrogen bond distance and modulate the observed electronic signal in the polarization-selective 1D VE experiment through orientation-dependent coupling with the νOH mode. Resolution of the excitation frequency axis with 2D VE experiments reveals that excitation frequency, detection frequency, and experimental delay affect the frequency and strength of the vibronic transitions observed. Our results demonstrate evidence of direct coupling of the high frequency νOH mode with the S1 ← S0 electronic transition in 10-Hydroxybenzo[h]quinoline (HBQ), and orientation-dependent couplings of the low-frequency 242 and 386 cm−1 modes to the νOH mode and the electronic transition. This demonstration of multidimensional VE spectroscopy on HBQ reveals the potential of using 1D and 2D VE spectroscopy to develop a quantitative understanding of the role of vibronic coupling in hydrogen bonding and ultrafast proton transfer for complex systems. [ABSTRACT FROM AUTHOR]

Carbonyl groups (C=O) play crucial roles in the photophysics and photochemistry of biological systems. O1s x-ray photoelectron spectroscopy allows for targeted investigation of the C=O group, and the coupling between C=O vibration and O1s ionization is reflected in the fine structures. To elucidate its characteristic vibronic features, systematic Franck–Condon simulations were conducted for six common biomolecules, including three purines (xanthine, caffeine, and hypoxanthine) and three pyrimidines (thymine, 5F-uracil, and uracil). The complexity of simulation for these biomolecules lies in accounting for temperature effects and potential tautomeric variations. We combined the time-dependent and time-independent methods to efficiently account for the temperature effects and to provide explicit assignments, respectively. For hypoxanthine, the tautomeric effect was considered by incorporating the Boltzmann population ratios of two tautomers. The simulations demonstrated good agreement with experimental spectra, enabling differentiation of two types of carbonyl oxygens with subtle local structural differences, positioned between two nitrogens (O1) or between one carbon and one nitrogen (O2). The analysis provided insights into the coupling between C=O vibration and O1s ionization, consistently showing an elongation of the C=O bond length (by 0.08–0.09 Å) upon O1s ionization. [ABSTRACT FROM AUTHOR]

Based upon global-minimum searches and first-principles electronic structural calculations, we present the perfectly planar B6O50/−/2− (1–3) systems. Notably, the C2v B6O50/− (1–2) clusters mark a groundbreaking advancement as the first boron oxide clusters to feature a five-membered ring. Both structures exhibit remarkable similarity, each characterized by a central B3O2 ring surrounded by three terminal BO groups. However, a notable difference lies in the bond lengths within the B3O2 core, which are uniform in the former but uneven in the latter. Detailed canonical molecular orbital and AdNDP analyses reveal the delocalized π bonds over the B3O2 ring, affirming the π aromaticity of B6O5−. The Cs B6O52− (3) cluster showcases a chain-like arrangement with three terminal BO groups and a typical OBO group linked to a B2 core, emphasizing the crucial role of electrons in shaping the structure. Detailed chemical bonding analyses reveal the presence of dual three-center four-electron (3c–4e) π hyperbonds in 3, shedding further light on its unique bonding characteristics. Moreover, this study comprehensively extends the B(BO)n0/− (n = 1–5) series and provides valuable perspectives on their structural evolution. [ABSTRACT FROM AUTHOR]

We present an efficient first-principles based method geared toward reliably predicting the structures of solid materials across the Periodic Table. To this end, we use a density functional theory baseline with a compact, near-minimal min+s basis set, yielding low computational costs and memory demands. Since the use of such a small basis set leads to systematic errors in chemical bond lengths, we develop a linear pairwise correction, available for elements Z = 1–86 (excluding the lanthanide series), parameterized for use with the Perdew–Burke–Ernzerhof exchange–correlation functional. We demonstrate the reliability of this corrected approach for equilibrium volumes across the Periodic Table and the transferability to differently coordinated environments and multi-elemental crystals. We examine relative energies, forces, and stresses in geometry optimizations and molecular dynamics simulations. [ABSTRACT FROM AUTHOR]

A new computational framework for spinor-based relativistic exact two-component (X2C) calculations is developed using contracted basis sets with a spin–orbit contraction scheme. Generally contracted, j-adapted basis sets of p-block elements using primitive functions in the correlation-consistent basis sets are constructed for the X2C Hamiltonian with atomic mean-field spin–orbit integrals (the X2CAMF scheme). The contraction coefficients are taken from atomic X2CAMF Hartree–Fock spinors, thereby following the simple concept of a linear combination of atomic orbitals. Benchmark calculations of spin–orbit splittings, equilibrium bond lengths, and harmonic vibrational frequencies demonstrate the accuracy and efficacy of the j-adapted spin–orbit contraction scheme. [ABSTRACT FROM AUTHOR]

Photoexcitation is a powerful way to induce phase transition of strongly correlated materials and dynamically control terahertz (THz) devices integrated with photoinduced phase transition (PIPT) materials. To clarify controversies over the physical mechanism between electronic insulator-metal transition (IMT) and structural phase transition (SPT) of photoexcited vanadium dioxide (VO 2), the underlying atomic and electronic state changes during photoinduced monoclinic-to-rutile phase transition are illustrated, and the separation with different thresholds between the quasi-instantaneous IMT and the ultrafast SPT is discovered. Below the SPT threshold, there exist metastable states exhibiting the metal-like monoclinic phases, i.e., the strongest metallicity and weak monoclinic phases, when the bond lengths of the V–V pairs are closest. By analyzing the electronic transport properties of these metal-like monoclinic phases, the THz response of the whole phase transition process can be characterized for first time through the quantum-electromagnetic dispersion modeling method. The THz properties of the practical VO 2 film are simulated and the great alignments between the measurements and the simulations verify the proposed analysis method, which provides a powerful exploration path and insights for the theoretical analysis and design verification of PIPT materials and their optoelectronic THz devices. [ABSTRACT FROM AUTHOR]

Photoexcitation is a powerful way to induce phase transition of strongly correlated materials and dynamically control terahertz (THz) devices integrated with photoinduced phase transition (PIPT) materials. To clarify controversies over the physical mechanism between electronic insulator-metal transition (IMT) and structural phase transition (SPT) of photoexcited vanadium dioxide (VO 2), the underlying atomic and electronic state changes during photoinduced monoclinic-to-rutile phase transition are illustrated, and the separation with different thresholds between the quasi-instantaneous IMT and the ultrafast SPT is discovered. Below the SPT threshold, there exist metastable states exhibiting the metal-like monoclinic phases, i.e., the strongest metallicity and weak monoclinic phases, when the bond lengths of the V–V pairs are closest. By analyzing the electronic transport properties of these metal-like monoclinic phases, the THz response of the whole phase transition process can be characterized for first time through the quantum-electromagnetic dispersion modeling method. The THz properties of the practical VO 2 film are simulated and the great alignments between the measurements and the simulations verify the proposed analysis method, which provides a powerful exploration path and insights for the theoretical analysis and design verification of PIPT materials and their optoelectronic THz devices. [ABSTRACT FROM AUTHOR]

The spectrum of localized hole states (traps) in amorphous silicon nitride, a-Si3N4, is experimentally studied using the method of thermally stimulated depolarization. The experiment is compared with theoretical calculations using three models of the energy spectrum of traps: discrete spectrum (monoenergetic trap), continuous spectrum, and Gaussian trap distribution. The experiment is quantitatively described by a model of a discrete spectrum of traps with an energy of 1.15 eV and a width of no more than 0.01 eV. In the case of a continuous and Gaussian spectrum of traps, the contribution to depolarization is made by the deepest traps. The blurring of the trap energy level in a-Si3N4 due to the absence of long-range order (fluctuations in the Si–N bond length and fluctuations in the N–Si–N and Si–N–Si angles) does not exceed 0.01 eV. [ABSTRACT FROM AUTHOR]

Acetaldehyde is a very relevant atmospheric species whose photodissociation has been extensively studied in the first absorption band both experimentally and theoretically. Very few works have been reported on acetaldehyde photodissociation at higher excitation energies. In this work, the photodissociation dynamics of acetaldehyde is investigated by means of high-level multireference configuration interaction ab initio calculations. Five different fragmentation pathways of acetaldehyde are explored by calculating the potential-energy curves of the ground and several excited electronic states along the corresponding dissociating bond distances. The excitation energy range covered in the study is up to 10 eV, nearly the ionization energy of acetaldehyde. We intend to rationalize the available experimental results and, in particular, to elucidate why some of the studied fragmentation pathways are experimentally observed in the different excitation energy regions and some others are not. Based on the shape of the calculated potential curves, we are able to explain the main findings of the available experiments, also suggesting possible dynamical dissociation mechanisms in the different energy regions. Thus, the reported potential curves are envisioned as a useful tool to interpret the currently available experiments as well as future ones on acetaldehyde photodissociation at excitation wavelengths in the range studied here. [ABSTRACT FROM AUTHOR]

This article investigates the relaxation dynamics of the ionized 2a1 state of a water molecule within a water dimer. The study was motivated by findings from two previous pieces of research that focused on the relaxation behaviors of the inner-valence ionized water dimer. The present study discloses an observation indicating that water dimers display specific fragmentation patterns following inner-valence ionization, depending on the position of the vacancy. Vacancies were created in the 2a1 state of the proton-donating water molecule (PDWM) and proton-accepting water molecule (PAWM). Utilizing Born–Oppenheimer molecular dynamics simulations, the propagation of the 2a1 ionized state was carried out for both scenarios. The results revealed proton transfer occurred when the vacancy resided in the PDWM, accompanied by the closing of decay channels for O–H bond distance (RO–H) > 1.187 Å (matching Richter et al.'s findings). Conversely, when vacancy was on PAWM, we observed no closing of decay channels (aligning with Jahnke et al.'s findings). This difference translates to distinct fragmentation pathways. In PDWM cases, 2a1 state ionization leads to H3O+ −OH• formation. In contrast, PAWM vacancies result in decay pathways leading to H2O+–H2O+ products. [ABSTRACT FROM AUTHOR]

We optimize the internuclear geometry and electronic structure of a model chiral system to achieve a maximal photoelectron circular dichroism (PECD) in its one-photon ionization by circularly polarized light. The electronic structure calculations are performed by the single center method, while the optimization is done using quantum alchemy employing a Taylor series expansion. Thereby, the effect of bond lengths and uncompensated charge distributions on the chiral response of the model is investigated theoretically in some detail. It is demonstrated that manipulating a chiral asymmetry of the ionic potential may enhance the dichroic parameter (i.e., the PECD) of the randomly oriented model system well beyond β1 = 25%. Furthermore, we demonstrate that quantum alchemy is applicable to PECD despite the unusually strong coupling of spatial and electronic degrees of freedom and discuss the relative impact of the individual degrees of freedom in this model system. We define the necessary conditions for the computational design of PECD for real (non-model) chiral molecules using our approach. [ABSTRACT FROM AUTHOR]

We propose a novel implementation of the extended-dynamics equations for isothermal–isobaric ensemble in molecular dynamics, as the Martyna–Tobias–Klein thermostat and barostat. This method is suitable for systems with constraints and the Verlet-family integrators. Instead of iterations or the Trotter-expansion-based methods, both velocities and box sizes (scaling of bond lengths) are predicted. The algorithm begins with force calculation, requiring neither quarter nor half time steps, and necessitating iterations only inside SHAKE. Several tests demonstrate that the quality is comparable to other implementations. It is found that the formula relating the extended barostat mass to the characteristic time of volume fluctuations is inaccurate for condensed systems, which has consequences for the parameter setup. Emphasis is also put on the verification of the precise isothermal–isobaric ensemble and finite-size effects. [ABSTRACT FROM AUTHOR]

The conformational space of cis-1,2-cyclohexanediol, a model molecule for cyclic vicinal diols, was investigated using rotational spectroscopy and density functional theory calculations. Four low energy conformers within an energy window of 5 kJ mol−1 were identified computationally. A rotational spectrum of jet-cooled cis-1,2-cyclohexanediol was recorded with a chirped pulse Fourier transform microwave spectrometer. Two sets of rotational transitions were observed and could be assigned to conformers of cis-1,2-cyclohexanediol. The non-observation of other low energy conformers was explained by conformational conversion barrier height calculations and results from experimental spectra recorded with different carrier gases. Eight isotopologues, including those with 13C and 18O, of the lowest energy conformer were observed, allowing the determination of the semi-experimental equilibrium structure, reSE. Interestingly, the structural analysis revealed that the C–O bond length of the intramolecular hydrogen-bond donor is shorter than that of the acceptor. This appears to be a general characteristic of vicinal diols and can be used as a novel hydrogen-bond marker in such compounds. [ABSTRACT FROM AUTHOR]

Density functional approximations to the exchange–correlation energy can often identify strongly correlated systems and estimate their energetics through energy-minimizing symmetry-breaking. In particular, the binding energy curve of the strongly correlated chromium dimer is described qualitatively by the local spin density approximation (LSDA) and almost quantitatively by the Perdew–Burke–Ernzerhof generalized gradient approximation (PBE-GGA), where the symmetry breaking is antiferromagnetic for both. Here, we show that a full Perdew–Zunger self-interaction-correction (SIC) to LSDA seems to go too far by creating an unphysical symmetry-broken state, with effectively zero magnetic moment but non-zero spin density on each atom, which lies ∼4 eV below the antiferromagnetic solution. A similar symmetry-breaking, observed in the atom, better corresponds to the 3d↑↑4s↑3d↓↓4s↓ configuration than to the standard 3d↑↑↑↑↑4s↑. For this new solution, the total energy of the dimer at its observed bond length is higher than that of the separated atoms. These results can be regarded as qualitative evidence that the SIC needs to be scaled down in many-electron regions. [ABSTRACT FROM AUTHOR]

The structure of grain boundaries (GBs) of metals is essential in determining the solute segregation at GBs; however, its complexity prohibits the understanding of the underlying mechanism. We propose a geometric descriptor of GB segregation based on the non-local coordination number of cut surfaces from GBs, which determines the segregation energies of solutes at the grain boundaries of metals across multidimensional GB space, different solutes, and different matrices. The effectiveness of the descriptor originates from the correlation between bonding strength, d-band width, and coordination number. This descriptor only depends on the bond length and angle of pre-segregation sites at GBs and can be readily used for description and prediction. Our scheme builds a novel picture for understanding the role of GB structures in segregation and provides a useful tool for the design of advanced alloys. [ABSTRACT FROM AUTHOR]

VO2 is renowned for its electric transition from an insulating monoclinic (M1) phase, characterized by V–V dimerized structures, to a metallic rutile (R) phase above 340 K. This transition is accompanied by a magnetic change: the M1 phase exhibits a non-magnetic spin-singlet state, while the R phase exhibits a state with local magnetic moments. Simultaneous simulation of the structural, electric, and magnetic properties of this compound is of fundamental importance, but the M1 phase alone has posed a significant challenge to the density functional theory (DFT). In this study, we show none of the commonly used DFT functionals, including those combined with on-site Hubbard U to treat 3d electrons better, can accurately predict the V–V dimer length. The spin-restricted method tends to overestimate the strength of the V–V bonds, resulting in a small V–V bond length. Conversely, the spin-symmetry-breaking method exhibits the opposite trends. Each of these two bond-calculation methods underscores one of the two contentious mechanisms, i.e., Peierls lattice distortion or Mott localization due to electron–electron repulsion, involved in the metal–insulator transition in VO2. To elucidate the challenges encountered in DFT, we also employ an effective Hamiltonian that integrates one-dimensional magnetic sites, thereby revealing the inherent difficulties linked with the DFT computations. [ABSTRACT FROM AUTHOR]

VO2 is renowned for its electric transition from an insulating monoclinic (M1) phase, characterized by V–V dimerized structures, to a metallic rutile (R) phase above 340 K. This transition is accompanied by a magnetic change: the M1 phase exhibits a non-magnetic spin-singlet state, while the R phase exhibits a state with local magnetic moments. Simultaneous simulation of the structural, electric, and magnetic properties of this compound is of fundamental importance, but the M1 phase alone has posed a significant challenge to the density functional theory (DFT). In this study, we show none of the commonly used DFT functionals, including those combined with on-site Hubbard U to treat 3d electrons better, can accurately predict the V–V dimer length. The spin-restricted method tends to overestimate the strength of the V–V bonds, resulting in a small V–V bond length. Conversely, the spin-symmetry-breaking method exhibits the opposite trends. Each of these two bond-calculation methods underscores one of the two contentious mechanisms, i.e., Peierls lattice distortion or Mott localization due to electron–electron repulsion, involved in the metal–insulator transition in VO2. To elucidate the challenges encountered in DFT, we also employ an effective Hamiltonian that integrates one-dimensional magnetic sites, thereby revealing the inherent difficulties linked with the DFT computations. [ABSTRACT FROM AUTHOR]

Molecular-level nucleation has not been clearly understood due to the complexity of multi-body potentials and the stochastic, rare nature of the process. This work utilizes molecular dynamics (MD) simulations, incorporating a first-principles-based deep neural network (DNN) potential model, to investigate homogeneous water vapor condensation. The nucleation rates and critical nucleus sizes predicted by the DNN model are compared against commonly used semi-empirical models, namely extended simple point charge (SPC/E), TIP4P, and OPC, in addition to classical nucleation theory (CNT). The nucleation rates from the DNN model are comparable with those from the OPC model yet surpass the rates from the SPC/E and TIP4P models, a discrepancy that could mainly arise from the overestimated bulk free energy by SPC/E and TIP4P. The surface free energy predicted by CNT is lower than that in MD simulations, while its bulk free energy is higher than that in MD simulations, irrespective of the potential model used. Further analysis of cluster properties with the DNN model unveils pronounced variations of O–H bond length and H–O–H bond angle, along with averaged bond lengths and angles that are enlarged during embryonic cluster formation. Properties such as cluster surface free energy and liquid-to-vapor density transition profiles exhibit significant deviations from CNT assumptions. [ABSTRACT FROM AUTHOR]

The theory of analytic gradients is presented for the projector-based density functional theory (DFT) embedding approach utilizing the Huzinaga-equation. The advantages of the Huzinaga-equation-based formulation are demonstrated. In particular, it is shown that the projector employed does not appear in the Lagrangian, and the potential risk of numerical problems is avoided at the evaluation of the gradients. The efficient implementation of the analytic gradient theory is presented for approaches where hybrid DFT, second-order Møller–Plesset perturbation theory, or double hybrid DFT are embedded in lower-level DFT environments. To demonstrate the applicability of the method and to gain insight into its accuracy, it is applied to equilibrium geometry optimizations, transition state searches, and potential energy surface scans. Our results show that bond lengths and angles converge rapidly with the size of the embedded system. While providing structural parameters close to high-level quality for the embedded atoms, the embedding approach has the potential to relax the coordinates of the environment as well. Our demonstrations on a 171-atom zeolite and a 570-atom protein system show that the Huzinaga-equation-based embedding can accelerate (double) hybrid gradient computations by an order of magnitude with sufficient active regions and enables affordable force evaluations or geometry optimizations for molecules of hundreds of atoms. [ABSTRACT FROM AUTHOR]

The energetic stability of positron–dianion systems [A−; e+; A−] is studied via many-body theory, where A− includes H−, F−, Cl−, and the molecular anions (CN)− and (NCO)−. Specifically, the energy of the system as a function of ionic separation is determined by solving the Dyson equation for the positron in the field of the two anions using a positron–anion self-energy as constructed in Hofierka et al. [Nature 606, 688 (2022)] that accounts for correlations, including polarization, screening, and virtual-positronium formation. Calculations are performed for a positron interacting with H 2 2 − , F 2 2 − , and Cl 2 2 − and are found to be in good agreement with previous theory. In particular, we confirm the presence of two minima in the potential energy of the [H−; e+; H−] system with respect to ionic separation: a positronically bonded [H−; e+; H−] local minimum at ionic separations r ∼ 3.4 Å and a global minimum at smaller ionic separations r ≲ 1.6 Å that gives overall instability of the system with respect to dissociation into a H2 molecule and a positronium negative ion, Ps−. The first predictions are made for positronic bonding in dianions consisting of molecular anionic fragments, specifically for (CN) 2 2 − and (NCO) 2 2 − . In all cases, we find that the molecules formed by the creation of a positronic bond are stable relative to dissociation into A− and e+A− (positron bound to a single anion), with bond energies on the order of 1 eV and bond lengths on the order of several ångstroms. [ABSTRACT FROM AUTHOR]

31P nuclear magnetic resonance (NMR) chemical shifts were shown to be very sensitive to the basis set used at the geometry optimization stage. Commonly used energy-optimized basis sets for a phosphorus atom containing only one polarization d-function were shown to be unable to provide correct equilibrium geometries for the calculations of phosphorus chemical shifts. The use of basis sets with at least two polarization d-functions on a phosphorus atom is strongly recommended. In this paper, an idea of creating the basis sets purposed for the geometry optimization that provide the least possible error coming from the geometry factor of accuracy in the resultant NMR shielding constants is proposed. The property-energy consisted algorithm with the target function in the form of the molecular energy gradient relative to P–P bond lengths was applied to create new geometry-oriented pecG-n (n = 1, 2) basis sets for a phosphorus atom. New basis sets have demonstrated by far superior performance as compared to the other commonly used energy-optimized basis sets in massive calculations of 31P NMR chemical shifts carried out at the gauge-including atomic orbital-coupled cluster singles and doubles/pecS-2 level of the theory by taking into account solvent, vibrational, and relativistic corrections. [ABSTRACT FROM AUTHOR]

The correlation consistent basis sets (cc-pVnZ with n = D, T, Q, 5) for the Ga–Br elements have been redesigned, tuning the sets for use for density functional approximations. Steps to redesign these basis sets for an improved correlation energy recovery and efficiency include truncation of higher angular momentum functions, recontraction of basis set coefficients, and reoptimization of basis set exponents. These redesigned basis sets are compared with conventional cc-pVnZ basis sets and other basis sets, which are, in principle, designed to achieve systematic improvement with respect to increasing basis set size. The convergence of atomic energies, bond lengths, bond dissociation energies, and enthalpies of formation to the Kohn–Sham limit is improved relative to other basis sets where convergence to the Kohn–Sham limit is typically not observed. [ABSTRACT FROM AUTHOR]

Experimental findings by others suggest that OH radicals are formed in unexpected abundance on or near surfaces of 1–50 µm microdroplets comprised of pure water, but the mechanism by which these radicals are generated is not yet fully resolved. In this work, we examine two possibilities using ab initio electronic structure methods: (1) electron transfer (ET) from a microdroplet surface-bound OH− anion to a nearby H3O+ cation and (2) proton transfer (PT) from such a H3O+ cation to a nearby OH− anion. Our findings suggest that both processes are possible but only if the droplet's underlying water molecules comprising the microdroplet provide little screening of the Coulomb interaction between the anion and cation once they reach ∼10 Å of one another. In the ET event, an OH radical is formed directly; for PT, the OH formation occurs because the new O–H bond formed by the transferred proton is created at a bond length sufficiently elongated to permit homolytic cleavage. Both the ET and PT pathways predict that H atoms will also be formed. Finally, we discuss the roles played by strong local electric fields in mechanisms that have previously been proposed and that occur in our two mechanisms. [ABSTRACT FROM AUTHOR]

Semi-experimental structures ( r e SE ) are derived from experimental ground state rotational constants combined with theoretical vibrational corrections. They permit a meaningful comparison with equilibrium structures based on high-level ab initio calculations. Typically, the vibrational corrections are evaluated with second-order vibrational perturbation theory (VPT2). The amount of error introduced by this approximation is generally thought to be small; however, it has not been thoroughly quantified. Herein, we assess the accuracy of theoretical vibrational corrections by extending the treatment to fourth order (VPT4) for a series of small linear molecules. Typical corrections to bond distances are on the order of 10−5 Å. Larger corrections, nearly 0.0002 Å, are obtained to the bond lengths of NCCN and CNCN. A borderline case is CCCO, which will likely require variational computations for a satisfactory answer. Treatment of vibrational effects beyond VPT2 will thus be important when one wishes to know bond distances confidently to four decimal places (10−4 Å). Certain molecules with shallow bending potentials, e.g., HOC+, are not amenable to a VPT2 description and are not improved by VPT4. [ABSTRACT FROM AUTHOR]

π-conjugated polymers have been used in a wide range of practical applications, partly due to their unique properties that originate in the delocalization of electrons through the polymer backbone. The level of delocalization can be characterized by the induced bond length alternation (BLA), with shorter BLA connected with strong delocalization and vice versa. The accurate description of this structural parameter can be considered a benchmark for testing the capability of different electronic structure methods for self-interaction error (SIE) removal and electron correlation inclusion. Density functional theory (DFT), in its local or semi-local flavors, suffers from SIE and, thus, underestimates the BLA compared to self-interaction-free methods. In this work, we utilize the Fermi–Löwdin orbital self-interaction correction (FLOSIC) method for one-electron self-interaction removal to characterize the BLA of five oligomers with increasing length extrapolated to the polymeric limit. We compare the self-interaction-free BLA to several DFT approximations, Møller–Plesset second-order perturbation theory (MP2), and the BLA obtained with the domain based local pair natural orbital CCSD(T) [DLPNO-CCSD(T)] approximation. Our findings show that FLOSIC corrects for the small BLA given by (semi-)local DFT approximations, but it tends to overcorrect with respect to CAM-B3LYP, MP2, and DLPNO-CCSD(T). [ABSTRACT FROM AUTHOR]

Pseudopotentials (PP) are extensively used in electronic structure calculations, particularly for molecules containing heavy elements. Parameters in PPs are mainly determined from ab initio results, and errors of such PPs in density functional theory (DFT) calculations have been studied previously. However, PP errors on results with spin–orbit coupling and those in time-dependent DFT (TDDFT) calculations have not been reported previously. In this work, we investigate the error of the small-core energy-consistent Stuttgart/Koln pseudopotentials in DFT and TDDFT calculations with and without spin–orbit coupling. Ground state bond lengths, harmonic frequencies, dissociation energies, and vertical excitation energies for a series of closed-shell diatomic heavy and superheavy p-block molecules are calculated using several popular exchange-correlation functionals. PP errors are estimated by comparing with results using the all-electron Dirac–Coulomb (-Gaunt) Hamiltonian. Our results show that the difference between ground state properties and most excitation energies in scalar-relativistic calculations with the PP and those of all-electron calculations is quite small. This difference becomes somewhat larger when spin–orbit coupling (SOC) is present, especially for properties that are affected by SOC to some extent. In addition, the errors of the PPs are insensitive to the employed exchange-correlation functionals in most cases. Our results indicate that reasonable DFT and TDDFT results can be obtained using the small-core energy-consistent Stuttgart/Koln pseudopotentials for heavy and super-heavy p-block molecules. [ABSTRACT FROM AUTHOR]

Motivated by the recent experimental study on hydrogen storage in MXene multilayers [Liu et al., Nat. Nanotechnol. 16, 331 (2021)], for the first time we propose a workflow to computationally screen 23 857 compounds of MXene to explore the general relation between the activated H2 bond length and adsorption distance. By using density functional theory we generate a dataset to investigate the adsorption geometries of hydrogen on MXenes, based on which we train physics-informed atomistic line graph neural networks (ALIGNNs) to predict adsorption parameters. To fit the results, we further derived a formula that quantitatively reproduces the dependence of H2 bond length on the adsorption distance from MXenes within the framework of Pauling's resonating valence bond theory, revealing the impact of transition metal's ligancy and valence on activating dihydrogen in H2 storage. [ABSTRACT FROM AUTHOR]

Supported Cu 3 cluster on N -doped graphene is a quite promising catalyst for CO electrochemical reduction to propanol. [Display omitted] Electroreduction of carbon monoxide into high-energy fuel is an excellent energy strategy for sustainable development, but the high yield of multi-carbon products remains a difficult challenge. Inspired by the successful synthesis of various trimer metal clusters and studies on electrocatalysis of CO to C3 products by Cu-based catalysts, Cu 3 supported on N -doped graphene structures (Cu 3 @NG) are investigated as electrocatalysts for CORR toward propanol in this paper. Due to the appropriate Cu-Cu bond length, the moderate charge of the Cu site, mild CO adsorption energy, and 100 % CO coverage, the absorbed 3*CO substance can form the critical *CO-CO-CO intermediate with a rather low kinetic barrier of 0.25 eV. The limiting potential of the whole process for the formation of propanol is just −0.15 V. Our work not only showed that Cu 3 @NG is an excellent catalyst for the formation of propanol with high selectivity at low potential but also indicated that the *CO coverage can greatly reduce the CO hydrogenation potential and bonding of some intermediates such as *CH 2 O. This research will provide a new idea for the material design of products tending to multi-carbon. [ABSTRACT FROM AUTHOR]

Mn4+/Tb3+ doped CaAl 2 Si 2 O 8 (CASO) phosphors were synthesized using high temperature solid-state method. Crystal structure and lattice parameters were acquired through X-ray diffraction technique and Rietveld refinement. Energy band and density of states of CASO host were given by density functional theory calculations. By calculating the bond lengths and total energies of Mn4+ doped CASO, the results show that Mn4+ occupies the Ca2+ site with six coordinating oxygen atoms in the phosphor. The CASO:Mn4+,Tb3+ phosphors exhibit two major emission bands originating from the red emitting of Mn4+ (2E g → 4A 2g) and the green emitting of Tb3+ (5D 4 → 7F 5). Due to the emitting intensity of Mn4+ and Tb3+ varies differently with temperature, the temperature sensing study was carried out adopting the fluorescence intensity ratio of Tb to Mn. The relative sensitivity of CASO:1.5%Mn4+,6%Tb3+ achieves a peak value of 3.38 % K−1 at 423 K. The repeatability was tested and the result showed that the sample exhibits excellent repeatability. The luminescence color of the employed phosphors transitions from orange to green as the temperature increases. These findings uncover potential applications of CASO:Mn4+,Tb3+ phosphors in temperature sensing and multicolor-tunable emitting. [ABSTRACT FROM AUTHOR]

Mg3Cr2Al18 (abbreviated in this report as MCA) is the parent phase for a large class of intermetallic compounds that belong to the cubic crystal space group, $Fd\bar{3}m$ Fd 3 ¯ m. The purpose of this paper is to introduce an ideal, unrelaxed crystal structure for compounds isostructural with MCA. There are five distinct atomic sublattices in MCA compounds, which can be denoted, A, B, C, D, and E. With this, a general description for MCA structures can be written as $A_1^{8a} B_2^{16c} C_2^{16d} D_6^{48f} E_{12}^{96g}$ A 1 8 a B 2 16 c C 2 16 d D 6 48 f E 12 96 g , where the superscripts represent the Wyckoff special equipoints associated with the various sublattices in MCA, and the subscripts indicate the contributions of each sublattice to the stoichiometry of one formula unit in an any given MCA-structured compound. Sublattices D and E are where deviations from ideality occur in real, MCA-like compounds. This paper examines MCA bond lengths, nearest-neighbour polyhedral arrangements, 3-D sublattice crystal structures, 2-D atom tessellation patterns, and crystal chemical effects associated with atomic relaxations on the D and E sublattices. The ideal MCA crystal structure developed in this report provides an appropriate initial structure for use as input to crystal structure refinements of diffraction data for MCA-like phases being examined experimentally, or as input for computational, atomistic simulations of the structures of such compounds. [ABSTRACT FROM AUTHOR]

Mg3Cr2Al18 (abbreviated in this report as MCA) is the parent phase for a large class of intermetallic compounds that belong to the cubic crystal space group, $Fd\bar{3}m$ Fd 3 ¯ m. The purpose of this paper is to introduce an ideal, unrelaxed crystal structure for compounds isostructural with MCA. There are five distinct atomic sublattices in MCA compounds, which can be denoted, A, B, C, D, and E. With this, a general description for MCA structures can be written as $A_1^{8a} B_2^{16c} C_2^{16d} D_6^{48f} E_{12}^{96g}$ A 1 8 a B 2 16 c C 2 16 d D 6 48 f E 12 96 g , where the superscripts represent the Wyckoff special equipoints associated with the various sublattices in MCA, and the subscripts indicate the contributions of each sublattice to the stoichiometry of one formula unit in an any given MCA-structured compound. Sublattices D and E are where deviations from ideality occur in real, MCA-like compounds. This paper examines MCA bond lengths, nearest-neighbour polyhedral arrangements, 3-D sublattice crystal structures, 2-D atom tessellation patterns, and crystal chemical effects associated with atomic relaxations on the D and E sublattices. The ideal MCA crystal structure developed in this report provides an appropriate initial structure for use as input to crystal structure refinements of diffraction data for MCA-like phases being examined experimentally, or as input for computational, atomistic simulations of the structures of such compounds. [ABSTRACT FROM AUTHOR]

This research presents, for the first time, a comprehensive and rigorous investigation of ruthenium(II) chalcogenonitrosyl bonding situations in two sets of coordination compounds: [Ru(NE)Cl2(LOEt)] (1a–4a) and [Ru(NE)Cl2(LOEt)]− (1b–4b), where E = O, S, Se, Te. Prior to and following the one-electron reduction, the Ru–NE bonding situations were subjected to analysis. The calculated geometric parameters indicate that both the Ru–NE and N–E bond lengths are susceptible to variation depending on the nature of the chalcogen employed. Furthermore, the results demonstrate that the one-electron reduction process serves to diminish the NE double bond character. The generalized Kohn–Sham energy decomposition analysis (GKS-EDA) was conducted to illustrate the Ru–NE bonding scenarios prior to and following the one-electron reduction. The results provide valuable insights into the nature of Ru(II)–NE (E = O, S, Se, Te) bonds, the influence of chalcogens on ruthenium compounds, as well as how the one-electron reduction affects the release of NE groups. The main findings indicate that the total interaction energy, prior to the one-electron reduction, is three to four times more stabilizing than in the reduced analogs, confirming that the reduction unequivocally enhances the lability of the Ru–NE bond even when heavier chalcogen analogues are employed. [ABSTRACT FROM AUTHOR]

The chemical bond lengths and angles of group 16 dihydrides were investigated. The relativistic effects are essential for heavy elements molecules calculations. Here we implement two relativistic effects, that is, scalar relativistic effects and spin‐orbit coupled zero‐order regular approximation. Concerning molecular symmetry, scalar relativistic effects and spin‐orbit relativistic effects show different descriptions. They are single group and double group, respectively. In addition, non‐relativistic effects were used for very weak relativistic effects on molecules and for comparing with and without relativistic effects for heavy element molecules. From H2O to PoH2, the bonding lengths and angles are due to sp hybridization orbitals, while LvH 2$$ {}_2 $$ is mainly due to p‐orbital bonding, resulting in a different configuration of bond lengths and angles than other group 16 dihydrides. The chemical bonding of group 16 dihydrides was analyzed from a single group point of view by operating the double group results to a single group. [ABSTRACT FROM AUTHOR]

The physical origins of composition-, temperature-, and size-motivated changes in refractive index in crystals have long been a puzzle. Combining the bond-order-length-strength theory, local bond average approach, and core–shell structural model, we investigated the refractive indexes in dependencies of composition, temperature, and size for the ternary wurtzite group-Ⅲ nitride alloys. The theoretical reproduction of the observations disclosed that (i) the doping of small atoms caused the contraction in bond length, the strengthening in bond energy, and the decrease of refractive index, whereas the doping of large atoms led to an elongation of bond length, a weakening of bond energy, and an increase of refractive index; (ii) the refractive index is inversely proportional to the cohesive energy and the cube of the Debye temperature; and (iii) with the gradual decrease in solid size, the coordination number lowers, the bond length contracts, the bond energy gains, the surface-to-volume ratio rises, and the refractive index decreases. The proposed formulation not only shows an in-depth comprehension of the physical essence of the stimuli impact on the refractive index but also is expected to be conducive to the exploitation, optimization, and operation of the new-type photonic, piezoelectric, and pyroelectric nanometer devices for the ternary wurtzite alloys. [ABSTRACT FROM AUTHOR]

Oxidative stress, associated with excessive production of reactive oxygen and nitrogen species (ROS, RNS), contributes to the development and progression of many ailments, such as aging, cardiovascular diseases, Alzheimer's disease, Parkinson's disease, diabetes, cancer, preeclampsia or multiple sclerosis. While phenols and polyphenols are the most studied antioxidants structurally similar compounds such as anilines or thiophenols are sporadically analyzed despite their radical scavenging potential. This work assesses the impact of structural features of phenols and thiophenols on their antioxidant activity. Seventeen pairs of phenol/thiophenol analogues, possessing both electron-donating and withdrawing groups were selected for this study. Several physicochemical properties of the compounds were determined by density functional theory (DFT) calculations at the (U)B3LYP/6-311++G(d,p) level of theory for gas phase calculations and at the (U)B3LYP/6-311++G(d,p) scrf = (smd, solvent = water) level for the solvated ones. Correlations between calculated properties and experimental radical scavenging activities were investigated to identify the pivotal physical characteristics contributing to antioxidant efficiency. These include S–H and O–H bond distances and bond dissociation enthalpies (BDE), dipole moments, logP values, highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) orbital energies, and the HOMO–LUMO gap energies that were calculated at the M06-2X/6-311++G(d,p) level of theory, and Fukui functions. The experimental activity was evaluated using the 2,2′-azino-bis(3-ethylbenzthiazoline-6-sulphonic acid) (ABTS) and 2,2-diphenyl-1picrylhydrazyl (DPPH) radical scavenging assays. Several compounds exhibited superior scavenging abilities, surpassing that of the reference antioxidant Trolox. The extensive DFT calculations revealed that in the gas phase, lower BDE values, compared to IP and PA, suggested that the HAT mechanism predominates in case of these compound groups. In contrast, in water, significant reductions in PA due to solvent effects suggested that the SPLET mechanism is dominant under aqueous conditions. [ABSTRACT FROM AUTHOR]

Lead‐free tin halide perovskites, particularly Cs2SnBr6, are gaining significant attention for their potential in optoelectronic applications. In this study, we investigate the material's electrical and vibrational properties under varying illumination conditions, providing novel insights into the impact of light on its conduction mechanisms. Cs2SnBr6 was synthesized via liquid‐phase and solid‐state reactions, and we report comprehensive analyses of its impedance, AC conductivity, and dielectric properties over a wide frequency range. Notably, our results reveal that light‐induced photogenerated charge carriers enhance conductivity, which is well explained by the overlapping large‐polaron tunneling (OLPT) model, while in darkness, the material follows the correlated barrier hopping (CBH) model. Moreover, Raman spectroscopy highlights structural changes, including a slight shortening of Sn–Br bond distances under illumination, which influences the vibrational frequencies and intensities. These findings emphasize the crucial role of illumination in tuning the electrical and dielectric responses of Cs2SnBr6, thus demonstrating its considerable promise for future optoelectronic applications, such as solar cells and photodetectors. [ABSTRACT FROM AUTHOR]

[Fe(L)2Cl(H2O)] (1), [Co(L)2(H2O)2] (2), [Ni(L)2] (3), and [Cu(L)2] (4) complexes of Schiff base ligand (L), (Z)‐2‐(1‐((furan‐2‐ylmethyl)imino)ethyl)phenol were synthesized and characterized using Spectro analytical techniques. Studies on DNA binding revealed an intercalation binding mode with a significant binding potential that is supported by intrinsic (Kb) and stern‐volmer (Ksq) quenching constants. Nuclease activity was assessed using gel‐electrophoresis with CT–DNA. In silico drug binding studies were carried out using the AutoDock Vina, and drug interactions were visualized using the BIOVIA Discovery Studio Visualizer. Pharmacokinetic properties were evaluated using the Swiss ADMET online server tool. The title compounds were modelled using the Gaussian 09 software, optimizing bond lengths and angles. The theoretical data from in silico docking studies are experimentally correlated with in vivo carrageenan‐induced anti‐inflammatory screening, which demonstrated a higher rate of inhibition in paw edema for complex 4. Studies on HeLA and A549 cell lines for anticancer effect indicated a promising activity for complex 4. Antimicrobial assays against bacterial and fungal strains showcased greater activity levels of the complexes in comparison with its ligand. [ABSTRACT FROM AUTHOR]

The hot corrosion behaviors of the Sc2O3–Gd2O3–Y2O3 co‐doped ZrO2 (SGYSZ) ceramics in Na2SO4 + V2O5 molten salts at high temperature were investigated in this work. The results show that the SGYSZ exhibits good hot corrosion resistance to the 50 wt% Na2SO4 + 50 wt% V2O5 molten salts at 1000°C. It is attributed to the doping of small Sc3+, which can decrease the M–O bonding length and increase the binding energy, and then reduce the dissolution and loss of the stabilizer as rare‐earth cations Gd3+ and Y3+ from the SGYSZ. The significant synergistic effects of Sc3+, Gd3+, and Y3+ are beneficial to enhance the phase stability and corrosion resistance of the zirconate ceramic due to the different ionic radius and basicity of rare‐earth oxides (RE2O3). More Gd3+ in SGYSZ can promote the forming and coarsening of the GdVO4 during mineralization, which can also decrease the loss of Y3+ from the SGYSZ. The corrosion mechanism is discussed in detail. [ABSTRACT FROM AUTHOR]

As the major hydration product of cement, hydrated calcium silicate (C‐S‐H) governs the overall performance of cement‐based materials. The molar ratio of CaO to SiO2 (Ca/Si ratio) significantly affects the structure and properties of C‐S‐H. This study analyzed the effect of Ca/Si ratios (0.83–2.0) on the structural morphology evolution, bond lengths and angles, polymerization process, and nanoporosity of amorphous C‐S‐H, with the help of the ReaxFF force field. The results showed that the reacted C‐S‐H tend to form a fibrous network‐like morphology at low Ca/Si ratios, while the silicate chains are prone to accumulating at high Ca/Si ratios, forming a dense granular ovoid structure. Meanwhile, the Ca/Si ratio has no effect on the bond lengths and angles. In addition, the Ca2+ ions can interrupt the silicate chains during hydration, which leads to a decrease in the average silicate chain length with increasing Ca/Si ratio. The porosity of C‐S‐H decreases from 59.3% to 54.3% when the Ca/Si ratio increases from 0.83 to 2.0. It can be deduced from these findings that the increase in the Ca/Si ratio decreases the compressive strength of cement‐based materials but increases their durability. [ABSTRACT FROM AUTHOR]

Topological constraint theory (TCT) and molecular dynamics simulations (MD) are utilized to study the effect of Y2O3/SiO2 substitution on the short‐ and medium‐range structures and properties in aluminosilicate glasses containing Y2O3. Experimental methods have been applied to characterize the structures of the as‐prepared glasses to verify the accuracy of MD simulations. It was found that the Y–O bond distance is around 2.38 Å and the Y average coordination number is around 5.78. The introduction of Y2O3 disrupts the Si–O and Al–O bonds in the network structure and elevates the proportion of non‐bridging oxygen. Meanwhile, the calculation results of Si–O and Al–O bond distances show Y shortens the distance between Si, Al cations and oxygen ions, improving the network structure tightness. Furthermore, the glass transition temperature (Tg), Vickers hardness (HV), and elastic modulus (E) of glass increase with the increase of Y2O3 content, while the viscosity decreases. TCT was used to analyze the relationship among glass compositions, structures, and properties, and the properties prediction models for Tg, HV, and E were established. The effectiveness of TCT prediction models was proved by the comparison between the results predicted by the models and the data reported in the literature. [ABSTRACT FROM AUTHOR]

Electrocatalytic hydrogen generation requires a cutting‐edge, multifunctional electrocatalyst with abundant active sites to drive multielectron transfer reactions. High entropy alloys (HEA) are complex five or more‐element alloy systems with high configurational entropy which makes them unique materials for next‐generation electrocatalysts. Here, in this work, based on new screening guidelines for catalyst selections that combine density‐functional theory calculated Gibbs formation‐enthalpy with bond length and electronegativity variance, a novel HEA electrocatalyst consisting of five elements, namely, Pt, Pd, Ni, Co, and Mn is designed. By simple room temperature electrodeposition, the designed catalyst is prepared and its hydrogen evolution reaction (HER) is explored and validated through experimental and theoretical approaches. The HEA demonstrated a superior HER activity with an overpotential of 22.6 mV at −10 mA cm−2 which outperforms Pt/C commercial catalyst. No evident degradation of the material is detected even after 100 hours of continuous operation under high current density. Moreover, the HEA has shown exceptional performance in harsh electrolyte conditions such as in simulated seawater and actual seawater. Remarkably, the density‐functional theory calculated small Gibbs formation‐enthalpy (≈0 eV) compared to Pt/C places the new HEA near the apex of Trasatti's model of Volcano plot, which is also suggestive of superior HER activity. [ABSTRACT FROM AUTHOR]

The search for stable compounds containing an antiaromatic cyclic 4π system is a challenge for inventive chemists that can look back on a long history. Here we report the isolation and characterization of the novel 4π-electron tetrasilacyclobutadiene, an analogue of a 4π neutral cyclobutadiene that exhibits surprising features of a Möbius-type aromatic ring. Reduction of RSiCl3 (R = (iPr)2PC6H4) with KC8 in the presence of cycloalkyl amino-carbene (cAAC) led to the formation of corresponding silylene 1. Compound 1 shows further reactivity when treated with RSiCl3 under reducing conditions resulting in the formation of unsymmetrical bis-silylene 2, which was isolated as a dark red crystalline solid. Compound 3 was obtained when chlorosilylene was reduced with potassium graphite in the presence of 2. Although 3 is, according to Hückel's rule an antiaromatic system it possesses significant aromatic character due to the unusual delocalization of the HOMO-1 and the loss of degeneracy of the higher-lying π MOs. The aromatic character of 3 is indicated by the structural stability of the compound by the very similar Si-Si bond lengths and by the NICS values. There is an unusual π conjugation between the 4 π electrons in the nearly square-planar Si4 ring where the delocalization of the HOMO-1 occurs at two opposite sides of the ring. The reversal of the π conjugation resembles the twisting of the π conjugation in Möbius aromatic systems and it contributes to the stability of the compound. Introducing main group elements into neutral and charged cyclobutadienes enables the fine-tuning of their stability and electronic structures by disrupting their aromatic and antiaromatic character. Here, the authors report the isolation and characterization of a 4π-electron tetrasilacyclobutadiene, an analogue of a 4π neutral cyclobutadiene that exhibits features of a Möbius-type aromatic ring. [ABSTRACT FROM AUTHOR]

Glass substrate for organic light‐emitting diodes (OLED) requires to undergo several heating processes and thus makes a thermal shrinkage of glass substrate. The shrinkage leads to a pixel shifting and becomes a key factor affecting OLED panel quality. To understand the thermal shrinkage process, it is crucial to study the structural relaxation and thermodynamic characteristics of glass substrate around the glass transition temperature. In this work, the network structure of alkali‐free aluminosilicate glasses and their structural relaxation in the glass transition temperature range were studied by both experiment and molecular dynamics (MD) simulations. The alkali‐free aluminosilicate glasses with 15.3–22.1 mol% alkali‐earth oxides (RO) were prepared, and their reheating shrinkage and enthalpy relaxation were measured to reflect the structural relaxation during thermal cycling. As the RO content increased from 15.3 to 22.1 mol%, the CTE25–300°C increased from 3.00 × 10−6 to 4.13 × 10−6/°C, and the strain point decreased from 780 to 744°C. The reheating shrinkage of 4 × 4 × 40 mm glass rods also increased from 4.01 to 9.51 ppm, and the volume shrinkage primarily occurred in the initial stages of holding at 600°C. The enthalpy relaxation of the glass was measured using differential scanning calorimetry, and there were very small enthalpy relaxation in alkali‐free aluminosilicate glasses due to their good thermal stability. Therefore, MD simulations were employed to calculate the potential energy and evaluate the structural relaxation. The potential energy of the glasses during thermal cycling revealed a significant energy hysteresis upon reheating, and the structural relaxation occurred at relatively low temperature with the increasing addition of RO contents. Further analysis confirmed that the structural relaxation of glasses was mainly due to the structural arrangement of alkali‐earth metal ions, including the changes in bond length, angles, coordination number, and atomic mean‐square distance. [ABSTRACT FROM AUTHOR]

The manifestation of mechanical phenomena in quantum materials at the macroscopic level is intricately linked to pronounced electron‐electron interactions within their lattices, a relationship that becomes especially evident in low‐dimensional materials. Violet phosphorous (VP), a nascent 2D material distinguished by its unique vertically aligned tubular structures, has garnered considerable attention owing to its layer‐dependent electronic bandgap, exceptional carrier mobility, and robust air stability. Herein, a comprehensive exploration of the phonon modes exhibited by few‐layer VP through an integrated experimental‐theoretical approach, focusing on the modulation of their Raman response under the uniaxial strain along a‐axis, b‐axis, and tube direction, respectively, is undertaken. Density functional theory calculations highlight when strain is applied along the a‐ or b‐axis direction, the strain is predominantly mitigated through tube rotational adaptations instead of the change of bond length and bond angle, culminating in a pronounced anisotropic Raman response. Moreover, the strain engineering can effectively optimize the photoelectric response performance of VP, including increase the responsivity ≈2500% and elevate the anisotropic ratio from 2.26 to 3.38. This investigation not only confirms the superior stretchability and impact resistance properties of cross‐structured VP but also establishes the groundwork for exploring the strain‐induced anisotropic optoelectric properties to VP. [ABSTRACT FROM AUTHOR]

Patterns created by buckling of thin films on substrates have important applications in many fields such as functional devices and surface engineering. Introducing hierarchy to such patterns may further enhance their functionalities. However, unpredictable geometrical and material defects may disturb the formation of patterns. It is thus attractive to construct robust hierarchical buckling patterns that are insensitive to defects. Here, inspired by the construction of topologically protected localized states insensitive to defects in wave mechanics, a thin film of finite length bonded on a simple Winkler substrate is quantitatively designed by just varying its width continuously. Under the action of specific compressive stresses within such a continuous thin film owing to thermal mismatch between the film and substrate, there really exists a hierarchical film-buckling mode generated by topologically protected localization. Conditions under which the topologically protected localization exists are given explicitly, and the underlying mechanism is intuitively demonstrated from the mechanical point of view. Moreover, an analytical solution for the decay rate of localization is obtained as a satisfactory second-order approximation, from which the associated influencing factors are clearly identified. This work is expected to provide a theoretical guidance for the design and optimization of hierarchical patterns in film/substrate systems via the concept of topology. [ABSTRACT FROM AUTHOR]