A finite-nucleus model for relativistic electronic structure calculations using a Douglas-Kroll-transformed Hamiltonian

Summary We report the implementation of a Gaussian finite-nucleus model in the framework of the spin-free no-pair Hamiltonian obtained from the Douglas-Kroll transformation of the no-pair operator with external-field projectors. The finite nucleus regularizes the weak singularity of the wavefunction at the locations of the nuclei and provides a means for efficient exponent optimization using a spinaveraged relativistic one-component operator. We report and discuss basis sets for the gold atom obtained from various optimization procedures, making use of a point nucleus as well as employing various finite-nucleus models.Dedicated to Prof. Dr. Werner Kutzelnigg on occasion of his sixtieth birthday     

Vibrational investigation, molecular orbital studies and molecular electrostatic potential map analysis on 3-chlorobenzoic acid using hybrid computational calculations

The FT-IR and FT-Raman spectra of 3-chlorobenzoic acid (3CBA) are recorded in the liquid state. The fundamental vibrational frequencies, intensity of vibrational bands and the optimized geometrical parameters of the compound are evaluated using HF and DFT (LSDA/B3LYP/B3PW91/MPW1PW91) methods with 6-311+G(d,p) basis set. The theoretical wave numbers are scaled down and compared with the experimental values which showed very good agreement. Comparison of stimulated spectra with the experimental spectra provides important information about the ability of the hybrid computational method to describe the vibrational modes. The HOMO, LUMO, chemical hardness (η), chemical potential (μ), electrophilicity values (ω) and maximum amount of electronic charge transfer (ΔN(max)) are calculated. The molecular electrostatic potential (MESP) is calculated and the corresponding graphs are drawn. Some thermodynamic parameters and physico-chemical properties are calculated and discussed.     

Vector coupling coefficients for calculations of transition-metal atoms and ions by the SCF coupling operator method

We derived the necessary conditions to which the vector coupling coefficients (VCC ) a and b describing atomic L,S-multiplets of the configurations d^N (1 ≤ N ≤ 9), should satisfy. Special attention is paid to the states of non-Roothaan type for which VCC depend on the choice of degenerate d-orbitals basis set determined within the accuracy up to an orthogonal transformation u. It is shown that for such states the direct sum of matrices ‖a‖ and ‖b‖ must be the non-symmetric matrix. Obtained VCC were used for the ab initio calculations (basis set (14s9p5d)/[8s4p2d] from [15]) on first-row transition atoms (from Sc to Cu) to compare to similar calculations [16], in which the Peterson's VCC have been used, and with calculations [15] carried out by the atomic SCF program [4] as well.     

Efficient generation of Heisenberg Hamiltonian matrices for VB calculations of potential energy surfaces

The spin‐Hamiltonian valence bond theory relies upon covalent configurations formed by singly occupied orbitals differing by their spin counterparts. This theory has been proven to be successful in studying potential energy surfaces of the ground and lowest excited states in organic molecules when used as a part of the hybrid molecular mechanics—valence bond method. The method allows one to consider systems with large active spaces formed by n electrons in n orbitals and relies upon a specially proposed graphical unitary group approach. At the same time, the restriction of the equality of the numbers of electrons and orbitals in the active space is too severe: it excludes from the consideration a lot of interesting applications. We can mention here carbocations and systems with heteroatoms. Moreover, the structure of the method makes it difficult to study charge‐transfer excited states because they are formed by ionic configurations. In the present work we tackle these problems by significant extension of the spin‐Hamiltonian approach. We consider (i) more general active space formed by n ± m electrons in n orbitals and (ii) states with the charge transfer. The main problem addressed is the generation of Hamiltonian matrices for these general cases. We propose a scheme combining operators of electron exchange and hopping, generating all nonzero matrix elements step‐by‐step. This scheme provides a very efficient way to generate the Hamiltonians, thus extending the applicability of spin‐Hamiltonian valence bond theory. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Theoretical calculations of hyperfine coupling constants for muoniated butyl radicals

The hyperfine coupling constants (HFCCs) of all the butyl radicals that can be produced by muonium (Mu) addition to butene isomers (1- and 2-butene and isobutene) have been calculated, to compare with the experimental results for the muon and proton HFFCs for these radicals reported in paper II (Fleming, D. G.; et al. J. Phys. Chem. A 2011, 10.1021/jp109676b) that follows. The equilibrium geometries and HFCCs of these muoniated butyl radicals as well as their unsubstituted isotopomers were treated at both the spin-unrestricted MP2/EPR-III and B3LYP/EPR-III levels of theory. Comparisons with calculations carried out for the EPR-II basis set have also been made. All calculations were carried out in vacuo at 0 K only. A C-Mu bond elongation scheme that lengthens the equilibrium C-H bond by a factor of 1.076, on the basis of recent quantum calculations of the muon HFCCs of the ethyl radical, has been exploited to determine the vibrationally corrected muon HFCCs. The sensitivity of the results to small variations around this scale factor was also investigated. The computational methodology employed was "benchmarked" in comparisons with the ethyl radical, both with higher level calculations and with experiment. For the β-HFCCs of interest, compared to B3LYP, the MP2 calculations agree better with higher level theories and with experiment in the case of the eclipsed C-Mu bond and are generally deemed to be more reliable in predicting the equilibrium conformations and muon HFCCs near 0 K, in the absence of environmental effects. In some cases though, the experimental results in paper II demonstrate that environmental effects enhance the muon HFCC in the solid phase, where much better agreement with the experimental muon HFCCs near 0 K is found from B3LYP than from MP2. This seemingly better level of agreement is probably fortuitous, due to error cancellations in the DFT calculations, which appear to mimic these environmental effects. For the staggered proton HFCCs of the butyl radicals exhibiting no environmental effect in paper II, the best agreement with experiment is consistently found from the B3LYP calculations, in agreement also with benchmark calculations carried out for the ethyl radical.     

Multicenter point charge model for high-quality molecular electrostatic potentials from AM1 calculations

The natural atomic orbital/point charge (NAO-PC) model based upon the AM1 wave function has been developed to calculate molecular electrostatic potentials (MEPs). Up to nine point charges (including the core charge) are used to represent heavy atoms. The positions and magnitudes of the eight charges that represent the atomic electron cloud are calculated from the natural atomic orbitals (NAOs) and their occupations. Each hybrid NAO is represented by two point charges situated at the centroid of each lobe. The positions of the centroids and the magnitudes of the charges were obtained by numerical integration of the Slater-type hybrids and the results used to set up polynomials and look-up tables that replace the integration step in the actual MEP calculation. The MEPs calculated using this method are found to be in better agreement with those obtained using RHF/6-31G* than those obtained from the AM1 wave function using Coulson charges or with MOPAC-ESP. The MEP calculations are extremely fast and have, for instance, been incorporated into an interactive graphics package. © 1993 John Wiley & Sons, Inc.     

DelPhiForce,a tool for electrostatic force calculations: Applications to macromolecular binding

Long‐range electrostatic forces play an important role in molecular biology, particularly in macromolecular interactions. However, calculating the electrostatic forces for irregularly shaped molecules immersed in water is a difficult task. Here, we report a new tool, DelPhiForce, which is a tool in the DelPhi package that calculates and visualizes the electrostatic forces in biomolecular systems. In parallel, the DelPhi algorithm for modeling electrostatic potential at user‐defined positions has been enhanced to include triquadratic and tricubic interpolation methods. The tricubic interpolation method has been tested against analytical solutions and it has been demonstrated that the corresponding errors are negligibly small at resolution 4 grids/Å. The DelPhiForce is further applied in the study of forces acting between partners of three protein–protein complexes. It has been demonstrated that electrostatic forces play a dual role by steering binding partners (so that the partners recognize their native interfaces) and exerting an electrostatic torque (if the mutual orientations of the partners are not native‐like). The output of DelPhiForce is in a format that VMD can read and visualize, and provides additional options for analysis of protein–protein binding. DelPhiForce is available for download from the DelPhi webpage at http://compbio.clemson.edu/downloadDir/delphiforce.tar.gz © 2017 Wiley Periodicals, Inc.     

NDDO MO calculations: isotropic hyperfine coupling constants and nuclear spin-spin coupling constants

The nonempirical NDDO MO method in its unrestricted form has been used to evaluate isotropic hyperfine coupling constants and nuclear spin-spin coupling constants. Satisfactory agreement with INDO and experimental results is obtained.     

A hybrid technique of orthonormality constrained orbital optimization in SCF calculations

A recently proposed orthonormality constrained orbital optimization technique is operationally modified further by coupling it to a gradient biased method, namely the steepest descent procedure of McWeeny. The hybrid technique developed in this way is shown to have better convergence properties in closed and unrestricted open-shell calculations. The technique can be adapted to MCSCF procedures as well. The important role played by "orbital symmetries" in the operation of the method is analysed. Similarities and differences of the present method with the orthogonal gradient method are pointed out. Possible avenues of circumventing convergence difficulty that one may encounter in pathological cases, particularly in ab initio calculations involving extended basis set, are suggested.     

Competition of hydrogen-bonding and electrostatic interactions within hybrid polymer multilayers

Using a layer-by-layer sequential adsorption technique, we report the construction of hybrid films in which layers of hydrogen-bonded polymers are embedded within electrostatically associated polyelectrolytes. The components of the hybrid film include a neutral hydrogen-bonding polymer, a weak polycarboxylic acid, and a strong polycation. Depending on the pH value used for the deposition of the electrostatic film, we found two distinctive regimes of film growth. At pHs lower than a critical value, deposition of electrostatic layers occurred on top of hydrogen-bonded stacks to produce hybrid, three-component films. At pHs higher than a critical value, neutral, hydrogen-bonded chains were displaced by the adsorbing chains of the polycation, producing two-component films. The property of the hydrogen-bonded stacks of hybrid films to be selectively dissolved by exposing them to a high pH makes these films promising candidates for producing free polyelectrolyte films.     

Local unitary transformation method for large-scale two-component relativistic calculations: case for a one-electron Dirac Hamiltonian

An accurate and efficient scheme for two-component relativistic calculations at the spin-free infinite-order Douglas-Kroll-Hess (IODKH) level is presented. The present scheme, termed local unitary transformation (LUT), is based on the locality of the relativistic effect. Numerical assessments of the LUT scheme were performed in diatomic molecules such as HX and X(2) (X = F, Cl, Br, I, and At) and hydrogen halide clusters, (HX)(n) (X = F, Cl, Br, and I). Total energies obtained by the LUT method agree well with conventional IODKH results. The computational costs of the LUT method are drastically lower than those of conventional methods since in the former there is linear-scaling with respect to the system size and a small prefactor.     

Poisson-Boltzmann calculations versus molecular dynamics simulations for calculating the electrostatic potential of a solvated peptide

We performed a very long molecular dynamics simulation of a peptide in explicit water molecules and ions and averaged the electrostatic potential caused by peptide, water and ions at eight points in the vicinity of the peptide. These electrostatic potential values were directly compared to the potential calculated by solving the non-linear Poisson-Boltzmann equation for the system, which describes the solvent using continuum electrostatics. We analyze the contribution of dielectric constant, conformational flexibility and solvation effects on the electrostatic potential at these eight points. Received: 24 April 1998 / Accepted: 4 August 1998 / Published online: 23 November 1998     

Inclusion of discrete-to-continuum coupling in multiphoton excitation and dissociation calculations

A nonperturbative approach to multiphoton excitation and dissociation of molecules is presented in which coupling to the continuum is treated explicitly. Transitions among continuum levels are not modeled directly, but something of their effect is represented by assuming that the continuum population density is so low as to be effectively zero at all times. Two trial applications are briefly discussed.     

Stability and adsorption properties of electrostatic complexes: design of hybrid nanostructures for coating applications

We report the presence of a correlation between the bulk and interfacial properties of electrostatic coacervate complexes. Complexes were obtained by co-assembly between cationic-neutral diblocks and oppositely charged surfactant micelles or 7 nm cerium oxide nanoparticles. Light scattering and reflectometry measurements revealed that the hybrid nanoparticle aggregates were more stable through both dilution and rinsing (from either a polystyrene or a silica surface) than their surfactant counterparts. These findings were attributed to a marked difference in critical association concentration between the two systems and to the frozen state of the hybrid structures.     

Ab initio calculations of electrostatic potentials and deformation densities for a series of choline ester model systems

Ab initio LCAO-MO-SCF calculations using a double zeta basis set have been performed for the methyl esters of acetic acid, carbamic acid, methylcarbonic acid, and trifluoracetic acid, in order to model the corresponding choline esters. The systems have been compared by means of population analyses, electron density differences, electrostatic potentials and potential differences. The significance of the electrostatic potential in connection with crystal structure and packing has been studied. The differences in the proton affinity of the compounds have been correlated to differences in the potentials.     

A simplified relativistic time-dependent density-functional theory formalism for the calculations of excitation energies including spin-orbit coupling effect

In the present work we have proposed an approximate time-dependent density-functional theory (TDDFT) formalism to deal with the influence of spin-orbit coupling effect on the excitation energies for closed-shell systems. In this formalism scalar relativistic TDDFT calculations are first performed to determine the lowest single-group excited states and the spin-orbit coupling operator is applied to these single-group excited states to obtain the excitation energies with spin-orbit coupling effects included. The computational effort of the present method is much smaller than that of the two-component TDDFT formalism and this method can be applied to medium-size systems containing heavy elements. The compositions of the double-group excited states in terms of single-group singlet and triplet excited states are obtained automatically from the calculations. The calculated excitation energies based on the present formalism show that this formalism affords reasonable excitation energies for transitions not involving 5p and 6p orbitals. For transitions involving 5p orbitals, one can still obtain acceptable results for excitations with a small truncation error, while the formalism will fail for transitions involving 6p orbitals, especially 6p1/2 spinors.