Real-Space Interpretation of Interatomic Charge Transfer and Electron Exchange Effects by Combining Static and Kinetic Potentials and Associated Vector Fields**

Intricate behaviour of one-electron potentials from the Euler equation for electron density and corresponding gradient force fields in crystals was studied. Channels of locally enhanced kinetic potential and corresponding saddle Lagrange points were found between chemically bonded atoms. Superposition of electrostatic and kinetic potentials and electron density allowed partitioning any molecules and crystals into atomic - and potential-based -basins; -basins explicitly account for the electron exchange effect, which is missed for -ones. Phenomena of interatomic charge transfer and related electron exchange were explained in terms of space gaps between zero-flux surfaces of - and -basins. The gap between - and -basins represents the charge transfer, while the gap between - and -basins is a real-space manifestation of sharing the transferred electrons caused by the static exchange and kinetic effects as a response against the electron transfer. The regularity describing relative positions of -, -, and - basin boundaries between interacting atoms was proposed. The position of -boundary between - and -ones within an electron occupier atom determines the extent of transferred electron sharing. The stronger an H⋅⋅⋅O hydrogen bond is, the deeper hydrogen atom's -basin penetrates oxygen atom's -basin, while for covalent bonds a -boundary closely approaches a -one indicating almost complete sharing of the transferred electrons. In the case of ionic bonds, the same region corresponds to electron pairing within the -basin of an electron occupier atom.     

On the accurate estimation of free energies using the jarzynski equality

The Jarzynski equality is one of the most widely celebrated and scrutinized nonequilibrium work theorems, relating free energy to the external work performed in nonequilibrium transitions. In practice, the required ensemble average of the Boltzmann weights of infinite nonequilibrium transitions is estimated as a finite sample average, resulting in the so-called Jarzynski estimator, . Alternatively, the second-order approximation of the Jarzynski equality, though seldom invoked, is exact for Gaussian distributions and gives rise to the Fluctuation-Dissipation estimator . Here we derive the parametric maximum-likelihood estimator (MLE) of the free energy considering unidirectional work distributions belonging to Gaussian or Gamma families, and compare this estimator to . We further consider bidirectional work distributions belonging to the same families, and compare the corresponding bidirectional to the Bennett acceptance ratio () estimator. We show that, for Gaussian unidirectional work distributions, is in fact the parametric MLE of the free energy, and as such, the most efficient estimator for this statistical family. We observe that and perform better than and , for unidirectional and bidirectional distributions, respectively. These results illustrate that the characterization of the underlying work distribution permits an optimal use of the Jarzynski equality. © 2018 Wiley Periodicals, Inc.     

Cover Feature: How Different are the Diamagnetic and Paramagnetic Contributions to Off-Nucleus Shielding in Aromatic and Antiaromatic Rings? (ChemPhysChem 9/2023)

The spatial variations in the diamagnetic and paramagnetic contributions to the off-nucleus isotropic shielding, , and to the zz component of the off-nucleus shielding tensor, , around benzene (C₆H₆) and cyclobutadiene (C₄H₄) are investigated using complete-active-space self-consistent field wavefunctions. Despite the substantial differences between and around the aromatic C₆H₆ and the antiaromatic C₄H₄, the diamagnetic and paramagnetic contributions to these quantities, and , and and , are found to behave similarly in the two molecules, shielding and deshielding, respectively, each ring and its surroundings. The different signs of the most popular aromaticity criterion, the nucleus-independent chemical shift (NICS), in C₆H₆ and C₄H₄ are shown to follow from a change in the balance between the respective diamagnetic and paramagnetic contributions. Thus, the different NICS values for antiaromatic and antiaromatic molecules cannot be attributed to differences in the ease of access to excited states only; differences in the electron density, which determines the overall bonding picture, also play an important role.     

A Density Functional Theory Study on Nonlinear Optical Properties of Double Cage Excess Electron Compounds: Theoretically Design M[Cu(Ag)@(NH3)n](M = Be,Mg and Ca; n = 1–3)

In this work, we investigated the nonlinear optical (NLO) properties of excess electron electride molecules of M[Cu(Ag)@(NH₃)_n](M = Be, Mg and Ca; n = 1–3) using density functional theory (DFT). This electride molecules consist of an alkaline-earth (Be, Mg and Ca) together with transition metal (Cu and Ag) doped in NH₃ cluster. The natural population analysis of charge and their highest occupied molecular orbital suggests that the M[Cu(Ag)@(NH₃)_n] compound has excess electron like alkaline-earth metal form double cage electrides molecules, which exhibit a large static first hyperpolarizability () (electron contribution part) and one of which owns a peak value of 216,938 (a.u.) for Be[Ag@(NH₃)₂] and vibrational harmonic first hyperpolarizability () (nuclear contribution part) values and the ratio of /, namely, η values from 0.02 for Be[Ag@(NH₃)] to 0.757 for Mg[Ag@(NH₃)₃]. The electron density contribution in different regions on values mainly come from alkaline-earth and transition metal atoms by first hyperpolarizability density analysis, and also explains the reason why values are positive and negative. Moreover, the frequency-dependent values β(−2ω,ω,ω) are also estimated to make a comparison with experimental measures. © 2018 Wiley Periodicals, Inc.     

Fluorine Substitution of TCNQ Alters the Redox-Driven Catalytic Pathway for the Ferricyanide-Thiosulfate Reaction

Mechanistic variation in catalysis through substituent-based redox tuning is well established. Fluorination of TCNQ (TCNQ=tetracyanoquinodimethane) provides ~850 mV variation in the redox potentials of the and (n=0, 2, 4) processes. With , catalysis of the kinetically very slow ferrocyanide-thiosulfate redox reaction in aqueous solution occurs via a mechanism in which the catalyst is reduced to when reacting with which is oxidised to . Subsequently, reacts with to form and reform the catalyst, in another thermodynamically favoured process. An analogous mechanism applies with as a catalyst. In contrast, since the reaction of with is thermodynamically unfavourable, an alternative mechanism is required to explain the catalytic activity observed in this non-fluorinated system. Here, upon addition of , reduction of to occurs with concomitant oxidation of to , which then acts as the catalyst for oxidation. Thermodynamic data explain the observed differences in the catalytic mechanisms. (n=0, 4) also act as catalysts for the ferricyanide-thiosulfate reaction in aqueous solution. The present study shows that homogeneous pathways are available following addition of these dissolved materials. Previously, these (n=0, 4) coordination polymers have been regarded as insoluble in water and proposed as heterogeneous catalysts for the ferricyanide-thiosulfate reaction. Details and mechanistic differences were established using UV-visible spectrophotometry and cyclic voltammetry.     

Fluxional bonds in quasi-planar and half-sandwich (M = K,Rb, and Cs)

Detailed molecular orbital and bonding analyses reveal the existence of both fluxional σ- and π-bonds in the global minima C_s ( 1 ) and C_s MB₁₈ ( 3 ) and transition states C_s ( 2 ) and C_s ( 4 ) of dianion and monoanions (M = K, Rb, and Cs). It is the fluxional bonds that facilitate the fluxional behaviors of the quasi-planar and half-sandwich which possess energy barriers smaller than the difference of the corresponding zero-point corrections. © 2019 Wiley Periodicals, Inc.     

Supramolecular Assembly of Metal Complexes by (Aryl)I⋅⋅⋅d [PtII] Halogen Bonds

The theoretical data for the half-lantern complexes [{Pt( )(μ- )}₂] [ 1 – 3 ; is cyclometalated 2-Ph-benzothiazole; is 2-SH-pyridine ( 1 ), 2-SH-benzoxazole ( 2 ), 2-SH-tetrafluorobenzothiazole ( 3 )] indicate that the Pt ⋅⋅⋅ Pt orbital interaction increases the nucleophilicity of the outer d orbitals to provide assembly with electrophilic species. Complexes 1 – 3 were co-crystallized with bifunctional halogen bonding (XB) donors to give adducts ( 1 – 3 )₂ ⋅ (1,4-diiodotetrafluorobenzene) and infinite polymeric [ 1⋅ 1,1′-diiodoperfluorodiphenyl]_n. X-ray crystallography revealed that the supramolecular assembly is achieved through (Aryl)I ⋅⋅⋅ d [Pt^II] XBs between iodine σ-holes and lone pairs of the positively charged (Pt^II)₂ centers acting as nucleophilic sites. The polymer includes a curved linear chain ⋅⋅⋅ Pt₂ ⋅⋅⋅ I(arene^F)I ⋅⋅⋅ Pt₂ ⋅⋅⋅ involving XB between iodine atoms of the perfluoroarene linkers and (Pt^II)₂ moieties. The ¹⁹⁵Pt NMR, UV/Vis, and CV studies indicate that XB is preserved in CH(D)₂Cl₂ solutions.     

Photoionization Bands of Cyanogen: Multi-Mode Vibronic Coupling and Renner-Teller Effects

Multi-mode vibronic coupling in the , , and electronic states of Cyanogen radical cation (C N ) is investigated with the aid of ab initio quantum chemistry and first principles quantum dynamics methods. The electronic degenerate states of Π symmetry of C N undergo Renner-Teller (RT) splitting along degenerate vibrational modes of π symmetry. The RT split components form symmetry allowed conical intersections with those from nearby RT split states or with non-degenerate electronic states of Σ symmetry. A parameterized vibronic Hamiltonian is constructed using standard vibronic coupling theory in a diabatic electronic basis and symmetry rules. The parameters of the Hamiltonian are derived from ab initio calculated adiabatic electronic energies. The vibronic spectrum is calculated, assigned and compared with the available experimental data. The impact of various electronic coupling on the vibronic structure of the spectrum is discussed.     

Development and Verification of a Diagnostic Technology for Waste Battery Deterioration Factors

We defined four major deterioration factors (electrolyte loss (EL), lithium loss (LL), lithium precipitation (LP), and compound deterioration (CD)). Then, we derived eleven key performance indicators (KPIs) for comparative analysis. After that, we fabricated three deteriorated cells for each of three deterioration factors (EL, LL, and LP) and one cell with CD (for verification) with four individual (dis)charging experiment manuals. The two major contributions of this study are the performance of 1) trend analysis to determine a suitable diagnostic metric by inspecting the eleven KPIs and 2) comparison analysis of and to verify the effectiveness of utilizing as a real-time deterioration diagnostic factor using a concept of model-in-the-loop simulation. The results show that 1) has the most conspicuous trendline tendency among the eleven comparison targets for all four major deterioration factors, and 2) the angle difference between the two trends of and lies within a minimum of 9° and a maximum of 43° (with a sscale on the x-axis and a scale on the y-axis for a clear trend line analysis). From this, we can conclude that the trendline-based real-time deterioration analysis employing may be practically applicable to a limited extent.     

Complexes Featuring a cis-[M U M] Core (M=Rh,Ir): A New Route to Uranium-Metal Multiple Bonds

Although examples of multiple bonds between actinide elements and main-group elements are quite common, studies of the multiple bonds between actinide elements and transition metals are extremely rare owing to difficulties associated with their synthesis. Here we report the first example of molecular uranium complexes featuring a cis-[M U M] core (M=Rh, Ir), which exhibits an unprecedented arrangement of two M U double dative bond linkages to a single U center. These complexes were prepared by the reactions of chlorine-bridged heterometallic complexes [{U{N(CH₃)(CH₂CH₂NPⁱPr₂)₂}(Cl)₂[(μ-Cl)M(COD)]₂}] (M=Rh, Ir) with MeMgBr or MeLi, a new method for the construction of species with U−M multiple bonds. Theoretical calculations including dispersion confirmed the presence of two U M double dative bonds in these complexes. This study not only enriches the U M multiple bond chemistry, but also provides a new opportunity to explore the bonding of actinide elements.     

Temperature Dependence of 1H Paramagnetic Chemical Shifts in Actinide Complexes,Beyond Bleaney's Theory: The AnVIO22+–Dipicolinic Acid Complexes (An=Np,Pu) as an Example

Actinide +VI complexes ( = , and ) with dipicolinic acid derivatives were synthesized and characterized by powder XRD, SQUID magnetometry and NMR spectroscopy. In addition, and complexes were described by first principles CAS based and two-component spin-restricted DFT methods. The analysis of the ¹H paramagnetic NMR chemical shifts for all protons of the ligands according to the X-rays structures shows that the Fermi contact contribution is negligible in agreement with spin density determined by unrestricted DFT. The magnetic susceptibility tensor is determined by combining SQUID, pNMR shifts and Evans’ method. The SO-RASPT2 results fit well the experimental magnetic susceptibility and pNMR chemical shifts. The role of the counterions in the solid phase is pointed out; their presence impacts the magnetic properties of the complex. The temperature dependence of the pNMR chemical shifts has a strong contribution, contrarily to Bleaney's theory for lanthanide complexes. The fitting of the temperature dependence of the pNMR chemical shifts and SQUID magnetic susceptibility by a two-Kramers-doublet model for the complex and a non-Kramers-doublet model for the complex allows for the experimental evaluation of energy gaps and magnetic moments of the paramagnetic center.     

The Reaction of Sulfur Dioxide Radical Cation with Hydrogen and its Relevance in Solar Geoengineering Models

SO₂ has been proposed in solar geoengineering as a precursor of H₂SO₄ aerosol, a cooling agent active in the stratosphere to contrast climate change. Atmospheric ionization sources can ionize SO₂ into excited states of , which quickly reacts with trace gases in the stratosphere. In this work we explore the reaction of with excited by tunable synchrotron radiation, leading to ( ), where H contributes to O₃ depletion and OH formation. Density Functional Theory and Variational Transition State Theory have been used to investigate the dynamics of the title barrierless and exothermic reaction. The present results suggest that solar geoengineering models should test the reactivity of with major trace gases in the stratosphere, such as H₂ since this is a relevant channel for the OH formation during the nighttime when there is not OH production by sunlight. OH oxides SO₂, triggering the chemical reactions leading to H₂SO₄ aerosol.     

Oxidation of methylamine

A detailed chemical kinetic model for oxidation of methylamine has been developed, based on theoretical work and a critical evaluation of data from the literature. The rate coefficients for the reactions of CHNH + O CHNH / CHNH + HO, CHNH + H CH + NH, CHNH CHNH, and CHNH + O CHNH + HO were calculated from ab initio theory. The mechanism was validated against experimental results from batch reactors, flow reactors, shock tubes, and premixed flames. The model predicts satisfactorily explosion limits for CHNH and its oxidation in a flow reactor. However, oxidation in the presence of nitric oxide, which strongly promotes reaction at lower temperatures, is only described qualitatively. Furthermore, calculated flame speeds are higher than reported experimental values; the model does not capture the inhibiting effect of the NH group in CHNH compared to CH. More work is desirable to confirm the products of the CHNH + NO reaction and to look into possible pathways to NH in methylamine oxidation.     

Potential of AMnO3 (A=Ca,Sr, Ba,La) as Active Layer in Inorganic Perovskite Solar Cells

Inorganic perovskite CaMnO was proposed as a substitution for the TiO anatase in electron transport layers of solar cells containing the hybrid perovskite CH NH PbI based on increased mobility of electrons and better optical matching. Due to a suitable band gap concerning the absorption of sunlight, we investigate the potential of CaMnO and similar manganite perovskites, where Ca is replaced by either Sr, Ba or La, as an absorber layer in inorganic perovskite solar cells. In this study, we have used optical measurements on the synthesized AMnO (A=Ca, Sr, Ba, La) samples to aid density functional theory calculations (DFT) in order to accurately simulate the electronic and optical properties of AMnO compounds and gauge their potential for the role of absorber layer. Both experimental measurements and theoretical calculations show suitable band gap of 1.1-1.5 eV, depending on the compound, and absorption coefficients of the order of cm in the visible part of the spectrum.     

Exponentially correlated Hylleraas-configuration interaction studies of atomic systems. II. Non-relativistic energies of the 1 through 6 states of the Li+ ion

The exponentially correlated Hylleraas-configuration interaction wave function (E-Hy-CI) is a generalization of the Hylleraas-configuration interaction (Hy-CI) in which the single of an Hy-CI wave function is generalized to the generic type . This type of correlation has the right behavior both in the vicinity of the cusp and as goes to infinity; this work shows that wave functions containing both linear and exponential factors converge more rapidly than either one alone for low-lying excited states of ¹S symmetry. E-Hy-CI variational calculations with up to 8568 configurations lead to a nonrelativistic energy of ?7.2799 1341 2669 3059 6491 6759 hartree for the 1 ¹S ground state of the Li+ ion.     

On the normalized Laplacian of Möbius phenylene chain and its applications

Let denote a molecular graph of linear [n] phenylene with n hexagons and n squares, and let the Möbius phenylene chain be the graph obtained from the by identifying the opposite lateral edges in reversed way. Utilizing the decomposition theorem of the normalized Laplacian characteristic polynomial, we study the normalized Laplacian spectrum of , which consists of the eigenvalues of two symmetric matrices ℒ _R and ℒ _Q of order 3n. By investigating the relationship between the roots and coefficients of the characteristic polynomials of the two matrices above, we obtain an explicit closed-form formula of the multiplicative degree-Kirchhoff index as well as the number of spanning trees of . Furthermore, we determine the limited value for the quotient of the multiplicative degree-Kirchhoff index and the Gutman index of .     

Size Evolution of Photoabsorption Spectra of Small Clusters: A Computational Study

Photoabsorption spectra of clusters, N=5–9, have been calculated using a diatomics-in-molecules like electronic structure model and a path-integral Monte Carlo sampling method. A qualitative change in the calculated spectra has been observed at N=9, which has been interpreted in terms of a structural transformation in the clusters consisting in a transition from trimer-like ionic cores observed for N≤7 to dimer-like ionic cores prevailing in through an intermediate state (comparable abundances of both types of ionic cores) observed in . The calculated spectra have been thoroughly compared with an earlier calculation on , , and reported from our group and data available for the same cluster sizes from an experiment.     

Effect of magnetic field on electroosmotic flow of viscoelastic fluids in a microchannel

Electroosmotic flow is an efficient transportation technology driven by applying an external electric field across the microchannel, which has a great potential for future application. This work is presented to study the unsteady electroosmotic flow of viscoelastic fluids combined with a constant pressure gradient and a vertical magnetic field through a parallel plate microchannel. For the reason that the upper and bottom walls of the parallel plate microchannel in microfluidic devices can be made of different materials, this leads to different hydrophobic properties, asymmetric zeta wall potentials, and different slip boundary conditions. The Navier slip model with different slip coefficients at walls is considered. The generalized Maxwell fluid with fractional derivative is adopted for the constitutive equation of the fluid. The analytical and numerical solutions of velocity are derived by employing the integral transform method and finite difference method, respectively. Excellent agreement is found between the numerical solutions and analytical solutions. Finally, the effects of fractional parameter , relaxation time , slip coefficients and , the ratio of wall zeta potentials , Hartmann number , and electrical field strength parameter on velocity profiles are interpreted graphically in detail.     

Analyzing Anisotropic Exchange in a Pentanuclear Os2Ni3 Complex

Spin Hamiltonian parameters of a pentanuclear Os Ni cyanometallate complex are derived from ab initio wave function based calculations, namely valence-type configuration interaction calculations with a complete active space including spin-orbit interaction (CASOCI) in a single-step procedure. While fits of experimental data performed so far could reproduce the data but the resulting parameters were not satisfactory, the parameters derived in the present work reproduce experimental data and at the same time have a reasonable size. The one-centre parameters (local matrices and single-ion zero field splitting tensors) are within an expected range, the anisotropic exchange parameters obtained in this work for an Os−Ni pair are not exceedingly large but determine the low-T part of the experimental χT curve. Exchange interactions (both isotropic and anisotropic) obtained from CASOCI have to be scaled by a factor of 2.5 to obtain agreement with experiment, a known deficiency of such types of calculation. After scaling the parameters, the isotropic Os−Ni exchange coupling constant is cm⁻¹ and the D parameter of the (nearly axial) anisotropic Os−Ni exchange is ⁻¹, so anisotropic exchange is larger in absolute size than isotropic exchange. The negative value of the isotropic J (indicating antiferromagnetic coupling) seemingly contradicts the large-temperature behaviour of the temperature dependent susceptibility curve, but this is caused by the negative g value of the Os centres. This negative g value is a universal feature of a pseudo-octahedral coordination with configuration and strong spin-orbit interaction. Knowing the size of these exchange interactions is important because Os(CN) is a versatile building block for the synthesis of / magnetic materials.     

Reversible Metal-Hydride Transformation in Mg-Ti-H Nanoparticles at Remarkably Low Temperatures

We study the kinetics of hydrogen sorption in Mg-Ti-H nanoparticles prepared by gas phase condensation of mixed Mg-Ti vapors under a H₂-containing atmosphere. Four samples with different Ti contents from 14 to 63 at.% Ti are examined in the 100–150 °C range. The hydrogen absorption kinetics coupled with the formation of MgH₂ can be described by a nucleation and growth model. The activation energy is in the range kJ/mol and the rate constant (at 150 °C) increases from s⁻¹ to s⁻¹ with increasing Ti content. Hydrogen desorption is well modeled by a sequence of surface-limited and contracting-volume kinetics, except at the highest Ti content where nucleation and growth is observed. The activation energy of surface-limited kinetics is /mol. The rate constant (at 150 °C) increases from s⁻¹ to s⁻¹ with the Ti content. These results open an unexplored kinetic window for Mg-based reversible hydrogen storage close to ambient temperature.