Quantum control of electron‐proton symmetry breaking in dissociative ionization of H2 by intense laser pulses

Forward and backward electron/proton ionization/dissociation spectra from one‐dimensional non‐Born‐Oppenheimer H₂ molecule exposed to ultrashort intense laser pulses ( W/cm², λ = 800 nm) have been computed by numerically solving the time‐dependent Schrödinger equation. The resulting above‐threshold ionization and above‐threshold dissociation spectra exhibit the characteristic forward‐backward asymmetry and sensitivity to the carrier‐envelope phase (CEP), particularly for high energies. A general framework for understanding CEP effects in the asymmetry of dissociative ionization of H₂ has been established. It is found that the symmetry breaking of electron‐proton distribution with π periodic modulation occurs for all CEPs except for ( integer) and the largest asymmetry coming from the CEP of . At least one of the electron and proton distributions is asymmetric when measured simultaneously. Inspection of the nuclear and electron wave packet dynamics provides further information about the relative contribution of the gerade and ungerade states of to the dissociation channel and the time delay of electrons in asymmetric ionization. © 2014 Wiley Periodicals, Inc.     

Theoretical prediction of the optical rotation of chiral molecules in ordered media: A computational study of (Ra)‐1,3‐dimethylallene, (2R)‐2‐methyloxirane,and (2R)‐N‐methyloxaziridine

A theoretical procedure has been developed and implemented to calculate the optical rotation of chiral molecules in ordered phase via origin‐independent diagonal components , of the optical activity tensor and origin‐independent components , for , of the mixed electric dipole‐electric quadrupole polarizability. Origin independence was achieved by referring these tensors to the principal axis system of the electric dipole dynamic polarizability at the same laser frequency ω. The approach has been applied, allowing for alternative quantum mechanical methods based on different gauges, to estimate near Hartree–Fock values for three chiral molecules, (2R)‐N‐methyloxaziridine C₂NOH₅, (2R)‐2‐methyloxirane (also referred to as propylene oxide) C₃OH₆, and ( )‐1,3‐dimethylallene C₅H₈, at two frequencies. The theoretical predictions can be useful for an attempt at measuring correspondent experimental values in crystal phase. © 2015 Wiley Periodicals, Inc.     

Additivity in kramers pairs symmetry: Multiplets with up to four unpaired electrons

Many fermions Kramers pairs formalism is considered from the prospective of the sum of individual single fermion time‐reversal operators. The obtained many fermions “pseudo Kramers pairs operator” ( ), as well as its square ( ), have formally the same structure as the many fermion spin operators and . Nevertheless, the shape of eigenfunctions with respect to and is different. Herein all Kramers adapted eigenfunctions of for cases of up to four unpaired fermions are compiled, and their properties with respect to further advocated. It will be shown that degeneracy of the multiplets recovers the proper behavior with respect to Pascal's triangle. A projection operator for obtaining the “high spin” Kramers adapted eigenfunctions is suggested. Noncommutation of with spin and angular momentum operators and degeneracy is discussed at last. © 2016 Wiley Periodicals, Inc.     

The electron delocalization range in stretched bonds

The electron delocalization range function EDR( (Janesko et al., J. Chem. Phys. 2014, 141, 144104) quantifies the width of the one‐particle density matrix about point , measuring aspects of delocalization. Here, we explore the EDR in stretched and compressed chemical bonds. The EDR illustrates how compressing chemical bonds localizes the one‐particle density matrix about the reference point, and captures aspects of fractional occupancy, and left‐right correlation in stretched covalent bonds.     

(Rg = He ∼ Rn,n = 1–4): In quest of the potential trapping ability of the aromatic ring

A new series of divalent boron‐rare gas cations (Rg = He ∼ Rn, n = 1–4) have been predicted theoretically at the B3LYP, MP2, and CCSD(T) levels to present the structures, stability, charge distributions, bond natures, and aromaticity. The Rg B bond energies are quite large for heavy rare gases and increase with the size of the Rg atom. Because of steric hindrance new Rg atoms introduced to the B₄ ring will weaken the Rg B bond. Thus in the Rg B bond has the largest binding energy 90–100 kcal/mol. p‐ has a slightly shorter Rg B bond length and a larger bond energy than o‐ . NBO and AIM analyses indicate that for the heavy Rg atoms Ar ∼ Rn the B Rg bonds have character of typical covalent bonds. The energy decomposition analysis shows that the σ‐donation from rare gases to the boron ring is the major contribution to the Rg B bonding. Adaptive natural density partitioning and nuclear‐independent chemical shift analyses suggest that both and have obvious aromaticity.     

General build up of basis and matrix in the diagonalization approach. Determination of Kramers configuration state functions

Algorithms to build the basis and matrix representation to obtain the Kramers configuration space functions (KCSFs) via diagonalization will be formally generalized to an arbitrary number of unpaired (open shell) fermions. Effective build up of the matrix representation will be outlined (including threading and graphical processing unit parallelism) to subsequently obtain the KCSFs via calling external/numerical library routines for diagonalization. The effective build up of the matrix representation relays on a binary tree search algorithm to allow evaluation the action on a given basis vector. The binary tree search avoids the treatment of zero matrix elements which leads to an exponential acceleration. The implementation ( basis creation, matrix representation, and matrix diagonalization) will be done in an all in core and all at once manner, hence the available core memory sets the physical limits in practical applications. Memory limitations, sparsity of the matrix, general case of n fermions in m spinors, and the application of KCSFs will be put into further perspective.     

A theoretical study on the stability difference of the borane and carborane C2Bn−2Hn (5 ≤ n ≤ 7) clusters

In order to study the electronic structure and structural stability of borane and carborane C₂B_n?2H_n (5 ≤ n ≤ 7) clusters, especially the stability difference between the borane and carborane C₂B₃H₅. The frontier orbital energy levels of the borane and carborane C₂B_n?2H_n (5 ≤ n ≤ 7) clusters are calculated at CCSD(T)/aug‐cc‐pVXZ//B3LYP/def2‐TZVPP level. The results are further analyzed by qualitative frontier orbital method based on the cap–ring interaction. The results reveal that: (1) the larger E_gap(HOMO‐LUMO energy gap) of carborane C₂B_n?2H_n (5 ≤ n ≤ 7) clusters than borane (5 ≤ n ≤ 7) clusters originates from the more effective cap–ring orbital overlap of carborane C₂B_n?2H_n (5 ≤ n ≤ 7) clusters than that of borane (5 ≤ n ≤ 7) clusters; (2) the smallest E_gap of the borane results from the highest energy level of the ring symmetry‐adapted linear combination orbital of cluster; and (3) the largest E_gap of the carborane C₂B₃H₅ is induced by the most effective cap–ring orbital interaction of C₂B₃H₅ cluster. © 2014 Wiley Periodicals, Inc.     

Carbon–sulfur chains: A high‐resolution infrared and quantum‐chemical study of C3S and SC7S

In the course of a 5 μm high‐resolution infrared study of laser ablation products from carbon–sulfur targets, the ν₁ vibrational mode region of linear C₃S has been studied continuously from 2046 to 2065 cm^?1. Besides the prominent vibrational fundamental, the region was found to feature the , and even hot bands, the latter two of which were observed for the first time. Owing to the high signal‐to‐noise ratio obtained, the ν₁ mode of S could also be observed in natural abundance for the first time at high spectral resolution in the infrared. At 2061 cm^?1, hidden inside the branch of the C₃S ν₁ fundamental mode, a weak new band is observed which exhibits very tight line spacing and stems from a heavy both carbon and sulfur containing carrier. On the basis of high‐level quantum‐chemical calculations of selected carbon–sulfur chains and other carbon‐rich cumulenes, this feature is attributed to the ν₅ vibrational fundamental of linear SC₇S, which stands for the first gas‐phase spectroscopic detection of this long cumulenic chain.     

Einstein A Coefficients,Oscillator Strengths,and Absolute Band Strengths for the First Negative System,and Franck-Condon Factors for the Second Negative System of

The Franck-Condon factors and r centroids for the first negative and the second negative band systems of the molecule, based on the Rydberg-Klein-Rees potential energy curves, have been computed. The variation of the electronic transition moment with the internuclear separation has been studied for the first negative bands and the Einstein A coefficients, the oscillator strengths, and the absolute band strengths for this system have been calculated by adopting the recent experimental data on the lifetimes of the various levels of the excited B state. On a calculé pour la molécule les facteurs de Franck-Condon et les r centroïds pour les deux premiers systèmes de bandes négatives avec des courbes de potential de Rydberg-Klein-Rees. On a étudié la variation du moment de transition éléctronique avec la séparation internucléaire pour les premières bandes négatives. Les constantes A d'Einstein, les forces d'oscillateur et les intensités absolues des bandes ont été calculé des données expérimentales récentes sur les durées de vie des niveaux différents de l'état excité B . Die Franck-Condon-Faktoren und die r-Zentroide der zwei ersten negativen Banden-systeme des -Moleküls, wurden mit den Rydberg-Klein-kees Potentialkurven berechnet. Die Variation des elektronischen Übergangsmoment mit dem Kernabstand würde für die ersten negativen Bänder studiert. Die Einsteinschen A-Koeffizienten, die Oszillatorstärken and die absoluten Bandenstärken wurden mit die neuen experimentellen Tatsachen über die Lebensdauern der verschiedenen Niveaus des angeregten B Zustands berechnet.     

The nature of the chemical bond in borazine (B3N3H6), boroxine (B3O3H3), carborazine (B2N2C2H6), and related species

The bonding problem in borazine (B₃N₃H₆), boroxine (B₃O₃H₃), and carborazine (B₂N₂C₂H₆) is successfully addressed through the consideration of the excited states of the constituent fragments, namely BH( ), NH( ), and CH( ). We propose the participation of resonant structures for all three species that help to explain the experimental findings. A discussion on the chemical pattern of the parental molecule benzene (C₆H₆) helps to make coherent the whole bonding analysis on the titled species.     

Thermodynamic properties of diatomic molecule systems under SO(2,1)‐anharmonic Eckart potential

Due to one of the most representative contributions to the energy in diatomic molecules being the vibrational, we consider the generalized Morse potential (GMP) as a typical interaction for one‐dimensional microscopic systems, which describes local anharmonic effects. From the Eckart potential (EP) model, it is possible to find a connection with the GMP model, as well as obtaining the analytical expression for the energy spectrum because it is based on algebras. This gives the macroscopic properties such as vibrational mean energy U, specific heat C, Helmholtz free energy F, and entropy S for a heteronuclear diatomic system, as well as with the exact partition function and its approximation for the high temperature region. Finally, a comparison is made between the graphs of some thermodynamic functions obtained with the GMP and the Morse potential (MP) for molecules.     

Rovibrational spectra of bounded diatomic molecules

Spectra of a bounded diatomic molecule is studied numerically. Shifted Deng–Fan oscillator potential has been used to model the molecule. The accurate five‐point finite difference method has been used to solve the Schrödinger equation for rovibrational motion of the molecule. The energies of the bound states as well as free states of the molecule have been calculated. In addition, radial matrix elements like , n = 1, 2, and 3 have been calculated. These have been used to calculate the ‐pole static polarizabilities. The variation of bound state energies, matrix elements and ‐pole static polarizabilities with the boundary radius has also been studied. The Stark effect in case of this bounded system has also been investigated.     

Full–dimensional quantum molecular dynamics calculations of the rovibrationally mediated X → 2 transition of nitrous oxide

A full dimensional time‐dependent quantum wavepacket approach is used to study the photodissociation dynamics of nitrous oxide for the X → 2 bound–bound transition based on new highly accurate potential energy and transition dipole moment surfaces. The computed 2 absorption spectra at room temperature are characterized by sharp vibrational structures that contribute slightly to the diffuse vibrational structures around the maximum peak at 180 nm of the first ultraviolet absorption band (from the contribution of 2 , 1 , and 2 states) of N₂O. Transitions from different initial rovibrational states reveal that the sharp structures arise mainly from N₂? O bending vibrations, whereas, at higher temperatures, the N₂? O and N? NO stretching vibrations are responsible for enhancing the intensity of the structures. At absorption wavelengths 166 nm and 179 nm, vibrational quantum state distributions of N₂ product fragments decrease monotonically with increasing vibrational quantum number v = 0, 1, 2. At 166 nm, rotational quantum state distributions of N₂ at fixed v = 0 and v = 1 display multimodal profiles with maximum peaks at j = 77 and j = 75, respectively, whereas, the distributions at the 179 nm absorption wavelength display bimodal profiles with maximum peaks at j = 73 and j = 71, respectively. Accordingly, the presence of rotationally hot N₂ from previous experimental and theoretical works in the first band strongly implies a significant influence of the 2 state in determining the final dissociation pathway of N₂ + O. © 2016 Wiley Periodicals, Inc.     

Ritz variational calculation for the singly excited states of compressed two‐electron atoms

A detailed analysis on the effect of spherical impenetrable confinement on the structural properties of two‐electron ions in states has been performed. The energy values of 1sns [ ] ( ) states of helium‐like ions ( ) are estimated within the framework of Ritz variational method using explicitly correlated Hylleraas‐type basis sets. The correlated wave functions used here are consistent with the finite boundary conditions due to spherical confinement. A comparative study between the singlet and triplet states originating from a particular electronic configuration shows incidental degeneracy and the subsequent level‐crossing phenomena. The thermodynamic pressure felt by the ion inside the sphere pushes the energy levels toward continuum. Critical pressures for the transition to strong confinement regime (where the singly excited two‐electron energy levels cross the corresponding one‐electron threshold) as well as for the complete destabilization are also estimated.     

Study of adiabatic connection in density functional theory with an accurate wavefunction for two‐electron spherical systems

An accurate semianalytic wavefunction is proposed for the Hookium and two‐electron atoms for varying strength of where is the strength parameter and is coulomb interaction between two electrons. The wavefunction leads to energies that are as accurate as those from the Coupled cluster singles and doubles (CCSD) calculations. Using this wavefunction, we construct the external potential such that the density of the system remains unchanged as is varied. The work thus gives a unified picture of adiabatic connection for these systems based on an easy to use wavefunction and complements the past investigations done in this direction. Using the potential obtained, we explicitly calculate the energy of the corresponding positive ions and show that the chemical potential—calculated as the difference between the energies of the two‐electron system and its positive ion—is equal to the experimental ionization energy and remains unchanged as is varied. Furthermore, using total energies of these systems as a function of , we provide a new perspective into a variety of hybrid functionals.     

Valence bonds in elongated boron clusters

A well‐defined class of planar or quasi‐planar elongated boron clusters, of type , serves as a basis to identify the valence bond picture of delocalized boron networks. The origin of the series is the cluster, which exhibits σ‐aromaticity. The cluster generating step is the repetitive expansion by three boron atoms in the direction of elongation. Specific electron counting rules are obtained for π‐bonding, peripheral σ‐bonding and multicenter inner σ‐bonding. A valence bond structure is introduced which explains the remarkable regularity in the bonding pattern. The analysis supports 4c‐2e bonds as an alternative to the common 3c‐2e bonds. The results are validated by symmetry induction and ab initio calculations.     

Geometry of the canonical Van Vleck transformation

Bloch's transformation from the zeroth‐order space for a perturbation problem to the corresponding space of exact eigenvectors, was found as a geometrically defined alternative to the algebraically constructed Van Vleck transformation. Klein's theorem of uniqueness transferred some of this geometrical interpretation to its canonical form . Quite recently Kvaal has taken a large step further by writing as a product of commuting planar rotations, obtained by describing and in terms of certain principal vectors and canonical angles. Kvaal's approach is now developed further, using a new commutation relation which simplifies algebraic manipulations substantially. It allows for a simple definition of an operator for the angle between and which has Kvaal's vectors and angles as eigenvectors and eigenvalues. Klein's theorem is refined in various ways. The impact of the approach on a number of previous results is considered. © 2015 Wiley Periodicals, Inc.     

A theoretical study on cyclacenes: Analytical tight‐binding approach

We present a theoretical study of cyclacene molecules performed at tight‐binding level. The orbital energies and eigenvectors have been analytically computed, and exact expressions for the axial component of the total position spread and polarizability tensors have been obtained. In absence of dimerization, the system has a D_nh symmetry, where n is the number of hexagonal units. The energy bands present no gap at the Fermi level, and to this fact it corresponds a diverging (per‐electron) polarizability for in the direction of the system symmetry axis. The two (degenerate) components of the polarizability on the σ_h symmetry plane, conversely, remain finite for . The total position spread tensor presents a qualitatively different behavior, since all the three components of the position spread per electron remain finite for . The results are analyzed and discussed for both axial and planar components separately as these are affected differently with respect to the increasing system size. Both dipole polarizability and total position spread have been computed using an ab initio approach for the smallest systems, to compare the analytical tight‐binding expressions with a higher‐level theory.     

Structural,magnetic, and vibrational properties of stoichiometric clusters of CrN

We performed density‐functional‐theoretic calculations to investigate the structural, magnetic and vibrational properties of the stoichiometric clusters (CrN)_n ( ). We show that the building block of the ground‐state structures of these clusters is a square CrNCrN unit; the only exception with n > 2 occurs for (CrN)₃, but this cluster has an isomer not very far in energy from the ground state consisting of a pair of CrNCrN squares sharing a CrN bond. In the smaller CrN, (CrN)₂, and (CrN)₃ clusters the magnetic moments of the N atoms are non‐negligible and antiparallel to those of the Cr atoms, but for the larger species (CrN)₄, (CrN)₅, (CrN)₆, and (CrN)₉ the cluster magnetic moments are almost entirely due to the Cr atoms. Lack of imaginary vibrational frequencies in the predicted ground‐state structures of (CrN)_n ( ) confirms that they are mechanically stable equilibrium states. © 2015 Wiley Periodicals, Inc.     

3D‐QSAR studies and molecular design on a novel series of pyrimidine benzimidazoles as Lck inhibitors

The development of selective lymphocyte‐specific kinase (Lck) inhibitors has attracted much attention for the research of the treatment of T‐cell mediated autoimmune and inflammatory diseases. In the present work, three‐dimensional quantitative structure–activity relationship (3D‐QSAR) analyses are performed on a novel series of 4‐amino‐6‐benzimidazole‐pyrimidines acting as Lck inhibitors. The established 3D‐QSAR models show significant statistical quality and satisfactory predictive ability, with high q² and R² values: the comparative molecular field analysis (CoMFA) model (q² = 0.802, R² = 0.991), and the comparative molecular similarity indexes analysis (CoMSIA) model (q² = 0.731, R² = 0.982). The systemic external validation indicates that both CoMFA and CoMSIA models are quite robust and possess high predictive abilities with values of 0.881 and 0.877, values of 0.897 and 0.847, values of 0.897 and 0.850, and values of 0.897 and 0.854, respectively. Several key structural features accounting for the inhibitory activities of these compounds are discussed. Based on established models and design considerations, six new compounds with significantly improved activities are theoretically designed, which still await experimental confirmation and evaluation. These theoretical results may provide a useful reference for understanding the action mechanism and designing novel potential Lck inhibitors. © 2014 Wiley Periodicals, Inc.