Gaussian operator bases for correlated fermions

We formulate a general multi-mode Gaussian operator basis for fermions, to enable a positive phase-space representation of correlated Fermi states. The Gaussian basis extends existing bosonic phase-space methods to Fermi systems and thus allows first-principles dynamical or equilibrium calculations in quantum many-body Fermi systems. We prove the completeness of the basis and derive differential forms for products with one- and two-body operators. Because the basis satisfies fermionic superselection rules, the resulting phase space involves only c-numbers, without requiring anticommuting Grassmann variables. Furthermore, because of the overcompleteness of the basis, the phase-space distribution can always be chosen positive. This has important consequences for the sign problem in fermion physics.

A comparative study of the energetics, structures, and mechanisms of the HCN H HNC and LiCN <-> LiNC isomerizations

The potential energy surfaces of the HCN<->HNC and LiCN<->LiNC isomerization processes were determined by ab initio theory using fully optimized triple-zeta double polarization types of basis sets. Both the MP2 corrections and the QCISD level of calculations were performed to correct for the electron correlation. Results show that electron correlation has a considerable influence on the energetics and structures. Analysis of the intramolecular bond rearrangement processes reveals that, in both cases, H (or Li+) migrates in an almost elliptic path in the plane of the molecule. In HCN<->HNC, the migrating hydrogen interacts with the in-plane pi,pi* orbitals of CN, leading to a decrease in the C-N bond order. In LiCN<->LiNC, Li+ does not interact with the corresponding pi,pi* orbitals of CN.

Structural Damage Detection Using Modal Curvature and Fuzzy Logic

A fuzzy logic system (FLS) with a new sliding window defuzzifier is proposed for structural damage detection using modal curvatures. Changes in the modal curvatures due to damage are fuzzified using Gaussian fuzzy sets and mapped to damage location and size using the FLS. The first four modal vectors obtained from finite element simulations of a cantilever beam are used for identifying the location and size of damage. Parametric studies show that modal curvatures can be used to accurately locate the damage; however, quantifying the size of damage is difficult. Tests with noisy simulated data show that the method detects damage very accurately at different noise levels and when some modal data are missing.

Magnetically induced currents and magnetic response properties of molecules

The magnetically induced currents in organic monoring and multiring molecules, in Möbius shaped molecules and in inorganic all-metal molecules have been investigated by means of the Gauge-including magnetically induced currents (GIMIC) method. With the GIMIC method, the ring-current strengths and the ring-current density distributions can be calculated. For open-shell molecules, also the spin current can be obtained. The ring-current pathways and ring-current strengths can be used to understand the magnetic resonance properties of the molecules, to indirectly identify the effect of non-bonded interactions on NMR chemical shifts, to design new molecules with tailored properties and to discuss molecular aromaticity. In the thesis, the magnetic criterion for aromaticity has been adopted. According to this, a molecule which has a net diatropic ring current might be aromatic. Similarly, a molecule which has a net paratropic current might be antiaromatic. If the net current is zero, the molecule is nonaromatic. The electronic structure of the investigated molecules has been resolved by quantum chemical methods. The magnetically induced currents have been calculated with the GIMIC method at the density-functional theory (DFT) level, as well as at the self-consistent field Hartree-Fock (SCF-HF), at the Møller-Plesset perturbation theory of the second order (MP2) and at the coupled-cluster singles and doubles (CCSD) levels of theory. For closed-shell molecules, accurate ring-current strengths can be obtained with a reasonable computational cost at the DFT level and with rather small basis sets. For open-shell molecules, it is shown that correlated methods such as MP2 and CCSD might be needed to obtain reliable charge and spin currents. The basis set convergence has to be checked for open-shell molecules by performing calculations with large enough basis sets. The results discussed in the thesis have been published in eight papers. In addition, some previously unpublished results on the ring currents in the endohedral fullerene Sc3C2@C80 and in coronene are presented. It is shown that dynamical effects should be taken into account when modelling magnetic resonance parameters of endohedral metallofullerenes such as Sc3C2@C80. The ring-current strengths in a series of nano-sized hydrocarbon rings are related to static polarizabilities and to H-1 nuclear magnetic resonance (NMR) shieldings. In a case study on the possible aromaticity of a Möbius-shaped [16]annulene we found that, according to the magnetic criterion, the molecule is nonaromatic. The applicability of the GIMIC method to assign the aromatic character of molecules was confirmed in a study on the ring currents in simple monocylic aromatic, homoaromatic, antiaromatic, and nonaromatic hydrocarbons. Case studies on nanorings, hexaphyrins and [n]cycloparaphenylenes show that explicit calculations are needed to unravel the ring-current delocalization pathways in complex multiring molecules. The open-shell implementation of GIMIC was applied in studies on the charge currents and the spin currents in single-ring and bi-ring molecules with open shells. The aromaticity predictions that are made based on the GIMIC results are compared to other aromaticity criteria such as H-1 NMR shieldings and shifts, electric polarizabilities, bond-length alternation, as well as to predictions provided by the traditional Hückel (4n+2) rule and its more recent extensions that account for Möbius twisted molecules and for molecules with open shells.

Ab initio molecular orbital calculations on ion-molecule and ion pair-molecule complexes of the water-lithium cyanide system

Ab initio molecular orbital (MO) calculations with the 3-21G and 6-31G basis sets were performed on a series of ion-molecule and ion pair-molecule complexes for the H2O + LiCN system. Stabilisation energies (with counter-poise corrections), geometrical parameters, internal force constants and harmonic vibrational frequencies were evaluated for 16 structures of interest. Although the interaction energies are smaller, the geometries and relative stabilities of the monohydrated contact ion pair are reminiscent of those computed for the complexes of the individual ions. Thus, interaction of the oxygen lone pair with lithium leads to a highly stabilised C2v structure, while the coordination of water to the cyanide ion involves a slightly non-linear hydrogen bond. Symmetrical bifurcated structures are computed to be saddle points on the potential energy surface, and to have an imaginary frequency for the rocking mode of the water molecule. On optimisation the geometries of the solvent shared ion pair structures (e.g. Li+cdots, three dots, centered OH2cdots, three dots, centered CN−) revealed a proton transfer from the water molecule leading to hydrogen bonded forms such as Li-O-Hcdots, three dots, centered HCN. The variation in the force constants and harmonic frequencies in the various structures considered are discussed in terms of ion-molecular and ion pair-molecule interactions.

AB-initio study of alpha-silyl and beta-silyl substituent effects on vinyl radicals

The geometries of alpha- and beta-silyl substituted vinyl radicals and of alpha,beta-disilylvinyl radical have been optimised with the STO-3G and the STO-3G* basis sets. The relative stabilities of various conformers have been determined at the UMP2/6-31G* level. The stabilisation of vinyl radicals through alpha-silyl substitution is larger than that due to corresponding alkyl groups. The presence of an alpha-silyl group also leads to a tendency towards linearisation of the vinyl radical centre and a corresponding reduction in the inversion barrier. In marked contrast, the beta-silyl effect is negligible. The geometric, conformational and energetic consequences are insignificant. Overall, the silyl substituent effect at vinyl radicals is very different from that computed earlier for the vinyl cations, but qualitatively similar to that found in carbanions.

SCF and electron correlation studies on structure, force constants and vibrational spectra of borane monoammoniate complex

Effects of basis set and electron correlation on the equilibrium geometry, force constants and vibrational spectra of BH3NH3 have been studied. A series of basis sets ranging from double zeta to triple zeta including polarization and diffuse functions have been utilized. All the SCF based calculations overestimate the dative B-N bond distance and considerable improvement occurs when the treatment for electron correlation is introduced. Detailed vibrational analysis for BH3NH3 has been carried out. The mean absolute percentage deviation of the ab initio predicted vibration frequencies of (BH3NH3)-B-11 from the experiment is about 10% for the SCF based calculations and the MP2 method shows better agreement, the overall deviation being 5-6%. The ground state effective force constants of BH3NH3 were obtained using RECOVES procedure. The RECOVES sets of force constants are found to be highly satisfactory for the prediction of the vibrational frequencies of different isotopomers of BH3NH3. The mean absolute percentage deviation of the calculated frequencies of different isotopomers from the experiment is much less than 1%. The RECOVES-MP2/augDZP set of force constants was found to be the best set among the different sets for this molecule. Theoretical infrared intensities are in fair agreement with the observed spectral features.

[n]peristylanes and [n]oxa[n]peristylanes (n = 3-6): A theoretical study

Theoretical studies at the HF and Becke3LYP levels using 6-31G* basis sets were carried out on a series of [n]peristylanes and [n]oxa[n]peristylanes (n = 3-6) to understand their structure and energetics. The structures of the [3]- and [4]peristylanes (1, 2) and their era-derivatives (5, 6) were calculated to have the anticipated high symmetry, C-nv. In contrast, a C-s structure (9) at HF/6-31G* and another (25) at the Becke3LYP/6-31G* level were calculated for the [5]oxa[5]peristylane. The energy difference between them is extremely small even though there are major differences in the structures indicating every soft potential energy surface: On the other hand, the potential energy surface of [6]oxa[6]peristylane is not as soft. Similar structures were also calculated for the top rings. Calculations on the seco-compounds 11-14 and 15-19 (Table 4) indicate that there is no unusual strain involved in the formation of 27 from 19. The Li+ interaction energies of the [n]oxa[n]peristylanes are 61.7 (n = 3), 72.8 (n = 4), 84.2 (n = 5) and 91.7 (n = 6) kcal mol(-1) at the Becke3LYP/6-3IG* level. Dramatic differences between the C-C bond lengths obtained from the solid state X-ray diffraction studies and those from the calculations for the [n]oxa[n]peristylanes were also observed.

Ab initio studies of solvent effect on conformational equilibria and vibrational spectra of dipropionamide

Systematic ab initio molecular orbital studies of the conformational equilibria and vibrational spectra of dipropionamide using the basis sets 6-31g(d) and 6-31++G(d,p) have been carried out. The vibrational spectra of dipropionamide have been satisfactorily interpreted taking into account the agreement between the calculated frequencies, infrared and Raman band intensities and the shifts in the spectra of deuterated molecules with those observed. The previous assignments of most of the vibrational bands are well confirmed, a few bands need reassignment, however. The solvent effects were investigated by self-consistent reaction field theory using dipole and self-consistent isodensity polarized continuum model methods. The introduction of a dielectric medium has only a marginal effect on the conformational equilibria and vibrational spectra. However, the calculated changes in geometry and vibrational spectra on going from the gas phase to the solution phase are in accord with the increasing weight of the dipolar resonance structure in polar solvents. (C) 2002 Elsevier Science B.V. All rights reserved.

Ab initio, DFT and transition state theory calculations on 1,2-HF, HCI and CIF elimination reactions from CH2F-CH2Cl

This paper reports ab intio, DFT and transition state theory (TST) calculations on HF, HCI and CIF elimination reactions from CH2Cl-CH2F molecule. Both the ground state and the transition state for HX elimination reactions have been optimized at HF, MP2 and DFT calculations with 6-31G*, 6-31G** and 6-311++G** basis sets. In addition, CCSD(T) single point calculations were carried out with MP2/6-311++G** optimized geometry for more accurate determination of the energies of the minima and transition state, compared to the other methods employed here. Classical barriers are converted to Arrhenius activation energy by TST calculations for comparisons with experimental results. The pre-exponential factors, A, calculated at all levels of theory are significantly larger than the experimental values. For activation energy, E-a DFT gives good results for HF elimination, within 4-8 W mol(-1) from experimental values. None of the methods employed, including CCSD(T), give comparable results for HCI elimination reactions. However, rate constants calculated by CCSD(T) method are in very good agreement with experiment for HCI elimination and they are in reasonable agreement for HF elimination reactions. Due to the strong correlation between A and E., the rate constants could be fit to a lower A and E-a (as given by experimental fitting, corresponding to a tight TS) or to larger A and E-a (as given by high level ab initio calculations, corresponding to a loose TS). The barrier for CIF elimination is determined to be 607 U mol(-1) at HF level and it is unlikely to be important for CH2FCH2Cl. Results for other CH2X-CH2Y (X,Y = F/Cl) are included for comparison.

Thermal decomposition of haloethanols: single pulse shock tube and ab initio studies

This paper reports single pulse shock tube and ab initio studies on thermal decomposition of 2-fluoro and 2-chloroethanol at T=1000–1200 K. Both molecules have HX (X = F/Cl) and H2O molecular elimination channels. The CH3CHO formed by HX elimination is chemically active and undergoes secondary decomposition resulting in the formation of CH4, C2H6, and C2H4. A detailed kinetic simulation indicates that the formation of C2H4 could not be quantitatively explained as arising exclusively from secondary CH3CHO decomposition. Contributions from primary radical processes need to be considered to explain C2H4 quantitatively. Ab initio calculations on HX and H2O elimination reactions from the haloethanols at HF, MP2, and DFT levels with various basis sets up to 6/311++G**are reported. It is pointed out that due to strong correlations between A and Eα, comparison of these two parameters between experimental and theoretical results could be misleading.

Study on the structure and vibrational spectra of efavirenz conformers using DFT: Comparison to experimental data

Efavirenz, (S)-6-chloro-4-(cyclopropylethynyl)-1,4-dihydro-4-(trifluoromethyl)-2H-3 ,1-benzoxazin-2-one, is an anti HIV agent belonging to the class of the non-nucleoside inhibitors of the HIV-1 virus reverse transcriptase. A systematic quantum chemical study of the possible conformations, their relative stabilities and vibrational spectra of efavirenz has been reported. Structural and spectral characteristics of efavirenz have been studied by vibrational spectroscopy and quantum chemical methods. Density functional theory (DFT) calculations for potential energy curve, optimized geometries and vibrational spectra have been carried out using 6-311++G(d,p) basis sets and B3LYP functionals. Based on these results, we have discussed the correlation between the vibrational modes and the crystalline structure of the most stable form of efavirenz. A complete analysis of the experimental infrared and Raman spectra has been reported on the basis of wavenumber of the vibrational bands and potential energy distribution. The infrared and the Raman spectra of the molecule based on OFT calculations show reasonable agreement with the experimental results. The calculated HOMO and LUMO energies shows that charge transfer occur within the molecule. (C) 2011 Elsevier B.V. All rights reserved.

Cu-II-azide polynuclear complexes of Cu-4 building clusters with Schiff-base co-ligands: synthesis, structures, magnetic behavior and DFT studies

Three new copper-azido complexes Cu-4(N-3)(8)(L-1)(2)](n) (1), Cu-4(N-3)(6)(L-2)(2)(H2O)(2)] (2), and Cu-4(N-3)(6)(L-3)(2)](n) (3) L-1 is the imine resulting from the condensation of pyridine-2-carboxaldehyde with N-methylethylenediamine, HL2 and HL3 are the condensation products of 2-hydroxy-3-methoxybenzaldehyde with N,N-diethylethylenediamine and N-ethylethylenediamine respectively] have been synthesized by using 0.5 molar equivalents of the Schiff base ligands with Cu(NO3)(2)center dot 3H(2)O and an excess of NaN3. Single crystal X-ray structures show that the basic unit of these complexes contains very similar Cu-4(II) building blocks. While 1 and 3 have overall 1D structures, 2 forms discrete tetranuclear clusters due to blocking of two coordination sites on the tetranuclear cluster by water molecules. Magnetic susceptibility measurements over a wide range of temperatures exhibit the presence of both antiferromagnetic and ferromagnetic exchanges within the tetranuclear unit structures. Density functional theory calculations (using B3LYP functional and two different basis sets) have been performed on the complexes 1-3 to provide a qualitative theoretical interpretation of their overall magnetic behavior.

Cu-II-Azide Polynuclear Complexes of Three Different Building Clusters with the Same Schiff-Base Ligand: Synthesis, Structures, Magnetic Behavior, and Density Functional Theory Studies

Three copper-azido complexes Cu-4(N-3)(8)(L-1)(2)(MeOH)(2)](n) (1), Cu-4(N-3)(8)(L-1)(2)] (2), and Cu-5(N-3)(10)(L-1)(2)](n) (3) L-1 is the imine resulting from the condensation of pyridine-2-carboxaldehyde with 2-(2-pyridyl)ethylamine] have been synthesized using lower molar equivalents of the Schiff base ligand with Cu(NO3)(2)center dot 3H(2)O and an excess of NaN3. Single crystal X-ray structures show that the basic unit of the complexes 1 and 2 contains Cu-4(II) building blocks; however, they have distinct basic and overall structures due to a small change in the bridging mode of the peripheral pair of copper atoms in the linear tetranudear structures. Interestingly, these changes are the result of changing the solvent system (MeOH/H2O to EtOH/H2O) used for the synthesis, without changing the proportions of the components (metal to ligand ratio 2:1). Using even lower proportions of the ligand, another unique complex was isolated with Cu-5(II) building units, forming a two-dimensional complex (3). Magnetic susceptibility measurements over a wide range of temperature exhibit the presence of both antiferromagnetic (very weak) and ferromagnetic exchanges within the tetranuclear unit structures. Density functional theory calculations (using B3LYP functional, and two different basis sets) have been performed on the complexes 1 and 2 to provide a qualitative theoretical interpretation of their overall magnetic behavior.

Conformational Stability and Intramolecular Hydrogen Bonding in 1,2-Ethanediol and 1,4-Butanediol

The gas-phase infrared spectra of 1,2-ED and 1,4-BD have been, recorded at three different temperatures using a multipass gas cell of 6 m optical path length. DFT calculation has also been carried out using 6-311++G** and aug-cc-pVDZ basis sets to look for the existence of intramolecular hydrogen bonding, in them from the red shift and infrared absorption intensity enhancement of the bonded O-H band compared to that of the free O-H hand. Equilibrium population analysis With 10 conformers of 1,2-ED and 1,4-BD at experimental temperatures were-carried out for the reconstruction of the Observed vibrational spectra at that temperature,using standard statistical relationships. The most abundant conformer at experimental temperatures, was identified. In 1,2-ED a red shift of 45 cm(-1) in the intramolecularly interacting O-H stretching vibrational band position and no significant intensity enhancement compared to that of-the free O-H have been observed. On the contrary, in one of the hydrogen bonded conformers of 1,4-BD, a 124 cm(-1) red shift in the O-H stretching frequency and 8.5 times-intensity enhancement for the ``bonded'' O-H compared to that of the ``free'' O-H is seen. On the basis of this comparative study, we have concluded that strong intramolecular hydrogen bonding exists in 1,4-BD. But there appears, to be weak intramolecular hydrogen bonding in 1,2-ED at temperatures of 303, 313, and 323 Km the gas phase We have found that most stable hydrogen-bonded conformers of 1,4-BD are less populated than some of the non-hydrogen-bonded conformers. Even for the 1,4-BD, the relative population of the g'GG'Gt conformer, which has a strong intramolecular,hydrogen bond, is less than what is predicted. Perhaps the intramolecular hydrogen bond plays a less Significant role in the relative stability,of the various Conformers than what has been predicted from calculations and prevails in the literature.