4-Component relativistic magnetically induced current density using London atomic orbitals

We present the implementation and application of 4-component relativistic magnetically induced current density using London atomic orbitals for self-consistent field models. We obtain a magnetically balanced basis by a simple scheme where orbitals obtained by imposing restricted kinetic balance are extended by their unrestricted kinetic balance complement. The presented methodology makes it possible to analyze the concept of aromaticity based on the ring current criterion for closed-shell molecules across the periodic table and is independent of the choice of gauge origin. As a first illustration of the methodology we study plots of the magnetically induced current density and its divergence in the series C(5)H(5)E (E = CH, N, P, As, Sb, Bi) at the Kohn-Sham level, as well as integrated ring current susceptibilities, which we compare to previous results (R. Bast et al., Chem. Phys., 2009, 356, 187) obtained using a common gauge origin approach. We find that the current strength decreases monotonically along the series, but that all molecules qualify as aromatic according to the ring current criterion.     

Spin–orbit effects on magnetically induced current densities in the () clusters

This study reports the spin–orbit effects on the aromaticity of the , , , , , and anionic clusters via the magnetically induced current‐density method. All‐electron density functional theory (DFT) calculations were carried out using the four‐component Dirac‐Coulomb (DC) hamiltonian, including scalar and spin–orbit relativistic effects. The magnetic index of aromaticity was calculated by numerical integration over the current flow between two atoms in the pentagonal ring. These values were compared to the spin‐free values (spin–orbit coupling switched off), in order to assess the spin–orbit effect on aromaticity. It was found that in the heavy anions, and , there is a significant influence of the spin–orbit coupling. © 2018 Wiley Periodicals, Inc.     

Theoretical investigation of photoelectron spectra and magnetically induced current densities in ring-shaped transition-metal oxides

The molecular structures of cyclic group 6 transition-metal (M = Cr, Mo, W, Sg) oxides (M ₃O ₉ ^0/1?/2? ) species have been optimized at density functional theory (DFT) levels. The photoelectron spectra (PES) of M ₃O ₉ ^? (M = Cr, Mo, W) were calculated at the time-dependent DFT and approximate coupled-cluster singles doubles (CC2) levels and compared with experimental results. The CC2 calculations did not yield any reliable PES, whereas the molecular structures can be identified by comparing PES obtained at the DFT level with experiment. Magnetically induced current densities were calculated at the DFT level using the gauge-including magnetically induced current (gimic) approach. The current strengths and current pathways of the neutral M₃O₉ and the dianionic M₃O ₉ ^2? (M = Cr, Mo) oxides were investigated and analyzed with respect to a previous prediction of d-orbital aromaticity for Mo₃O₉ anions. Current-density calculations provide ring-current strengths that are used to assess the degree of aromaticity. Comparison of current-density calculations and calculations of nucleus-independent chemical shifts (NICS) shows that NICS calculations are not a reliable tool for determining the degree of aromaticity of the metal oxides.     

An investigation of aromaticity in hydroxybenzenes based on the study of magnetically induced current density

Hydroxybenzenes are derivatives of benzene in which one or more H atoms are replaced by OH groups. This work considers the effect of such replacement on the aromaticity of the ring, by assessing the aromaticity through the values and patterns of the magnetically induced current density. The results show that aromaticity is quenched and the quenching depends on the number and positions of the substituting OH groups. The isosurfaces of the magnitude of the current density also highlight connections between the donor and acceptor atoms forming an intramolecular hydrogen bond.     

Hydrogen-bond strengths by magnetically induced currents

We present here a "non-invasive" computational method to estimate the strength of individual hydrogen bonds using magnetically induced currents. The method is calibrated using H-bonding dimers, and applied on Watson-Crick DNA base pairs and proton wires in carbonic anhydrase.     

The detection of mercury vapour by magnetically induced optical rotation

Equations to describe magnetically-induced optical rotation signals obtained from mercury atoms in the vapour phase are described and tested. The method gives good analytical sensitivity and selectivity. Its advantages over the conventional techniques of analytical spectrometry include the automatic tolerance of continuum background absorption and of source drift.     

Different mutual positions of double bonds in open carbon chains and corresponding information from magnetically induced current densities

Magnetically induced current densities have been extensively utilized to investigate aromaticity. However, an applied external magnetic field induces a current in all chemical bonds. The present work investigates the values and patterns of magnetically induced current densities in open hydrocarbon chains containing two or more CC double bonds, considering the three possibilities for the mutual positions of the double bonds (conjugated, isolated, and consecutive) and comparing the corresponding current patterns. The patterns for chains containing three conjugated double bonds are also compared with the patterns in the aromatic benzene ring. The study highlights different current patterns for the three categories of double-bonds arrangements. The differences can be interpreted in terms of presence or absence of electron delocalization (present only in the case of conjugated double bonds) and, for consecutive double bonds, also in terms of the mutual perpendicularity of the π orbitals of adjacent double bonds.     

Observation of magnetically induced polarization in a ferrofluid

Light transmitted through a magnetite colloid subjected to a pulsed magnetic field is analyzed to determine particle size. Absorption measurements of light give a measure of the magnetic dipole moment. A comparison with theoretical estimates for particle coagulation is made and implications for medical research are noted.     

The magnetically induced circular emission spectrum of Cr(III) in guanidinium aluminium sulphate

The magnetically induced, circular emission spectrum (MICE) of Cr(III) in hexagonal guanidinium aluminium sulphate hexahydrate has been measured together with the axial emission spectrum. Four lines are observed in the ²E region suggesting that these are the result of trigonal field and site symmetry splittings since each component exhibits an A term type dispersion. The MICE spectrum enables several of the vibrational features to be identified.     

An approximate MO-LCAO-SCF method including overlap

A molecular orbital method is described which can be applied to molecules where the restrictions of π-electron theory are not fulfilled. It has the following main characteristics: 1. Atomic SCF functions are used as a basis, 2. Core-valence interactions are treated by means of perturbation theory, 3. Mulliken type approximations are used for many-center integrals.     

The direct current method for measuring charged membrane conductance

This paper deals with a method for measuring electrical resistance in charged membranes. The method is based on the application of a step change in direct current and on the analysis of the potential transient subsequent to the application of the current step. p]Membrane electrical resistance was determined by an extrapolation to zero time of potential differences measured after the current step. p]Experimental results obtained with commercial ion-exchange membranes were in good agreement with those computed from the Fick equation. The method developed gives more accurate values with a standard deviation lower than traditional techniques and allows the resistance of an asymmetrical membrane to be determined in both current directions.     

The contributions of magnetically induced dichroism and the resonant Voigt effect to the detection of silver by atomic spectroscopy

A simple theoretical analysis, based on the classical theory of the transmission of polarized light through a polarizing medium, has been developed to examine the factors determining the intensity of resonance radiation transmitted by a magnetized atomic vapour interposed between a pair of crossed polarizers. In common with other studies, this analysis predicts a quadratic dependence of the transmitted intensity on atom concentration but in addition, indicates that even at low atom concentrations differential absorption (dichroism) may make a significant addition to the transmitted intensity arising from differential refraction (birefringence). Further, by introducing a small angular off-set between the crossed polarizers, the concentration dependence of the transmitted intensity at low concentration may be increased and linearized. Experiments using a flame-generated atomic vapour of silver with a silver hollow cathode lamp or deuterium lamp as light source confirmed these predictions. When a strong magnetic field (9.75 kG) was used, a detection limit of 0.3 ppm (= μg/ml) for silver, with a near linear response up to 20 ppm was obtained.     

Implementation of divide-and-conquer method including Hartree-Fock exchange interaction

The divide-and-conquer (DC) method, which is one of the linear-scaling methods avoiding explicit diagonalization of the Fock matrix, has been applied mainly to pure density functional theory (DFT) or semiempirical molecular orbital calculations so far. The present study applies the DC method to such calculations including the Hartree-Fock (HF) exchange terms as the HF and hybrid HF/DFT. Reliability of the DC-HF and DC-hybrid HF/DFT is found to be strongly dependent on the cut-off radius, which defines the localization region in the DC formalism. This dependence on the cut-off radius is assessed from various points of view: that is, total energy, energy components, local energies, and density of states. Additionally, to accelerate the self-consistent field convergence in DC calculations, a new convergence technique is proposed.     

The role of induced current density in Steroelectronic effects: Perlin effect

The normal and reverse Perlin effect is usually explained by the redistribution of electron density produced by hyperconjugative mechanisms, which increases the electron population within axial or equatorial proton in normal or reverse effect, respectively. Here an alternative explanation for the Perlin effect is presented on the basis of the topology of the induced current density, which directly determines the nuclear magnetic shielding. Current densities around the C? H bond critical point and intra‐atomic and interatomic contributions to the magnetic shielding explain the observed Perlin effect. The balance between intra‐atomic and interatomic contributions determines the difference in the total atomic shielding. Normal Perlin effect is dominated by intra‐atomic part, whereas reverse effect is dominated by interatomic contribution. © 2015 Wiley Periodicals, Inc.     

The induced current strengths and aromatic pathways of heteroporphyrins and their antiaromatic derivatives

The magnetically induced current strengths as well as nucleus independent chemical shifts of aromatic heteroporphyrins and antiaromatic 22,24‐dideazaheteroporphyrins have been studied using the density functional theory method. The induced current strengths and pathways are obtained by numerical integration of the induced current densities following the specific chemical bonds. The total induced current strengths of antiaromatic 22,24‐dideazaheteroporphyrins is about 6 nA/T weaker than the one for the heteroporphyrins in absolute value. The substitution of pyrrole NH groups by O and S atoms does not change the total induced current strengths. The induced currents around the molecular macroring split at the heterocycles (pyrrole, furan, and thiophene) into the inner and outer routes. The heteroatoms (N, O, and S) have high resistance and consequently lead to a weaker induced current strength than the one passing the outer route in aromatic heteroporphyrins. For antiaromatic 22,24‐dideazaheteroporphyrins, the heteroatoms enhance the current strength and change the main current pathway into the inner route. The induced current strength following the NH moiety is stronger than the one passing the oxygen moiety of furan ring and the sulfur moiety of the thiophene ring in both heteroporphyrins and 22,24‐dideazaheteroporphyrins. © 2015 Wiley Periodicals, Inc.     

Magnetic properties and induced current density in acetylene

Coupled Hartree–Fock perturbation theory, within the framework of an accurate calculation, has been employed to visualize the electronic current density vector field induced by an external uniform magnetic field in the acetylene molecule. The current regime in a plane containing the molecular axis is rather different from those usually accounted for, evidencing the presence of a toroidal vortex in the C? H bond. The maps are useful in rationalizing the experimentally observed proton and carbon NMR chemical shifts.     

Proposed phytic acid standard including a method for its analysis

A method is described that accurately and rapidly quantifies the free and total phosphorous content of a commercially available, purified, phytic acid preparation. This allows its use as a standard for phytic acid determinations in foods. The method involves a wet ashing step followed by phosphorous measurement with a 1-amino-2-naphthol-4-sulfonic acid-molybdate reagent in a microplate reader at 660 nm. The procedure can be performed in 3 h with as little as 50 mg sample.     

Calculation of chemical potentials of chain molecules by the incremental gauge cell method

The gauge cell Monte Carlo method is extended to calculations of the incremental chemical potentials and free energies of linear chain molecules. The method was applied to chains of Lennard-Jones beads with stiff harmonic bonds up to 500 monomers in length. We show that the suggested method quantitatively reproduces the modified Widom particle insertion method of Kumar et al. [S. K. Kumar, I. Szleifer, and A. Z. Panagiotopoulos, Phys. Rev. Lett. 66(22), 2935 (1991)], and is by an order of magnitude more efficient for long chains in terms of the computational time required for the same accuracy of chemical potential calculations. The chain increment ansatz, which suggests that the incremental chemical potential is independent of the chain length, was tested at different temperatures. We confirmed that the ansatz holds only for coils above the θ temperature. Special attention is paid to the effects of the magnitude of adsorption potential and temperature on the behavior of single chains in confinements that are comparable in size with the free chain radius of gyration. At sufficiently low temperatures, the dependence of the incremental chemical potential on the chain length in wetting pores is superficially similar to a capillary condensation isotherm, reflecting monolayer formation following by pore volume filling, as the chain length increases. We find that the incremental gauge cell method is an accurate and efficient technique for calculations of the free energies of chain molecules in bulk systems and nanoconfinements alike. The suggested method may find practical applications, such as modeling polymer partitioning on porous substrates and dynamics of chain translocation into nanopores.     

Superoperator method for calculating configuration interaction,including doubly excited configurations

A superoperator algorithm that had been proposed previously [Teor. Éksp. Khim.,11, No. 1, 3 (1975)] for taking all doubly excited configurations into account has been worked out in detail for singlet states and has been realized numerically in theab initio version of STO-4G. The method does not require the construction of the configuration interaction matrix. This matrix is replaced by superoperators, which are rules for the transformation of the variational parameters, stacked, in a transition-type matrix. Calculations have been performed for small molecules of the second period, and also for linear chains of hydrogen atoms (for h₁₄, 1275 configurations were taken into account), for which additivity of the correlation energy is observed in the approximation that was used.Presented at Conference on Quantum Chemistry, Dnepropetrovsk, September, 1983.Translated from Teoreticheskaya i Éksperimental'naya Khimiya, Vol. 21, No. 4, pp. 385–391, July–August, 1985.     

The calculation of vibrational energy levels of polyatomic molecules including anharmonic effect using contact transformation perturbation method

Using contact transformation perturbation method based on the Taylor expansion of the potential energy function in terms of dimensionless normal coordinates up to sixth‐order, the vibrational energy levels in terms of force constants are derived. The contact transformation theory has been applied to simplify the calculation of perturbation effects. To calculate the second‐order vibrational energy correction, the third and fourth‐order terms of potential function have been placed in the first‐order perturbation Hamiltonian and the second‐order Hamiltonian contains hexatic ones. We present expressions which give relations between the fourth‐ and sixth‐order terms in dimensionless normal coordinates of the potential and the anharmonicity coefficients. For illustration, a set of vibrational energies levels of SO₂, and H₂O molecules including anharmonic effects has been calculated. © 2013 Wiley Periodicals, Inc.