Errors in the application of the JWKB method to calculating the survival factor in dissociative electron attachment using the local complex potential

The accuracy of the JWKB method for determining the survival factor defined for dissociative electron attachment (DEA) processes is examined for a range of electronic resonance lifetimes within the local complex potential approximation. The author concludes that the accuracy is inadequate for molecules with properties commonly found for shape resonance induced DEA. More accurate methods using the uniform Airy function approximation give much better results, but the direct numerical integration of Schrodinger's equation appears simpler still.     

A numerical integration scheme for the evaluation of correlation energy functionals

A numerical integration scheme is presented for three-dimensional integrals occurring in electronic structure calculations, concentrating attention on the evaluation of the correlation energy through a density-functional expression. The scheme is based on the choice of density-based weight functions that naturally partition the space into “atomic” volumes (in which the integration is performed in terms of spherical coordinates) and “diatomic” volumes (in which the integration is performed in terms of confocal elliptical coordinates). Such a choice is justified on the basis of the analytical behavior of the integrand. The attainable accuracy and the required computational effort within the proposed scheme are discussed in detail in a test application on the C₆₀ molecule in the symmetrical configuration. Finally, a comparison with previously proposed schemes is presented. © 1993 John Wiley & Sons, Inc.     

Development of the cyclic cluster model formalism for Kohn-Sham auxiliary density functional theory methods

The development of the cyclic cluster model (CCM) formalism for Kohn-Sham auxiliary density functional theory (KS-ADFT) methods is presented. The CCM is a direct space approach for the calculation of perfect and defective systems under periodic boundary conditions. Translational symmetry is introduced in the CCM by integral weighting. A consistent weighting scheme for all two-center and three-center interactions appearing in the KS-ADFT method is presented. For the first time, an approach for the numerical integration of the exchange-correlation potential within the cyclic cluster formalism is derived. The presented KS-ADFT CCM implementation was applied to covalent periodic systems. The results of cyclic and molecular cluster model (MCM) calculations for trans-polyacetylene, graphene, and diamond are discussed as examples for systems periodic in one, two, and three dimensions, respectively. All structures were optimized. It is shown that the CCM results represent the results of MCM calculations in the limit of infinite molecular clusters. By analyzing the electronic structure, we demonstrate that the symmetry of the corresponding periodic systems is retained in CCM calculations. The obtained geometric and electronic structures are compared with available data from the literature.     

Automatized Parameterization of the Density‐functional Tight‐binding Method. II. Two‐center Integrals

We present an efficient numerical integration scheme (TWOCENT) to be used in the context of automatized parameterization of the density‐functional tight‐binding (DFTB) method. The accuracy of the integration process is assessed and its range of applicability is discussed. The functionality of the developed code is tested by reproducing the electronic portion of the existing mio parameter sets and by reproducing a series of reference DFT band structures of elemental solids.     

Geometric integration in Born-Oppenheimer molecular dynamics

Geometric integration schemes for extended Lagrangian self-consistent Born-Oppenheimer molecular dynamics, including a weak dissipation to remove numerical noise, are developed and analyzed. The extended Lagrangian framework enables the geometric integration of both the nuclear and electronic degrees of freedom. This provides highly efficient simulations that are stable and energy conserving even under incomplete and approximate self-consistent field (SCF) convergence. We investigate three different geometric integration schemes: (1) regular time reversible Verlet, (2) second order optimal symplectic, and (3) third order optimal symplectic. We look at energy conservation, accuracy, and stability as a function of dissipation, integration time step, and SCF convergence. We find that the inclusion of dissipation in the symplectic integration methods gives an efficient damping of numerical noise or perturbations that otherwise may accumulate from finite arithmetics in a perfect reversible dynamics.     

An analysis of the integrand occurring in correlation energy functionals

An analysis of the integrand occurring in current density functional calculations is presented, concentrating attention on correlation energy functionals and the atomic regions, i.e., the regions of space surrounding atomic centers. The analysis follows the structure of a previously proposed numerical integration scheme for three-dimensional integrals occurring in electronic structure calculations. The scheme is based on the choice of density-based weight functions that naturally partition the space into «atomic» volumes (in which the integration is performed in terms of spherical coordinates) and «diatomic» volumes (in which the integration is performed in terms of confocal elliptical coordinates). From this analysis, a simplified procedure for the atomic «internal» integrations is developed, whereas preliminary results are discussed for the atomic and diatomic «external» integrations. The numerical tests are performed on the C₆₀ molecule in the symmetrical configuration. © John Wiley & Sons, Inc.     

On the ground state of Pd13

First-principles electronic structure calculations within a gradient corrected density functional formalism have been carried out to investigate the electronic structure and magnetic properties of Pd(13) clusters. It is shown that a bilayer ground-state structure that can be regarded as a relaxed bulk fragment is most compatible with the experimental results from Stern-Gerlach measurements. An icosahedral structure, considered to be the ground state in numerous previous studies, is shown to be around 0.14 eV above the ground state. A detailed analysis of the molecular orbitals reveals the near degeneracy of the bilayer or icosahedral structures is rooted in the stabilization by p- or d-like cluster orbitals. The importance of low-lying spin states in controlling the electronic and magnetic properties of the cluster is highlighted.     

Modeling the (HI)2 photodissociation dynamics through a nonadiabatic wave packet study of the I*-HI complex

The nonadiabatic photodissociation dynamics of (HI)2 is simulated by applying a wave packet approach which starts from the I*-HI complex (where I* denotes the I(2P1/2) excited electronic state) produced after the photodissociation of the first HI moiety within (HI)2. In the model, two excited electronic potential surfaces corresponding to I*-HI(A 1Pi1) and I-HI(A 1Pi1), which interact through spin-rotation coupling, are considered. The simulations show that upon photodissociation of HI within I*-HI, the dissociating H fragment undergoes intracluster collisions with the I* atom. Some of these collisional events induce an electronically nonadiabatic transition which causes the deactivation of I* to the I ground electronic state. The probability of such nonadiabatic process is found to be 0.37%. Most of the photodissociation process takes place in the upper excited electronic surface [that of the I*-HI(A 1Pi1) complex], where H dissociation is found to be mainly direct or involving weak H/I* intracluster collisions. These weak collisions with high collisional angular momentum, and therefore high collisional impact parameters associated, are responsible for most of the probability of nonadiabatic transitions found. The type of H/I* collisions leading to nonadiabatic transitions appears to be closely related to the nature of the spin-rotation coupling between the two excited electronic states involved.     

Cooperative and Orthogonal Switching in the Solid State Enabled by Metal-Organic Framework Confinement Leading to a Thermo-Photochromic Platform

Cooperative behavior and orthogonal responses of two classes of coordinatively integrated photochromic molecules towards distinct external stimuli were demonstrated on the first example of a photo-thermo-responsive hierarchical platform. Synergetic and orthogonal responses to temperature and excitation wavelength are achieved by confining the stimuli-responsive moieties within a metal–organic framework (MOF), leading to the preparation of a novel photo-thermo-responsive spiropyran-diarylethene based material. Synergistic behavior of two photoswitches enables the study of stimuli-responsive resonance energy transfer as well as control of the photoinduced charge transfer processes, milestones required to advance optoelectronics development. Spectroscopic studies in combination with theoretical modeling revealed a nonlinear effect on the material electronic structure arising from the coordinative integration of photoresponsive molecules with distinct photoisomerization mechanisms. Thus, the reported work covers multivariable facets of not only fundamental aspects of photoswitch cooperativity, but also provides a pathway to modulate photophysics and electronics of multidimensional functional materials exhibiting thermo-photochromism.     

Demonstration of molecular assembly on Si (100) for CMOS-compatible molecule-based electronic devices

In this work, we establish the potential of a UV-promoted direct attachment of alkanes with alcohol and thiol linkers to the silicon (100) surfaces for use in molecular electronic devices with increased potential for integration with existing CMOS technologies. Characterization of the self-assembled monolayers via Fourier transform infrared spectroscopy, spectroscopic ellipsometry, and X-ray photoemission spectroscopy shows that the films assembled on the Si (100) are comparable in quality, aliphatic monolayer coverage, and extent of substrate oxidation to those assembled on the more extensively studied Si (111) crystal face. Simple Si (100)-based electronic devices fabricated with the monolayers exhibited molecule-dependent electrical characteristics. These data highlight the effectiveness of the assembly on Si (100), the ability to fabricate enclosed Si (100)-based molecular devices, and the potential for the future integration of these devices with more conventional technologies.     

Electronic coupling matrix elements from charge constrained density functional theory calculations using a plane wave basis set

We present a plane wave basis set implementation for the calculation of electronic coupling matrix elements of electron transfer reactions within the framework of constrained density functional theory (CDFT). Following the work of Wu and Van Voorhis [J. Chem. Phys. 125, 164105 (2006)], the diabatic wavefunctions are approximated by the Kohn-Sham determinants obtained from CDFT calculations, and the coupling matrix element calculated by an efficient integration scheme. Our results for intermolecular electron transfer in small systems agree very well with high-level ab initio calculations based on generalized Mulliken-Hush theory, and with previous local basis set CDFT calculations. The effect of thermal fluctuations on the coupling matrix element is demonstrated for intramolecular electron transfer in the tetrathiafulvalene-diquinone (Q-TTF-Q(-)) anion. Sampling the electronic coupling along density functional based molecular dynamics trajectories, we find that thermal fluctuations, in particular the slow bending motion of the molecule, can lead to changes in the instantaneous electron transfer rate by more than an order of magnitude. The thermal average, (<|H(ab)|(2)>)(1/2)=6.7 mH, is significantly higher than the value obtained for the minimum energy structure, |H(ab)|=3.8 mH. While CDFT in combination with generalized gradient approximation (GGA) functionals describes the intermolecular electron transfer in the studied systems well, exact exchange is required for Q-TTF-Q(-) in order to obtain coupling matrix elements in agreement with experiment (3.9 mH). The implementation presented opens up the possibility to compute electronic coupling matrix elements for extended systems where donor, acceptor, and the environment are treated at the quantum mechanical (QM) level.     

On-the-fly localization of electronic orbitals in Car-Parrinello molecular dynamics

The ab initio molecular-dynamics formalism of Car and Parrinello is extended to preserve the locality of the orbitals. The supplementary term in the Lagrangian does not affect the nuclear dynamics, but ensures "on the fly" localization of the electronic orbitals within a periodic supercell in the Gamma-point approximation. The relationship between the resulting equations of motion and the formation of a gauge-invariant Lagrangian combined with a gauge-fixing procedure is briefly discussed. The equations of motion can be used to generate a very stable and easy to implement numerical integration algorithm. It is demonstrated that this algorithm can be used to compute the trajectory of the maximally localized orbitals, known as Wannier orbitals, in ab initio molecular dynamics with only a modest increase in the overall computer time. In the present paper, the new method is implemented within the generalized gradient approximation to Kohn-Sham density-functional theory employing plane wave basis sets and atomic pseudopotentials. In the course of the presentation, we briefly discuss how the present approach can be combined with localized basis sets to design fast linear scaling ab initio molecular-dynamics methods.     

Light-controlled organic/inorganic P-N junction nanowires

P-N junctions are of great importance both in modern electronic applications and in understanding other semiconductor devices. Organic/inorganic P-N junction nanowires composed of functional organic molecules and inorganic molecules may be able to realize new or improved chemical and physical properties that were not observed in the individual component on nanosize and their bulk materials. We report herein the fabrication of the organic/inorganic semiconductor P-N junction nanowire and the remarkable performance on the light-controlled diode within a single hybrid P-N junction nanowire. Controlling the conductivity of the P-N junction nanowire by the light irradiation simply to achieve diode work indicates a new way to realize the photoelectric integration in a single nanowire device.     

Optimized Atomic Partial Charges and Radii Defined by Radical Voronoi Tessellation of Bulk Phase Simulations

We present a novel method for the computation of well-defined optimized atomic partial charges and radii from the total electron density. Our method is based on a two-step radical Voronoi tessellation of the (possibly periodic) system and subsequent integration of the total electron density within each Voronoi cell. First, the total electron density is partitioned into the contributions of each molecule, and subsequently the electron density within each molecule is assigned to the individual atoms using a second set of atomic radii for the radical Voronoi tessellation. The radii are optimized on-the-fly to minimize the fluctuation (variance) of molecular and atomic charges. Therefore, our method is completely free of empirical parameters. As a by-product, two sets of optimized atomic radii are produced in each run, which take into account many specific properties of the system investigated. The application of an on-the-fly interpolation scheme reduces discretization noise in the Voronoi integration. The approach is particularly well suited for the calculation of partial charges in periodic bulk phase systems. We apply the method to five exemplary liquid phase simulations and show how the optimized charges can help to understand the interactions in the systems. Well-known effects such as reduced ion charges below unity in ionic liquid systems are correctly predicted without any tuning, empiricism, or rescaling. We show that the basis set dependence of our method is very small. Only the total electron density is evaluated, and thus, the approach can be combined with any electronic structure method that provides volumetric total electron densities—it is not limited to Hartree–Fock or density functional theory (DFT). We have implemented the method into our open-source software tool TRAVIS.     

Resonance raman scattering in the region of a forbidden non-degenarate electronic state by a doubly degenerate vibration

It is shown that the resonance Raman scattering from a forbidden non-degenerate electronic state by a doubly degenerate mode depends strongly on Jahn—Teller interactions within an allowed doubly degenerate electronic state that lends its intensity to the scattering forbidden state.     

Spin-orbit vibronic coupling in 3Pi states of linear triatomic molecules

The Renner-Teller vibronic-coupling problem of a 3Pi electronic state of a linear molecule is analyzed with the inclusion of the spin-orbit coupling of the 3Pi electronic state, employing the microscopic (Breit-Pauli) spin-orbit coupling operator for the two unpaired electrons. The 6x6 Hamiltonian matrix in a diabatic spin-electronic basis is obtained by an expansion of the molecular Hamiltonian in powers of the bending amplitude. The symmetry properties of the Hamiltonian with respect to the time-reversal operator and the relativistic vibronic angular momentum operator are analyzed. It is shown that there exists a linear vibronic-coupling term of spin-orbit origin, which has not been considered so far in the Renner-Teller theory of 3Pi electronic states. While two of the six adiabatic electronic wave functions do not exhibit a geometric phase, the other four carry nontrivial topological phases which depend on the radius of the integration contour. The spectroscopic effects of the linear spin-orbit vibronic-coupling mechanism have been analyzed by numerical calculations of the vibronic spectrum for selected model examples.     

Estimation of electronic coupling in pi-stacked donor-bridge-acceptor systems: correction of the two-state model

Comparison of donor-acceptor electronic couplings calculated within two-state and three-state models suggests that the two-state treatment can provide unreliable estimates of V(da) because of neglecting the multistate effects. We show that in most cases accurate values of the electronic coupling in a pi stack, where donor and acceptor are separated by a bridging unit, can be obtained as V(da) = (E(2)-E(1))mu(12)R(da) + (2E(3)-E(1)-E(2))2mu(13)mu(23)R(da) (2), where E(1), E(2), and E(3) are adiabatic energies of the ground, charge-transfer, and bridge states, respectively, mu(ij) is the transition dipole moments between the states i and j, and R(da) is the distance between the planes of donor and acceptor. In this expression based on the generalized Mulliken-Hush approach, the first term corresponds to the coupling derived within a two-state model, whereas the second term is the superexchange correction accounting for the bridge effect. The formula is extended to bridges consisting of several subunits. The influence of the donor-acceptor energy mismatch on the excess charge distribution, adiabatic dipole and transition moments, and electronic couplings is examined. A diagnostic is developed to determine whether the two-state approach can be applied. Based on numerical results, we showed that the superexchange correction considerably improves estimates of the donor-acceptor coupling derived within a two-state approach. In most cases when the two-state scheme fails, the formula gives reliable results which are in good agreement (within 5%) with the data of the three-state generalized Mulliken-Hush model.     

Three-dimensional orientation of the Qy electronic transition dipole moment within the chlorophyll a molecule determined by femtosecond polarization resolved VIS pump-IR probe spectroscopy

Chlorophyll a (Chl a) is the most abundant pigment on earth. In all plants, algae, and cyanobacteria, it plays a pivotal role as an antenna and reaction center pigment in the primary steps of photosynthesis. In the past, a true three-dimensional (3D) experimental determination of the Qy electronic transition dipole moment orientation could not be obtained. With combined femtosecond polarization resolved VIS pump-IR probe experiments and theoretical calculations of the infrared transition dipole moments (tdm's) in the electronic ground state, we determined the 3D orientation of the Qy electronic tdm of Chl a within the molecular structure. Polarization resolved experiments provided angles of the Qy electronic tdm with three different infrared tdm's, whose orientations within the molecular structure were taken from our theoretical calculations. The orientation of the Qy tdm results from the intersection of all three angles and was found to have an angle of (78 +/- 3)degrees with the x-axis, (12 3)degrees with the y-axis, and (86 +/- 2)degrees with the z-axis.