Quadratic response functions in the relativistic four-component Kohn-Sham approximation

Frequency-Dependent Quadratic Response Properties and Two-Photon Absorption from Relativistic Equation-of-Motion Coupled Cluster Theory

We present the implementation of quadratic response theory based upon the relativistic equation-of-motion coupled cluster method. We showcase our implementation, whose generality allows us to consider both time-dependent and time-independent electric and magnetic perturbations, by considering the static and frequency-dependent hyperpolarizability of hydrogen halides (HX, X = F-At), providing comprehensive insights into their electronic response characteristics. Additionally, we evaluated the Verdet constant for noble gases Xe and Rn and discussed the relative importance of relativistic and electron correlation effects for these magneto-optical properties. Finally, we calculate the two-photon absorption cross sections of transition [ns1S0 → (n + 1)s1S0] of Ga+ and In+, which are suggested as candidates for new ion clocks. As our implementation allows for the use of nonrelativistic Hamiltonians as well, we have compared our EOM-QRCC results to the QR-CC implementation in the DALTON code and show that the differences between CC and EOMCC response are in general smaller than 5% for the properties considered. Collectively, the results underscore the versatility of our implementation and its potential as a benchmark tool for other approximated models, such as density functional theory for higher-order properties.

Exploring quantum computational, molecular docking, and molecular dynamics simulation with MMGBSA studies of ethyl-2-amino-4-methyl thiophene-3-carboxylate

2-aminothiophenes derivative, Ethyl-2-amino-4-methyl thiophene-3-carboxylate (EAMC) has been synthesized, characterized, and investigated quantum chemically. It was experimentally investigated by different spectroscopic methods like- NMR (1H-NMR and 13C-NMR), FT-IR, and UV-Visible. B3LYP method and 6-311++G(d,p) basis set were employed for optimization of molecular structure and calculation of wave numbers of normal modes of vibrations and various other important parameters. Calculated bond lengths and angles were compared with the experimental bond lengths and Bond Angle Parameters. Optimized bond parameters and experimental bond parameters were found in good agreement. Complete potential energy distribution assignments were done successfully by VEDA. The HOMO/LUMO energy gap emphasizes adequate charge transfer happening within the molecule. A study of donor-acceptor interconnections was done via NBO analysis. MEP surface analysis was done to demonstrate charge distribution and reactive areas qualitatively in the molecule. The degree of relative localization of electrons was analyzed via ELF Diagram. The Fukui function analysis showed possible sites for attacks by different substituents. By using the TD-DFT method and PCM solvent model, the UV-Vis spectrum (gas, methanol, DMSO) and the maximum absorption wavelength was computed and compared with experimental data. 3D and 2D intermolecular interactions in the crystal were analyzed via Hirshfeld surface analysis and fingerprint plots reveal that the EAMC crystal was stabilized by H--H/H--H/C--H bond formation. The molecular docking was done with 7 different protein receptors on the molecule to find the best ligand-protein interactions. Molecular dynamic simulations and MMGBSA calculations were also carried out to find out the best binding of the ligand with the protein.Communicated by Ramaswamy H. Sarma.

Structural, vibrational spectroscopy, molecular docking, DFT studies and antibacterial activity of (E)-N1-(3-chlorobenzylidene)benzene-1,4-diamine

In this work, (E)-N1-(3-chlorobenzylidene)benzene-1,4-diamine (CBD) compound was synthesized with good yield. The spectral studies were recorded by FT-IR, FT-Raman, NMR and UV-Vis to determine structural parameters. The geometrical parameters were optimized using DFT calculations at 6-311++G(d,p) basis set. The calculated structural parameters of the molecule were in line with the experimental data. The molecular orbitals of the compound were investigated through highest occupied molecular orbital (HOMO), lowest unoccupied molecular orbital (LUMO) analysis. The hyper conjugative interaction energy E(2) along with donor, acceptor electron densities (EDs) were determined by natural bond orbital (NBO) analysis. The molecular electrostatic potential (MEP), mulliken atomic charges, non-linear optical (NLO) properties and potential energy surface (PES) scan were also calculated. The 1H and 13C NMR chemical shifts calculated using Gauge invariant atomic orbital (GIAO) method were compared with the experimental NMR chemical shifts. Thermogravimetry (TG) and Differential Scanning Calorimetry (DSC) were carried out to characterise the thermal behaviour and stability of CBD molecule. In addition, PreADMET tool was also used to estimate ADME and Toxicity of CBD compound. The compound screened against four pathogens two gram positive and two gram negative had shown good anti-bacterial behaviour. The molecular docking studies executed against anti-bacterial target topoisomerase DNA gyrase enzyme (2XCT) emphasized good binding behaviour over the standard drug.Communicated by Ramaswamy H. Sarma.

Diffuse Basis Functions for Relativistic s and d Block Gaussian Basis Sets

Diffuse s, p, and d functions have been optimized for use with previously reported relativistic basis sets for the s and d blocks of the periodic table. The functions were optimized on the 4:1 weighted average of the s² and p² configurations of the anion, with the d shell in the dⁿ⁺¹ configuration for the d blocks. Exponents were extrapolated for groups 2 and 12, which have unstable or weakly bound anions. The diffuse basis sets have been tested by application to calculations of electron affinities of the group 11 elements (Cu, Ag, and Au), double electron affinities of the group 11 monocations, and potential energy curves of Mg₂ and Ca₂ van der Waals dimers, as well as some response properties of the group 1 anions (Rb^-, Cs^-, and Fr^-), the group 2 elements (Sr, Ba, and Ra), and RbLi, CsLi, and FrLi molecules.

Nonlinear light absorption in many-electron systems excited by an instantaneous electric field: a non-perturbative approach

Applications of low-cost non-perturbative approaches in real time, such as time-dependent density functional theory, for the study of nonlinear optical properties of large and complex systems are gaining increasing popularity. However, their assessment still requires the analysis and understanding of elementary dynamical processes in simple model systems. Motivated by the aim of simulating optical nonlinearities in molecules, here exemplified by the case of the quaterthiophene oligomer, we investigate light absorption in many-electron interacting systems beyond the linear regime by using a single broadband impulse of an electric field; i.e. an electrical impulse in the instantaneous limit. We determine non-pertubatively the absorption cross section from the Fourier transform of the time-dependent induced dipole moment, which can be obtained from the time evolution of the wavefunction. We discuss the dependence of the resulting cross section on the magnitude of the impulse and we highlight the advantages of this method in comparison with perturbation theory by working on a one-dimensional model system for which numerically exact solutions are accessible. Thus, we demonstrate that the considered non-pertubative approach provides us with an effective tool for investigating fluence-dependent nonlinear optical excitations.

Efficient implementation of isotropic cubic response functions for two-photon absorption cross sections within the self-consistent field approximation

Within the self-consistent field approximation, computationally tractable expressions for the isotropic second-order hyperpolarizability have been derived and implemented for the calculation of two-photon absorption cross sections. The novel tensor average formulation presented in this work allows for the evaluation of isotropic damped cubic response functions using only ∼3.3% (one-photon off-resonance regions) and ∼10% (one-photon resonance regions) of the number of auxiliary Fock matrices required when explicitly calculating all the needed individual tensor components. Numerical examples of the two-photon absorption cross section in the one-photon off-resonance and resonance regions are provided for alanine-tryptophan and 2,5-dibromo-1,4-bis(2-(4-diphenylaminophenyl)vinyl)-benzene. Furthermore, a benchmark set of 22 additional small- and medium-sized organic molecules is considered. In all these calculations, a quantitative assessment is made of the reduced and approximate forms of the cubic response function in the one-photon off-resonance regions and results demonstrate a relative error of less than ∼5% when using the reduced expression as compared to the full form of the isotropic cubic response function.

The DIRAC code for relativistic molecular calculations

DIRAC is a freely distributed general-purpose program system for one-, two-, and four-component relativistic molecular calculations at the level of Hartree-Fock, Kohn-Sham (including range-separated theory), multiconfigurational self-consistent-field, multireference configuration interaction, electron propagator, and various flavors of coupled cluster theory. At the self-consistent-field level, a highly original scheme, based on quaternion algebra, is implemented for the treatment of both spatial and time reversal symmetry. DIRAC features a very general module for the calculation of molecular properties that to a large extent may be defined by the user and further analyzed through a powerful visualization module. It allows for the inclusion of environmental effects through three different classes of increasingly sophisticated embedding approaches: the implicit solvation polarizable continuum model, the explicit polarizable embedding model, and the frozen density embedding model.

Quantum mechanical, spectroscopic and docking studies of 2-Amino-3-bromo-5-nitropyridine by Density Functional Method

Experimental and theoretical investigations on the molecular structure, electronic and vibrational characteristics of 2-Amino-3-bromo-5-nitropyridine are presented. The vibrational frequencies were obtained by DFT/B3LYP calculations employing 6-311++G (d, p) basis set. This was compared with experimental FT-IR and FT-Raman spectral data. Simulated FT-IR (4000-400cm^-1) and FT-Raman spectra (4000-100cm^-1) showed good agreement with the observed spectra. The molecular equilibrium geometry of the title compound was fully optimized. Quantum chemical calculations of the equilibrium geometry and the complete vibrational assignments of wavenumbers using potential energy distribution (PED) were calculated with scaled quantum mechanics. HOMO-LUMO energies, energy gap (ΔE), electronegativity (χ), chemical potential (μ), global hardness (η), softness (S) and the Fukui function were calculated for the title molecule. The title compound has a low softness value (0.239) and the calculated value of electrophilicity index (5.905) describes the biological activity. The stability and charge delocalization of the title molecule were studied by Natural Bond Orbital (NBO) analysis, Non-Linear Optical (NLO) behaviour in terms of first order hyperpolarizability, dipole moment and anisotropy of polarizability and Molecular Electrostatic Potential (MEP) were accounted. The computed values of μ, α and β for the title molecule are 1.851 Debye, 1.723×10^-23esu and 7.428×10^-30esu respectively. The high β value and non-zero value of μ indicate that the title compound might be a good candidate for NLO material. Thermodynamic properties of the title molecule were studied for different temperatures thereby revealing the correlations between heat capacity (C), entropy (S) and enthalpy changes (H) with temperatures. Docking studies of the title compound were scrutinized to predict the preferred binding orientation, affinity and activity of the given compound. The title compound was docked into the active site of the protein 5FCT which belongs to the class of proteins exhibiting the property as a Dihydrofolate synthase inhibitor. A minimum binding energy of -5.9kcal/mol and intermolecular energy of -6.5kcal/mol is seen in the interaction.

A combined experimental and theoretical studies on FT-IR, FT-Raman and UV-vis spectra of 2-chloro-3-quinolinecarboxaldehyde

In the present study, the FT-IR and FT-Raman spectra of 2-chloro-3-quinolinecarboxaldehyde (2Cl3QC) have been recorded in the region 4000-400 and 3500-50 cm(-1), respectively. The fundamental modes of vibrational frequencies of 2Cl3QC are assigned. Theoretical information on the optimized geometry, harmonic vibrational frequencies, infrared and Raman intensities were obtained by means of density functional theory (DFT) gradient calculations with complete relaxation in the potential energy surface using 6-31G(d,p) basis set. The vibrational frequencies which were determined experimentally from the spectral data are compared with those obtained theoretically from DFT calculations. A close agreement was achieved between the observed and calculated frequencies by refinement of the scale factors. The infrared and Raman spectra were also predicted from the calculated intensities. Thermodynamic properties like entropy, heat capacity, zero point energy, have been calculated for the molecule. The predicted first hyperpolarizability also shows that the molecule might have a reasonably good non-linear optical (NLO) behavior. The calculated HOMO-LUMO energy gap reveals that charge transfer occurs within the molecule. Stability of the molecule arising from hyper conjugative interactions, charge delocalization have been analyzed using natural bond orbitals (NBO) analysis. The results show that charge in electron density (ED) in the π(∗) antibonding orbitals and E((2)) energies confirms the occurrence of ICT (intra-molecular charge transfer) within the molecule. UV-visible spectrum of the title molecule has also been calculated using TD-DFT/CAM-B3LYP/6-31G(d,p) method. The calculated energy and oscillator strength almost exactly reproduces reported experimental data.

Molecular structure analysis and spectroscopic characterization of 9-methoxy-2H-furo[3,2-g]chromen-2-one with experimental (FT-IR and FT-Raman) techniques and quantum chemical calculations

Experimental and theoretical investigations on the molecular structure, electronic and vibrational characteristics of 9-methoxy-2H-furo[3,2-g]chromen-2-one (9M2HFC) were presented. The vibrational frequencies were obtained by DFT/B3LYP calculations employing 6-311++G(d,p) basis set, were compared with experimental FT-IR and FT-Raman spectral data. The FT-IR spectrum (4000-400 cm(-1)) and FT-Raman spectrum (4000-100 cm(-1)) in solid phase were recorded for 9M2HFC. The geometry of the title compound was fully optimized. Quantum chemical calculations of the equilibrium geometry, the complete vibrational assignments of wavenumbers using potential energy distribution (PED) calculated with scaled quantum mechanics infrared intensities, Raman activities of the title molecule was reported. HOMO-LUMO energies, molecular electrostatic potential, Mulliken population analysis on atomic charges, natural bond orbital (NBO) analysis, non linear optical behavior in terms of first order hyperpolarizability, and thermodynamic properties of the title molecule were carried out. Finally, simulated FT-IR and FT-Raman spectra showed good agreement with the observed spectra.

Vibrational spectra, molecular structure, natural bond orbital, first order hyperpolarizability, thermodynamic analysis and normal coordinate analysis of Salicylaldehyde p-methylphenylthiosemicarbazone by density functional method

The thiosemicarbazone compound, Salicylaldehyde p-methylphenylthiosemicarbazone (abbreviated as SMPTSC) was synthesized and characterized by FTIR, FT-Raman and UV. Density functional (DFT) calculations have been carried out for the title compound by performing DFT level of theory using B3LYP/6-31++G(d,p) basis set. The molecular geometry and vibrational frequencies were calculated and compared with the experimental data. The detailed interpretation of the vibrational spectra has been carried out with aid of normal coordinate analysis (NCA) following the scaled quantum mechanical force field methodology. The electronic dipole moment (μD) and the first hyperpolarizability (βtot) values of the investigated molecule were computed using density functional theory (DFT/B3LYP) with 6-311++G(d,p) basis set. The stability and charge delocalization of the molecule was studied by natural bond orbital (NBO) analysis. Thearomaticities of the phenyl rings were studied using the standard harmonic oscillator model of aromaticity (HOMA) index. Mulliken population analysis on atomic charges is also calculated. The molecule orbital contributions are studied by density of energy states (DOSs).

Normal coordinate analysis and vibrational spectroscopy (FT-IR and FT-Raman) studies of 5-methyl-N-[4-(trifluoromethyl) phenyl]-isoxazole-4-carboxamide using density functional method

Vibrational spectral analysis of 5-methyl-N-[4-(trifluoromethyl) phenyl]-isoxazole-4-carboxamide is (5MN4TPI4C) molecule was carried out using FT-IR and FT-Raman spectroscopic techniques. The equilibrium geometry, harmonic vibrational wavenumbers, various bonding features have been computed using density functional B3LYP method with 6-311G(d,p) as basis set. The assignments of the vibrational spectra have been carried out with the aid of normal coordinate analysis (NCA) following the scaled quantum mechanical force field methodology (SQMFFM). Stability of the molecule arising from hyper conjugative interactions, charge delocalization has been analyzed using natural bond orbital (NBO) analysis. The non-linear optical (NLO) behavior of 5MN4TPI4C has been studied by determination of the electric dipole moment (μ) and hyperpolarizability (β) by using B3LYP/6-311G(d,p) method. The molecular orbital compositions and their contributions to the chemical bonding are studied by Total density of energy states (TDOS), sum of α and β electron (αβDOS) density of states. Thermodynamic properties (heat capacity, entropy and enthalpy) of the title compound at different temperatures are calculated.

FTIR, FT-RAMAN, NMR, spectra, normal co-ordinate analysis, NBO, NLO and DFT calculation of N,N-diethyl-4-methylpiperazine-1-carboxamide molecule

The Fourier Transform Infrared (FT-IR) and FT-Raman of N,N-diethyl-4-methylpiperazine-1-carboxamide (NND4MC) have been recorded and analyzed. The structure of the compound was optimized and the structural characteristics were determined by density functional theory (DFT) using B3LYP method with 6-31G(d,p) and 6-311G(d,p) basis sets. The difference between the observed and scaled wavenumber values of most of the fundamentals is very small. The theoretically predicted FT-IR and FT-Raman spectra of the title molecule have been constructed. The detailed interpretation of the vibrational spectra has been carried out with aid of normal coordinate analysis (NCA) following the scaled quantum mechanical force field methodology. Stability of the molecule arising from hyperconjugative interactions and charge delocalization has been analyzed using natural bond orbital (NBO) analysis. The results show that electron density (ED) in the σ(*) and π(*) antibonding orbitals and second order delocalization energies (E2) confirm the occurrence of intramolecular charge transfer (ICT) within the molecule. The electronic dipole moment (μD) and the first hyperpolarizability (βtot) values of the investigated molecule were computed using Density Functional Theory (DFT/B3LYP) with 6-31G(d,p) and 6-311G(d,p) basis sets. The calculated results also show that the NND4MC molecule may have microscopy nonlinear optical (NLO) behavior with non zero values. Mulliken atomic charges of NND4MC were calculated. The (13)C nuclear magnetic resonance (NMR) chemical shifts of the molecule were calculated by the gauge independent atomic orbital (GIAO) method and compared with experimental results. The UV-Vis spectrum of the compound was recorded. The theoretical electronic absorption spectra have been calculated by using CIS, TD-DFT methods. A study on the electronic properties, such as HOMO and LUMO energies, molecular electrostatic potential (MEP) were also performed.

Synthesis, structural, spectroscopic studies, NBO analysis, NLO and HOMO-LUMO of 4-methyl-N-(3-nitrophenyl)benzene sulfonamide with experimental and theoretical approaches

The sulfonamide compound, 4-methyl-N-(3-nitrophenyl)benzene sulfonamide (abbreviated as 4M3NPBS) has been synthesized and characterized by FTIR, FT-Raman NMR, single crystal X-ray diffraction and thermal analysis. Density functional (DFT) calculations have been carried out for the title compound by performing DFT level of theory using B3LYP/6-31G(d,p) basis set. The calculated results show that the predicted geometry can well reproduce the structural parameters. Predicted vibrational frequencies have been assigned and compared with the experimental IR spectra and they support each other. The UV-Vis spectrum of the compound was recorded. The theoretical electronic absorption spectra have been calculated by using CIS, TD-DFT, ZINDO methods. Stability of the molecule arising from hyper conjugative interactions, charge delocalization has been analyzed using natural bond orbital (NBO) analysis. The first order hyperpolarizability (β0) of 4M3NPBS was calculated using B3LYP/6-31G(d,p) method on the basis of finite-field approach. In addition, frontier molecular orbitals and molecular electrostatic potential were carried out by using density functional theory (DFT/B3LYP) 6-31G(d,p) basis set. The calculated HOMO and LUMO energies show that charge transfer occur in the molecule. The thermal stability of the title compound was determined with the aid of thermo gravimetric analysis (TGA) and differential thermal analysis (DTA).

Vibrational spectra, molecular structure, natural bond orbital, first order hyperpolarizability, TD-DFT and thermodynamic analysis of 4-amino-3-hydroxy-1-naphthalenesulfonic acid by DFT approach

Vibrational spectral analysis of 4-amino-3-hydroxy-1-naphthalenesulfonicacid (4A3HNSA) molecule were carried out using FT-IR and FT-Raman spectroscopic techniques. The equilibrium geometry, harmonic vibrational wavenumbers, various bonding features have been computed using density functional B3LYP method with 6-31G(d,p) as basis set. The Non-Linear Optical (NLO) behavior of 4A3HNSA has been studied by determination of the electric dipole moment (μ) and hyperpolarizability β using HF/6-31G(d,p) method. Stability of the molecule arising from hyperconjugative interactions, charge delocalization have been analyzed using natural bond orbital (NBO) analysis. The results show that charge in electron density (ED) in σ(*) and π(*) antibonding orbitals and second order delocalization energies (E2) confirms the occurrence of Intramolecular Charge Transfer (ICT) within the molecule. The energy and oscillator strength calculated by Time-Dependent Density Functional Theory (TD-DFT) complements with the experimental findings. The simulated spectra satisfactorily coincide with the experimental spectra.

Experimental (FT-IR, FT-Raman and UV-Vis) spectra and theoretical DFT investigations of 2,3-diaminophenazine

Vibrational analysis of the planar electron-rich heterocyclic 2,3-diaminophenazine (DAP) molecule was carried out using FT-IR and FT-Raman spectroscopic techniques. The equilibrium geometry, harmonic vibrational wavenumbers, various bonding features have been computed using density functional method. The calculated molecular geometry parameters have been compared with XRD data. The detailed interpretation of the vibrational spectra has been carried out. The first order hyperpolarizability (β(0)) of the molecular system and related properties (β, α(0) and Δα) of DAP are calculated using HF/6-31G(d,p) method on the basis of finite-field approach. Vibrational analysis reveals that the simultaneous IR and Raman activation of the C-C stretching mode in the molecule provide the evidence for the charge transfer interaction takes place from electron donating group to the ring. The energy and oscillator strength calculated by time-dependent density functional theory (TD-DFT) results complements with the experimental findings. The simulated spectra satisfactorily coincide with the experimental spectra.

Structural, vibrational (FT-IR and FT-Raman) and UV-Vis spectral analysis of 1-phenyl-3-(1,2,3-thiadiazol-5-yl)urea by DFT method

The Fourier Transform Infrared (FT-IR) and FT-Raman spectra of 1-phenyl-3-(1,2,3-thiadiazol-5-yl)urea (TDZ) have been recorded and analyzed. The molecular structure and vibrational spectra (harmonic and anharmonic) of TDZ were calculated by the Density Functional Theory (DFT) method using the B3LYP function with 6-31G(d,p) as the basis set. A detailed interpretation of the Infrared and Raman spectra of TDZ was reported based on Potential Energy Distribution (PED). The first order hyperpolarizability (β(0)) of this novel molecular system and related properties (β, α(0) and Δα) were calculated using HF/6-31G(d,p) method on the finite-field approach. Stability of the molecule arising from hyperconjugative interactions and charge delocalization have been analyzed using Natural Bond Orbital analysis. The results show that charge in electron density in the σ(∗) and π(∗) antibonding orbitals and second order delocalization energies E(2) confirm the occurrence of Intramolecular Charge Transfer (ICT) within the molecule. UV-Vis spectrum of the title molecule was recorded and the electronic properties such as HOMO and LUMO energies were determined by Time-Dependent Density Functional Theory approach.