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 .     

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.     

Ionized water clusters,n = 2 to 6: A high-accuracy study of structures and energetics

Ionized water clusters, , have been of remarkable interest owing to their crucial roles in many chemical and biological processes. Small cationic water clusters , n = 2 to 6 serve as reasonable models for understanding the nature of the ionized water. In this study, employing high-level ab initio quantum chemical methods, such as the density-fitted orbital-optimized linearized coupled-cluster doubles (DF-OLCCD), coupled-cluster singles and doubles (CCSD), and coupled-cluster singles and doubles with perturbative triples [CCSD(T)], a high-accuracy study of structures and energetics for cationic water clusters [, n = 2-6] is presented. In this study, 2 dimer, 8 trimer, 18 tetramer, 23 pentamer, and 25 hexamer clusters are reported. Most of the structures considered are reported for the first time. Relative, binding, and vertical attachment energies (VAEs), for the first time, are presented at the complete basis set (CBS) limit, extrapolating energies of the aug-cc-pVTZ and aug-cc-pVQZ basis sets, to provide the most accurate energetics to date. Our results demonstrate that as cluster size increases, the VAE value decreases, which indicates that large-size clusters better compensate for the electron deficiency compared with small-size clusters. The VAE values for pentamer and hexamer clusters are 118.5 to 165.5 and 121.9 to 153.7 kcal mol⁻¹, respectively. Further, our binding energy results, at the CCSD(T)/CBS level, indicate strong bindings in cationic clusters due to hydrogen bond interactions. The average binding energy per water molecule varies from −16.6 to −21.8 kcal mol⁻¹ for the clusters considered. Hence, we present the most extensive and accurate study on ionized water clusters to date. Further, our results indicate that the DF-OLCCD method is very promising for ionic molecular clusters, and its accuracy approaches the CCSD(T) quality. The inexpensive analytic gradients of DF-OLCCD, compared with CCSD(T), make it very helpful for high-accuracy studies of molecular geometries.     

On normalized Laplacians,multiplicative degree-Kirchhoff indices,and spanning trees of the linear [n]phenylenes and their dicyclobutadieno derivatives

Let be the molecular graph of the linear [n] phenylene with n hexagons and n − 1 squares, and let be the graph obtained by attaching four-membered rings to the terminal hexagons of . In this article, the normalized Laplacian spectrum of consisting of the eigenvalues of two symmetric tridiagonal matrices of order 3n is determined. An explicit closed-form formula of the multiplicative degree-Kirchhoff index (respectively the number of spanning trees) of is derived. Similarly, explicit closed-form formulas of the multiplicative degree-Kirchhoff index and the number of spanning trees of are obtained. It is interesting to see that the multiplicative degree-Kirchhoff index of (respectively ) is approximately to one half of its Gutman index.     

Is it possible to synthesize bulk salt compounds?

The existence and stability of bulk salt compounds are theoretically investigated in this study. This undertaking is carried out to address the following challenge: synthesizing a bulk salt compound containing a noble gas lighter than Kr. The reliability of theoretical calculations on systems is assessed by benchmark calculations of the well-known salt. In the benchmark calculations, a two-pronged evaluation strategy, including direct and indirect evaluation methods, is used to theoretically investigate the spectroscopic constants of cation and the existence and stability of the salt. The validity of the theoretical calculation methods in the benchmark calculations of salt allows us to adopt a similar methodology to effectively predict the existence and stability of salt compounds. Calculations based on the Born-Haber cycle using estimated lattice energies and some necessary ancillary thermochemical data show that salt compounds can be synthesized, and their upper-limit stable temperatures are estimated to be −237.589, −197.76, and −80.539°C. The salt compound is the most promising candidate. Calculations also show that the salt compounds cannot be stabilized.     

DFT and CIS(D) theoretical study on the ground and excited states of pentanitrogen cation

Different isomers of N₅⁺ were modeled at DFT(PBE0)/aug-cc-pV(Q + d)Z, and their ground(transition) state characteristics were assessed through frequency calculations. Single-point energies were accomplished at PBE0/aug-cc-pV(5 + d)Z. Nonlinear optical susceptibilities (NLO) of isomers were accomplished using Firefly, while the linear optical invariant was examined using the finite-field method, Firefly, and modified dipole field tensor in the presence of two different screening factors. The excited states, singlets and triplets, of were modeled at the CIS and CIS(D) and then their optical parameters were estimated at TDFT(PBE0)/aug-cc-pV(Q + d)Z using Firefly. The singlet is found the most stable isomer, with the inversional rate constant larger than that of the Cs isomer and high energy barrier with the triplet counterpart. Isomers 2 , 3 , and 4 are found local minima, while 5 and 6 are saddle points: transition states between equivalent invertomers. Energy calculations of the singlet and triplet isomers were in excellent agreement with the literature. An excellent correlation is found between the average polarizability and the impulse factor. Substantial variations were found between the singlet and triplet excited states in terms of energy, geometry, and optical properties from one side and with from the other side. Reactivity indices showed that N1 and N5 are the optimum nucleophilic and electrophilic reactivity sites.     

Dynamical Jahn-Teller effect in the first excited

Jahn-Teller effect of C₆₀ monoanion in the first electronically excited states was theoretically investigated. The orbital vibronic coupling parameters for t_1g next lowest unoccupied molecular orbitals were derived from the Kohn-Sham orbital levels calculated using a frozen phonon approach with both hybrid B3LYP and CAM-B3LYP functionals, which take long-range interaction correction into consideration. With these coupling parameters, the vibronic states of first excited were derived by exactly diagonalizing dynamical Jahn-Teller Hamiltonian. The results showed that dynamical Jahn-Teller effects are more significant in the first excited than those in the ground electronic states. This work also clarified that CAM-B3LYP gives results closer to experimental data than B3LYP.     

Extremal Wiener and Kirchhoff indices of globular caterpillars

The Wiener and Kirchhoff indices of a graph G are two of the most important topological indices in mathematical chemistry. A graph G is called to be a globular caterpillar if G is obtained from a complete graph K _s with vertex set {v₁,v₂,…, v _s} by attaching n _i pendent edges to each vertex v _i of K _s for some positive integers s and n₁,n₂,…,n _s, denoted by . Let be the set of globular caterpillars with n vertices (). In this article, we characterize the globular caterpillars with the minimal and maximal Wiener and Kirchhoff indices among , respectively.     

Study on the normalized Laplacian of a penta-graphene with applications

Let L _n denote a linear pentagonal chain with 2n pentagons. The penta-graphene (penta-C), denoted by R _n is the graph obtained from L _n by identifying the opposite lateral edges in an ordered way, whereas the pentagonal Möbius ring is the graph obtained from the L _n by identifying the opposite lateral edges in a reversed way. In this paper, through the decomposition theorem of the normalized Laplacian characteristic polynomial and the relationship between its roots and the coefficients, an explicit closed-form formula of the multiplicative degree-Kirchhoff index (resp. Kemeny's constant, the number of spanning trees) of R _n is obtained. Furthermore, it is interesting to see that the multiplicative degree-Kirchhoff index of R _n is approximately of its Gutman index. Based on our obtained results, all the corresponding results are obtained for .     

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.     

How does the pH influences the Ru-NO coordination compounds?

The pH influence has important role in the bioavailability of coordination compounds. fac-[Ru(NO)Cl₂(κ³N⁴,N⁸,N¹¹[1-carboxypropyl]cyclam)]⁺, 1 , and the species found at different pHs, 2 - 4 , were investigated. One series of computational methodologies has been used to investigate these compounds. One special highlight is to interacting quantum atoms method, where the total interaction energy, , between two atoms has been used as base to estimate the chemical bonds strength. The deprotonation of -CO₂H, 1 ➔ 2 (pK_a = 3.3), creates a hydrogen bond in the complex 2 , N( 3 )-H⋯ ·OCO⁻, with a more favorable than the presents in 1 , N( 3 )-H⋯ ·OCOH. There are no changes in in Ru-NO bond. The second deprotonation occurs in the N(2) atom of the cyclam group, 2 ➔ 3 (pK_a = 8.0). It promotes an increase in the covalent character of Ru-N( 2 ). In contrast, there is no changes in Ru-N( 5 )O bond. For higher pHs, there is a 3 ➔ 4 equilibrium (pK_a = 11.5) and the conversion of Ru-N( 5 )O for Ru-N( 5 )O₂. The Ru-N( 5 ) of 4 shows a larger ionic character than 3 . Thus, Ru-NO in 1 - 4 has worthy stability about a large pH range, showing potential application as NO scavengers.     

Virial relations in density embedding

The accuracy of charge-transfer excitation energies, solvatochromic shifts, and other environmental effects calculated via various density-embedding techniques depend critically on the approximations employed for the nonadditive noninteracting kinetic energy functional, . Approximating this functional remains an important challenge in electronic-structure theory. To assist in the development and testing of approximations for , we derive two virial relations for fragments in molecules. These establish separate connections between the nonadditive kinetic energies of the noninteracting and interacting systems of electrons, and quantities such as the electron-nuclear attraction forces, the partition (or embedding) energy and potential, and the Kohn-Sham potentials of the system and its parts. We numerically verify both relations on diatomic molecules.     

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.     

Comparison of triel bonds with different chalcogen electron donors: Its dependence on triel donor and methyl substitution

The complexes between R₃Tr (Tr = B, Al, and Ga; R = H, F, Cl, and Br) and H₂X (X = O, S, and Se) were theoretically studied. The interaction energies of R₃Al⋯H₂X and R₃Ga⋯H₂X are consistent with the electronegativity of the halogen atom R (R ≠ H), but an opposite dependence is found for R₃B⋯H₂X. The triel bond of R₃Tr⋯H₂X is weaker for the heavier chalcogen donor. The dependence of triel bonding strength on the triel atom is complicated, depending on the nature of R and X. The methyl substitution of H₂X causes a substantial increase in the interaction energy from −5.74 kcal/mol to −22.88 kcal/mol, and its effect is relevant to the nature of Tr, X, and R groups. For the S and Se donors, the increased percentage of interaction energy is almost the same due to the methyl substitution, which is larger than that of the O analogue. In most triel-bonded complexes, electrostatic dominates and polarization has comparable contribution. However, polarization plays a dominant role in R₃B⋯ and R₃B⋯ (R = Cl and Br; R′ = H and Me).     

Molecular Oxygen Trimer: Multiplet Structures and Stability

We present a detailed theoretical study of the molecular oxygen trimer where the potential energy surfaces of the seven multiplet states have been calculated by means of a pair approximation with very accurate dimer ab initio potentials. In order to obtain all the states a matrix representation of the potential using the uncoupled spin representation has been applied. The and states are nearly degenerate and low-lying isomers appear for most multiplicities. A crucial point in deciding the relative stabilities is the zero-point energy which represents a sizable fraction of the electronic well-depth. Therefore, we have performed accurate diffusion Monte Carlo studies of the lowest state in each multiplicity. Analysis of the wavefunction allows a deeper interpretation of the cluster structures, finding that they are significantly floppy in most cases.     

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.     

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.     

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.