Accurate Hartree-Fock-Slater calculations on small diatomic molecules with the finite-element method

We report on the self-consistent field solution of the Hartree-Fock-Slater equations using the finite-element method for the three small diatomic molecules N_2, BH and CO as examples. The quality of the results is not only better by two orders of magnitude than the fully numerical finite difference method of Laaksonen et al. but the method also requires a smaller number of grid points.

Electronic structure of UF_5

Non-relativistic and relativistic self-consistent Hartree- Fock-Slater and Dirac-Slater models have been used to calculate one-electron energy levels and ionization energies for UF_5. The calculations were performed in an assumed structure of C_4v symmetry with the uranium atom at the center of mass of the molecule. The spacing and level ordering are compared with earlier results obtained with the MS X\alpha method using the muffin-tin approximation. Connections with the multiphoton isotope separation scheme of UF_6 are discussed.

Relativistic atomic level calculations with nuclei carrying additional fractional charges

One-electron energy levels and wavelengths have been calculated for Na-like ions whose nuclei carry quarks with additional charges ±e/3, ±2e/3. The calculations are based on relativistic self-consistent field procedures. The deviations from experimental values exhibit regularities which allow an extrapolation for the wavelengths of 3s - 3p, 3s - 4p, 3p - 3d, and 3p - 4s transitions for the nuclear charge Z = 11± 1/3, ±2/3. A number of transitions are found in the region of visible light which could be used in an optical search for quark atoms.

Electronic configuration in the ground state of atomic lawrencium

Self-consistent relativistic Dirac-Hartree-Fock calculations have been made of some lowlying electronic energies for the atoms of all elements in ground-state ds^2 electron configurations. The results indicate that, contrary to some previous estimates, the ground electronic state of atomic Lr could be in either the 5f^14 6d7s^2 or the 5f^14 7p 7s^2 electron configuration. The separation between the lowest energy level of the 5f^14 6d7s^2 configuration and the lowest energy level of the 5f^14 7p7s^2 configuration is estimated to be (0 ± 3) x 10^3 cm^-1 for atomic Lr.

Relativity and screening effects in heavy-ion collisions

Diatomic correlation diagrams are the main basis for the description of heavy-ion collisions. We have constructed the first realistic relativistic many-electron correlation diagrams based on nonrelativistic self-consistent-field, Hartree-Fock calculations of diatomic molecules plus relativistic corrections. We discuss the relativistic influences as well as the many-electron screening effects in the I-Au system with a combined charge of Z = 132 as an example.

Molecular effects on the cross section for pair production and their importance in the evaluation of the total photonuclear cross section of {^16}O in the {\delta-resonance} region

The screening correction to the coherent pair-production cross section on the oxygen molecule has been calculated using self-consistent relativistic wave functions for the one-center and two-center Coulomb potentials. It is shown that the modification of the wave function due to molecular binding and the interference between contributions from the two atoms have both sizeable effects on the screening correction. The so-obtained coherent pair-production cross section which makes up the largest part of the total atomic cross section was used to evaluate the total nuclear absorption cross section from photon attenuation measurements on liquid oxygen. The result agrees with cross sections for other nuclei if A-scaling is assumed. The molecular effect on the pair cross section amounts to 15 % of the nuclear cross section in the {\delta-resonance} region.

Solution of the time dependent Dirac-Fock-Slater equation for many-electron collisions systems using a time window method

We present a new scheme to solve the time dependent Dirac-Fock-Slater equation (TDDFS) for heavy many electron ion-atom collision systems. Up to now time independent self consistent molecular orbitals have been used to expand the time dependent wavefunction and rather complicated potential coupling matrix elements have been neglected. Our idea is to minimize the potential coupling by using the time dependent electronic density to generate molecular basis functions. We present the first results for 16 MeV S{^16+} on Ar.

Non-relativistic and relativistic molecular calculations for the chalcogen hexafluorides: SF_6, SeF_6, TeF_6, PoF_6

Non-relativistic Hartree-Fock-Slater and relativistic Dirac-Slater self-consistent orbital models are applied for the analysis of the electronic structure of the chalcogen hexafluorides: SF_6, SeF_6, TeF_6 and PoF_6. The molecular eigenfunctions and eigenvalues are generated using the discrete variational method (DVM) with numerical basis functions. The results obtained for SF_6 are compared with other ab initio calculations. Information about relativistic level shifts and spin-orbit splitting has been obtained by comparison between the non-relativistic and relativistic results.

Automatic generation of active coordinates for quantum dynamics calculations: application to the dynamics of benzene photochemistry

A new practical method to generate a subspace of active coordinates for quantum dynamics calculations is presented. These reduced coordinates are obtained as the normal modes of an analytical quadratic representation of the energy difference between excited and ground states within the complete active space self-consistent field method. At the Franck-Condon point, the largest negative eigenvalues of this Hessian correspond to the photoactive modes: those that reduce the energy difference and lead to the conical intersection; eigenvalues close to 0 correspond to bath modes, while modes with large positive eigenvalues are photoinactive vibrations, which increase the energy difference. The efficacy of quantum dynamics run in the subspace of the photoactive modes is illustrated with the photochemistry of benzene, where theoretical simulations are designed to assist optimal control experiments

A variational approach for calculating Franck-Condon factors including mode-mode anharmonic coupling

We have implemented our new procedure for computing Franck-Condon factors utilizing vibrational configuration interaction based on a vibrational self-consistent field reference. Both Duschinsky rotations and anharmonic three-mode coupling are taken into account. Simulations of the first ionization band of Cl O2 and C4 H4 O (furan) using up to quadruple excitations in treating anharmonicity are reported and analyzed. A developer version of the MIDASCPP code was employed to obtain the required anharmonic vibrational integrals and transition frequencies

Finite-stretching corrections to the Milner-Witten-Cates theory for polymer brushes

This paper investigates finite-stretching corrections to the classical Milner-Witten-Cates theory for semi-dilute polymer brushes in a good solvent. The dominant correction to the free energy originates from an entropic repulsion caused by the impenetrability of the grafting surface, which produces a depletion of segments extending a distance $\mu \propto L^{-1}$ from the substrate, where $L$ is the classical brush height. The next most important correction is associated with the translational entropy of the chain ends, which creates the well-known tail where a small population of chains extend beyond the classical brush height by a distance $\xi \propto L^{-1/3}$. The validity of these corrections is confirmed by quantitative comparison with numerical self-consistent field theory.

Polydispersity-induced macrophase separation in diblock copolymer melts

The effect of A-block polydispersity on the phase behavior of AB diblock copolymer melts is examined using a complete self-consistent field theory treatment that allows for fractionation of the parent molecular-weight distribution. In addition to observing the established shift in phase boundaries, we find the emergence of significant two-phase coexistence regions causing, for instance, the disappearance of the complex phase window. Furthermore, we find evidence that polydispersity relieves packing frustration, which will reduce the tendency for long-range order.

Converting the nanodomains of a diblock-copolymer thin film from spheres to cylinders with an external electric field

We investigate the ability of an applied electric field to convert the morphology of a diblock-copolymer thin film from a monolayer of spherical domains embedded in the matrix to cylindrical domains that penetrate through the matrix. As expected, the applied field increases the relative stability of cylindrical domains, while simultaneously reducing the energy barrier that impedes the transition to cylinders. The effectiveness of the field is enhanced by a large dielectric contrast between the two block-copolymer components, particularly when the low-dielectric contrast component forms the matrix. Furthermore, the energy barrier is minimized by selecting sphere-forming diblock copolymers that are as compositionally symmetric as possible. Our calculations, which are the most quantitatively reliable to date, are performed using a numerically precise spectral algorithm based on self-consistent-field theory supplemented with an exact treatment for linear dielectric materials.

Electric field alignment in thin films of cylinder-forming diblock copolymer

We investigate thin films of cylinder-forming diblock copolymer confined between electrically charged parallel plates, using self-consistent-field theory ( SCFT) combined with an exact treatment for linear dielectric materials. Our study focuses on the competition between the surface interactions, which tend to orient cylinder domains parallel to the plates, and the electric field, which favors a perpendicular orientation. The effect of the electric field on the relative stability of the competing morphologies is demonstrated with equilibrium phase diagrams, calculated with the aid of a weak-field approximation. As hoped, modest electric fields are shown to have a significant stabilizing effect on perpendicular cylinders, particularly for thicker films. Our improved SCFT-based treatment removes most of the approximations implemented by previous approaches, thereby managing to resolve outstanding qualitative inconsistencies among different approximation schemes.

Undulation instability in block-copolymer lamellae subjected to a perpendicular electric field

We examine the stability of lamellar stacks in the presence of an electric field, E-0, applied normal to the lamellae. Calculations are performed with self-consistent field theory (SCFT) supplemented by an exact treatment of the electrostatic energy for linear dielectric materials. The calculations identify a critical electric field, E-0*, beyond which the lamellar stack becomes unstable with respect to undulations. This E-0* rapidly decreases towards zero as the number of lamellae in the stack diverges. Our quantitative predictions for E-0* are consistent with previous experimental measurements by Xu and co-workers.