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.     

Theoretical study of metal‐free organic dyes based on different configurations for efficient dye‐sensitized solar cells

Considering different solar dyes configuration, four novel metal‐free organic dyes based on phenoxazine as electron donor, thiophene and cyanovinylene linkers as the ‐conjugation bridge and cyanoacrylic acid as electron acceptor were designed to optimize open circuit voltage and short circuit current parameters and theoretically inspected. Density functional theory and time‐dependent density functional theory calculations were used to study frontier molecular orbital energy states of the dyes and their optical absorption spectra. The results indicated that D2‐4 dyes can be suitable candidates as sensitizers for application in dye sensitized solar cells and among these three dyes, D3 showed a broader and more bathochromically shifted absorption band compared to the others. The dye also showed the highest molar extinction coefficient. This work suggests optimizing the configuration of metal‐free organic dyes based on simple D‐ ‐A configuration containing alkyl chain as substitution, starburst conformation, and symmetric double D‐ ‐A chains would produce good photovoltaic properties.     

Nonstationarity and related measures for time‐dependent hartree–fock and multiconfigurational models

Based on an earlier article (Eberly and Singh, Phys. Rev. D 1973 , 7, 359) and related works on short‐time evolution, this article proposes a many‐electron formulation for the nonstationarity degree which can be assigned to quantum system at each time point. The key measure introduced, , is a nonstationarity index that can be thought of as an inverse nominal lifetime at each instance of time. The index is directly computed from the time derivative of one‐electron density matrix and is a size‐consistent quantity. In this article, the approach is developed for the time‐dependent Hartree–Fock (TDHF), single‐excitation (TDCIS), and time‐dependent full configuration interaction (TDFCI) models. As a rule, nonstationarity effects are more pronounced in correlated electron systems, and a joint analysis of and the multiconfigurational character of wave functions apparently provide a deeper insight into dynamical molecular processes. The performed calculations on small molecules in laser fields show a preference for the TDCIS model when comparing TDCIS and TDHF with the “exact” TDFCI model. © 2013 Wiley Periodicals, Inc.     

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.     

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.     

Critical parameters and spherical confinement of H atom in screened Coulomb potential

Critical parameters in three screened potentials, namely, Hulthén, Yukawa, and exponential cosine screened Coulomb potential are reported. Accurate estimates of these parameters are given for each of these potentials, for all states having . Comparison with literature results is made, wherever possible. Present values compare excellently with reference values; for higher n, ?, our results are slightly better. Some of these are presented for first time. Further, we investigate the spherical confinement of H atom embedded in a dense plasma modeled by an exponential cosine screened potential. Accurate energies along with their variation with respect to box size and screening parameter are calculated and compared with reference results in literature. Sample dipole polarizabilities are also provided in this case. The generalized pseudospectral method is used for accurate determination of eigenvalues and eigenfunctions for all calculations. © 2016 Wiley Periodicals, Inc.     

Ro‐vibrational studies of diatomic molecules in a shifted Deng–Fan oscillator potential

Bound‐state spectra of shifted Deng‐Fan oscillator potential are studied by means of a generalized pseudospectral method. Very accurate results are obtained for both low as well as high states by a nonuniform optimal discretization of the radial Schrödinger equation. Excellent agreement with literature data is observed in both s‐wave and rotational states. Detailed variation of energies with respect to potential parameters is discussed. Application is made to the ro‐vibrational levels of representative diatomic molecules (H₂, LiH, HCl, CO). Nine states having are calculated with good accuracy along with 15 other higher states for each of these molecules. Variation of energies with respect to state indices n, show behavior similar to that in the Morse potential. Many new states are reported here for the first time. In short, a simple, accurate, and efficient method is presented for this and other similar potentials in molecular physics. © 2013 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.     

Design of neutral organic superacids using fulvene derivatives with di‐enol substituent

A new class of superacids was designed using enolic derivatives of fulvene. After deprotonation, bond rearrangement leads in a stable conjugate base with an aromatic cyclopentadienyl ring and a carboxyl group. The gas phase enthalpies ( ) of the deprotonation, as an index of acidity, were calculated employing the B3LYP method and 6‐311++G(d,p) and aug‐cc‐PVDZ basis sets. The acidity of these compounds without any electron withdrawing groups was more than H₂SO₄ in gas phase. The acidity increased by substituting electron withdrawing groups (? F, ? CN, and ?O) into the molecules so that we could achieve a cyano derivative of fulvene with = 250 kcal/mol.     

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.     

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.     

Dealing with the shifted and inverted Tietz–Hua oscillator potential using the J‐matrix method

The tridiagonal J‐matrix approach has been used to calculate the low and moderately high‐lying eigenvalues of the rotating shifted Tietz–Hua (RSTH) oscillator potential. The radial Schrödinger equation is solved efficiently by means of the diagonalization of the full Hamiltonian matrix, with the Laguerre or oscillator basis. Ro–vibrational bound state energies for 11 diatomic systems, namely , , , NO, CO, , , , , , and NO⁺, are calculated with high accuracy. Some of the energy states for molecules are reported here for the first time. The results of the last four molecules have been introduced for the first time using the oscillator basis. Higher accuracy is achieved by calculating the energy corresponding to the poles of the S‐matrix in the complex energy plane using the J‐matrix method. Furthermore, the bound states and the resonance energies for the newly proposed inverted Tietz–Hua IRSTH‐potential are calculated for the H₂‐molecule with scaled depth. A detailed analysis of variation of eigenvalues with n, quantum numbers is made. Results are compared with literature data, wherever possible. © 2015 Wiley Periodicals, Inc.     

Ritz variational method for the high‐lying nonautoionizing doubly excited 1,3Fe states of two‐electron atoms

Energy eigenvalues of nonautoionizing doubly excited states originating from 2pnf ( ) configuration of two‐electron atoms have been calculated by expanding the basis set in explicitly correlated Hylleraas coordinates under the framework of Ritz variational method. A detailed discussion on the evaluation of correlated basis integrals is given. The energy eigenvalues of a number of these doubly excited states are being reported for the first time especially for the high lying states. The effective quantum numbers ( ) for the states mentioned above have been calculated by using the theory of quantum defect.     

Improved description of the orbital relaxation effect by practical use of the self‐energy

The self‐energy shift in the orbital relaxation (OR) term of the polarization propagator complete through the second‐order is presented. In combination with the optimal damping parameter in the OR term, the modified propagator produces the excitation energy of the coupled‐cluster with singles and doubles (CCSD) accuracy. The self‐energy shift requires the floating‐point operation of , where N refers to the magnitude of the molecular size. Because the second‐order polarization propagator requires the floating‐point operation of , the additional computational effort to construct the self‐energy is negligibly small. Numerical results are shown for several molecules including glycine, 2,3,5,6‐tetrafluorobenzene, and naphthalene, and promising agreements with those of CCSD are confirmed within less than 0.2 eV. The basis set dependence is also tested for the water molecule using aug‐cc‐pV NZ (N = D–7), where this newly developed approach mimics the behavior of the CCSD values. The self‐energy shifting for the second‐order response matrix in combination with the use of a dumping parameter is efficiently implemented for calculations of medium‐sized molecular systems, including glycine and naphthalene. The developed approach provides CCSD‐like accuracy at a more affordable computational expense. © 2014 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.     

Photodetachment of hydrogen negative ion near inelastic surfaces: Arbitrary laser polarization direction

The photodetachment of hydrogen negative ion near different inelastic surfaces is investigated by the semiclassical closed orbit theory for arbitrary laser polarization direction . A two‐term formula of photodetachment cross section consisting of a smooth background term and an oscillatory term is derived. The oscillatory term contains an extra angular factor that describes the dependence of oscillations in total cross section on the laser polarization direction. It is observed that the amplitude of oscillations in cross section reaches maximum at when laser polarization is parallel to the z‐axis and it approaches zero as the laser polarization direction becomes perpendicular to the z‐axis. It is also observed that as the reflection coefficient , which accounts for the inelastic behavior of the surfaces, increases the amplitude of oscillation also increases. © 2015 Wiley Periodicals, Inc.     

Electron‐density delocalization in many‐electron atoms confined by penetrable walls: A Hartree–Fock study

The electronic structure of several many‐electron atoms, confined within a penetrable spherical box, was studied using the Hartree–Fock (HF) method, coupling the Roothaan's approach with a new basis set to solve the corresponding one‐electron equations. The resulting HF wave‐function was employed to evaluate the Shannon entropy, , in configuration space. Confinements imposed by impenetrable walls induce decrements on when the confinement radius, R_c, is reduced and the electron‐density is localized. For confinements commanded by penetrable walls, exhibits an entirely different behavior, because when an atom starts to be confined, delivers values less than those observed for the free system, in the same way that the results presented by impenetrable walls. However, from a confinement radius, shows increments, and precisely in these regions, the spatial restrictions spread to the electron density. Thus, from results presented in this work, the Shannon entropy can be used as a tool to measure the electron density delocalization for many‐electron atoms, as the hydrogen atom confined in similar conditions.     

A DFT study of reactions of methyldiazonium ion with DNA/RNA nucleosides: Investigating effect of sugar moiety on methylation pattern of bases

Methyldiazonium ion ( ) is an ultimate carcinogen that can methylate multiple sites in DNA/RNA. In present contribution, density functional theory calculations using the B3LYP and M06‐2X functionals and the 6‐31G(d,p) and aug‐cc‐pVDZ basis sets are carried out to study methylation reactions of at the different nucleophilic sites of DNA/RNA bases and their nucleosides. Total 12 nucleophilic sites, that is, the N2, N3, N7, and O6 sites of guanine; the N1, N3, N6, and N7 sites of adenine; O2 and N3 sites of cytosine and the O2 and O4 sites of thymine and uracil have been considered for study. Thus, a total of 30 reactions have been studied here. The polarizable continuum model is used for solvation calculations. The N7 site of guanine, N7(G), is found to be most reactive in all the reactions studied here, which is in agreement with experiment. However, the calculated reactivity of toward the N7(G) site in aqueous media follows the order: guanine > deoxyguanosine > guanosine. The reactivities of many other sites including the O6(G), O2(C), and N3(A) sites are also modified in going from DNA/RNA bases to their nucleosides and from DNA to RNA nucleosides. Thus, we note that the presence of sugar moiety significantly modifies the methylation pattern of bases caused by . © 2014 Wiley Periodicals, Inc.     

Presentation of a new index for estimation of aromaticity in halo‐ and cyanobenzenes: The role of potential energy in aromaticity

A linear correlation has been obtained between average values of Hamiltonian kinetic energy ( ) and potential energy ( ) calculated at the bond critical points using atoms in molecules method. This relation was used to introduce a new index ( ) for estimation of aromaticity in halo‐ and cyanobenzenes. Potential energy has different terms such as attraction between nuclei and electrons, also repulsion of electrons which affect the inertia and mobility of electrons, respectively. Therefore, contribution of potential energy in this relation must be controlled. Contribution of potential energy in aromaticity has been managed using a fitting parameter. This parameter was obtained by fitting the aromaticity stabilization energy data with values of aromaticity calculated by index for halo‐ and cyanobenzenes. The contribution of potential energy in index is complete when molecule is nonaromatic and is negligible when molecule is antiaromatic. Indeed, molecule is aromatic when contribution of potential energy in index lies between above limits. © 2013 Wiley Periodicals, Inc.