Accelerating the density-functional tight-binding method using graphical processing units.

Acceleration of the density-functional tight-binding (DFTB) method on single and multiple graphical processing units (GPUs) was accomplished using the MAGMA linear algebra library. Two major computational bottlenecks of DFTB ground-state calculations were addressed in our implementation: the Hamiltonian matrix diagonalization and the density matrix construction. The code was implemented and benchmarked on two different computer systems: (1) the SUMMIT IBM Power9 supercomputer at the Oak Ridge National Laboratory Leadership Computing Facility with 1-6 NVIDIA Volta V100 GPUs per computer node and (2) an in-house Intel Xeon computer with 1-2 NVIDIA Tesla P100 GPUs. The performance and parallel scalability were measured for three molecular models of 1-, 2-, and 3-dimensional chemical systems, represented by carbon nanotubes, covalent organic frameworks, and water clusters.

Dopant-Free Hole-Transport Materials with Germanium Compounds Bearing Pseudohalide and Chalcogenide Moieties for Perovskite Solar Cells.

Hole-transport materials (HTMs) are key electronic components for the functioning of perovskite solar cells (PSCs) as they extract the photogenerated holes from the perovskite to be transported subsequently to the back electrode while minimizing the loss from electron recombination. Herein, we report the synthesis and characterization of novel germanium-based compounds with [{HC(CMeNAr)2}GeNCS] (2), [{HC(CMeNAr)2}Ge(S)NCS] (3), and [{HC(CMeNAr)2}Ge(Se)NCS] (4) compositions, with Ar = 2,6-iPr2C6H3 and the photovoltaic performance of 3 and 4 that is the same as for HTM in PSC. All compounds displayed excellent thermal properties and an appropriate alignment of energy levels for the perovskite with maximum optical absorption in the near-UV region. As revealed by space-charge limited-current (SCLC) measurements, compounds 3 and 4 have competing hole mobilities of about 1.37 × 10-4 and 4.88 × 10-4 cm2 V-1 s-1, respectively. Upon assessing PSC devices using 3 and 4 with triple-cation perovskite absorber Cs0.05(MA0.17FA0.83)0.95Pb(I0.83Br0.17)3, the power conversion efficiencies (PCEs) were about 13.03 and 9.23%, respectively, both without doping and additives, and were compared with benchmark HTM spiro-OMeTAD (2,2',7,7'-tetrakis(N,N-di-p-methoxyphenylamine)-9,9'-spirobifluorene). Quantum chemical calculations with DFT showed that the optoelectronic properties are strongly influenced by the combined contributions of the germanium atom, the pseudohalide moiety (NCS-), and chalcogenides (S2- or Se2-). Fine tuning the electronic properties of germanium is thus a good strategy for the targeted synthesis of potential conducting molecules in PSCs.

Theoretical Prediction and Analysis of the UV/Visible Absorption and Emission Spectra of Chiral Carbon Nanorings.

UV/vis absorption and emission spectra of recently synthesized chiral carbon nanorings were simulated using first-principles-based molecular dynamics and time-dependent density functional theory (TD-DFT). The chiral carbon nanorings are derivatives of the [ n]cycloparaphenylene ([ n]CPP) macrocycles, containing an acene unit such as naphthalene, ([ n]CPPN), anthracene ([ n]CPPA), and tetracene ([ n]CPPT), in addition to n paraphenylene units. In order to study the effect of increasing molecular size on absorption and emission spectra, we investigated the cases where n = 6 and 8. Frontier molecular orbital analysis was carried out to give insight into the degree of excitation delocalization and its relationship to the predicted absorption spectra. The lowest excited singlet state S1 corresponds to a HOMO-LUMO π-π* transition, which is allowed in all chiral carbon nanorings due to lack of molecular symmetry, in contrast to the forbidden HOMO-LUMO transition in the symmetric [ n]CPP molecules. The S1 absorption peak exhibits a blue-shift with increasing number of paraphenylene units in particular for [ n]CPPN and [ n]CPPA and less so in the case of [ n]CPPT. In the case of CPPN and CPPA, the transition density is mainly localized over a semicircle of the macrocycle with the acene unit in its center but is strongly localized on the tetracene unit in the case of CPPT. Molecular dynamics simulations performed on the excited state potential energy surfaces reveal red-shifted emission of these chiral carbon nanorings when the size of the π-conjugated acene units is increased, although the characteristic [ n]CPP blue-shift with increasing paraphenylene unit number n remains apparent. The anomalous emission blue-shift is caused by the excited state bending and torsional motions that stabilize the π HOMO and destabilize the π* LUMO, resulting in an increasing HOMO-LUMO gap.

Light-melt adhesive based on dynamic carbon frameworks in a columnar liquid-crystal phase.

Liquid crystal (LC) provides a suitable platform to exploit structural motions of molecules in a condensed phase. Amplification of the structural changes enables a variety of technologies not only in LC displays but also in other applications. Until very recently, however, a practical use of LCs for removable adhesives has not been explored, although a spontaneous disorganization of LC materials can be easily triggered by light-induced isomerization of photoactive components. The difficulty of such application derives from the requirements for simultaneous implementation of sufficient bonding strength and its rapid disappearance by photoirradiation. Here we report a dynamic molecular LC material that meets these requirements. Columnar-stacked V-shaped carbon frameworks display sufficient bonding strength even during heating conditions, while its bonding ability is immediately lost by a light-induced self-melting function. The light-melt adhesive is reusable and its fluorescence colour reversibly changes during the cycle, visualizing the bonding/nonbonding phases of the adhesive.

Hybridization of a flexible cyclooctatetraene core and rigid aceneimide wings for multiluminescent flapping π systems.

The hybridization of flexible and rigid π-conjugated frameworks is a potent concept for producing new functional materials. In this article, a series of multifluorescent flapping π systems that combine a flexible cyclooctatetraene (COT) core and rigid aceneimide wings with various π-conjugation lengths has been designed and synthesized, and their structure/properties relationships have been investigated. Whereas these molecules have a V-shaped bent conformation in the ground state, the bent structure changes to a planar conformation in the lowest excited singlet (S1 ) state irrespective of the lengths of the aceneimide wings. However, the fluorescence behavior in solution is distinct between the naphthaleneimide system and the anthraceneimide system. The former has a nonemissive S1 state owing to the significant contribution of the antiaromatic character of the planar COT frontier molecular orbitals, thereby resulting in complete fluorescence quenching in solution. In contrast, the latter anthraceneimide system shows an intense emission, which is ascribed to the planar but distorted S1 state that shows the allowed transition between the π-molecular orbitals delocalized over the COT core and the acene wings. The other characteristic of these π systems is the significantly redshifted fluorescence in the crystalline state relative to their monomer fluorescence. The relationship between the packing structures and the fluorescence properties was investigated by preparing a series of hybrid π systems with different sizes of substituents on the imide moieties, which revealed the effect of the twofold π-stacked structure of the V-shaped molecules on the large bathochromic shift in emission.

Quantum Dynamics Simulations Reveal Vibronic Effects on the Optical Properties of [n]Cycloparaphenylenes.

The size-dependent ultraviolet/visible photophysical property trends of [n]cycloparaphenylenes ([n]CPPs, n = 6, 8, and 10) are theoretically investigated using quantum dynamics simulations. For geometry optimizations on the ground- and excited-state Born-Oppenheimer potential energy surfaces (PESs), we employ density functional theory (DFT) and time-dependent DFT calculations. Harmonic normal-mode analyses are carried out for the electronic ground state at Franck-Condon geometries. A diabatic Hamiltonian, comprising four low-lying singlet excited electronic states and 26 vibrational degrees of freedom of CPP, is constructed within the linear vibronic coupling (VC) model to elucidate the absorption spectral features in the range of 300-500 nm. Quantum nuclear dynamics is simulated within the multiconfiguration time-dependent Hartree approach to calculate the vibronic structure of the excited electronic states. The symmetry-forbidden S0 → S1 transition appears in the longer wavelength region of the spectrum with weak intensity due to VC. It is found that the Jahn-Teller and pseudo-Jahn-Teller effects in the doubly degenerate S2 and S3 electronic states are essential in the quantitative interpretation of the experimental observation of a broad absorption peak around 340 nm. The vibronic mixing of the S1 state with higher electronic states is responsible for the efficient photoluminescence from the S1 state. The fluorescence properties are characterized on the basis of the stationary points of the excited-state PESs. The findings reveal that vibronic effects become important in determining the photophysical properties of CPPs with increased ring size.

Photochemical double 5-exo cyclization of alkenyl-substituted dithienylacetylenes: efficient synthesis of diarylated dithienofulvalenes.

Smooth and selective: Upon photoirradiation, bis(3-alkenyl-2-thienyl)acetylenes smoothly and selectively undergo double 5-exo-dig cyclization to produce a series of thiophene-fused pentafulvalenes with various aryl substituents. In this fused π-conjugated skeleton, the fused thiophene rings and the aryl substituents significantly modulate the electronic structure of the pentafulvalene skeleton.

A π-conjugated system with flexibility and rigidity that shows environment-dependent RGB luminescence.

We have designed and synthesized a π-conjugated system that consists of a flexible and nonplanar π joint and two emissive rigid and planar wings. This molecular system exhibits respectively red, green, and blue (RGB) emission from a single-component luminophore in different environments, namely in polymer matrix, in solution, and in crystals. The flexible unit gives rise to a dynamic conformational change in the excited state from a nonplanar V-shaped structure to a planar structure, leading to a dual fluorescence of blue and green colors. The rigid and planar moieties favor the formation of a two-fold π-stacked array of the V-shaped molecules in the crystalline state, which produces a red excimer-like emission. These RGB emissions are attained without changing the excitation energy.

Infrared absorption of methanol clusters (CH3OH)n with n = 2-6 recorded with a time-of-flight mass spectrometer using infrared depletion and vacuum-ultraviolet ionization.

We investigated IR spectra in the CH- and OH-stretching regions of size-selected methanol clusters, (CH(3)OH)(n) with n = 2-6, in a pulsed supersonic jet by using the IR-VUV (vacuum-ultraviolet) ionization technique. VUV emission at 118 nm served as the source of ionization in a time-of-flight mass spectrometer. The tunable IR laser emission served as a source of predissociation or excitation before ionization. The variations of intensity of protonated methanol cluster ions (CH(3)OH)(n)H(+) and CH(3)OH(+) and (CH(3)OH)(2)(+) were monitored as the IR laser light was tuned across the range 2650-3750 cm(-1). Careful processing of these action spectra based on photoionization efficiencies and the production and loss of each cluster due to photodissociation yielded IR spectra of the size-selected clusters. Spectra of methanol clusters in the OH region have been extensively investigated; our results are consistent with previous reports, except that the band near 3675 cm(-1) is identified as being associated with the proton acceptor of (CH(3)OH)(2). Spectra in the CH region are new. In the region 2800-3050 cm(-1), bands near 2845, 2956, and 3007 cm(-1) for CH(3)OH split into 2823, 2849, 2934, 2955, 2984, and 3006 cm(-1) for (CH(3)OH)(2) that correspond to proton donor and proton acceptor, indicating that the methanol dimer has a preferred open-chain structure. In contrast, for (CH(3)OH)(3), the splitting diminishes and the bands near 2837, 2954, and 2987 cm(-1) become narrower, indicating a preferred cyclic structure. Anharmonic vibrational wavenumbers predicted for the methanol open-chain dimer and the cyclic trimer with the B3LYP∕VPT2∕ANO1 level of theory are consistent with experimental results. For the tetramer and pentamer, the spectral pattern similar to that of the trimer but with greater widths was observed, indicating that the most stable structures are also cyclic.

Theoretical interpretation of the UV-vis spectrum of the CS2/Cl complex in the spectral region 320-550 nm.

Accurate multireference configuration interaction and time-dependent density functional calculations have been performed to interpret the experimental UV-vis spectrum of the CS(2)/Cl complex in the spectral region 320-550 nm. The molecular structure of the complex responsible for the previously observed UV-vis spectrum is recognized as ClSCS, not ClCS(2). Two low-lying excited states of ClSCS, responsible for its optical absorption, have been identified and analyzed. Optical excitation of ClSCS leads to the excitation-specific bond elongation that may lead to photofragmentation of the molecule. In addition, experimental conditions for verifying the presence of ClCS(2) are identified and detailed characterization of its optically active excited states with possible photofragmentation pathways is given.

The low-lying states of the scandium dimer.

A systematic investigation of low-lying states of Sc(2) using multireference perturbation theory (NEVPT2 and NEVPT3) indicates that the ground state of this system is (5)Sigma(u) (-) with r(e)=2.611 A, omega(e)=241.8 cm(-1), and D(e)=1.78 eV. This state is closely followed by other low-lying states of Sc(2): (3)Sigma(u) (-), (5)Delta(u), (3)Pi(g), (1)Pi(g), and (1)Sigma(u) (-). Our energy ordering of the (5)Sigma(u) (-) and (3)Sigma(u) (-) states confirms the recent MRCI results of Kalemos et al. [J. Chem. Phys. 132, 024309 (2010)] and is at variance with the earlier diffusion Monte Carlo predictions of Matxain et al. [J. Chem. Phys. 128, 194315 (2008)]. An excellent agreement between the second- and third-order NEVPT results and between the computed and experimental values of omega(e) (241.8 versus 238.9 cm(-1)) for the (5)Sigma(u) (-) state suggests high accuracy of our predictions.

Multireference perturbation theory can predict a false ground state.

Prediction of a false ground state with popular variants of multireference perturbation theory (CASPT2 and MRMP) is reported for a remarkably simple chemical system: the Sc(2) molecule.

Intruder states in multireference perturbation theory: the ground state of manganese dimer.

A detailed analysis of a severe intruder state problem in the multistate multireference perturbation theory (MS-MRPT) calculations on the ground state of manganese dimer is presented. An enormous number of detected intruder states (> 5000) do not permit finding even an approximate shape of the X(1)Sigma(g) (+) potential energy curve. The intruder states are explicitly demonstrated to originate from quasidegeneracies in the zeroth-order Hamiltonian spectrum. The electronic configurations responsible for appearance of the quasidegeneracies are identified as single and double excitations from the active orbitals to the external orbitals. It is shown that the quasidegeneracy problem can be completely eliminated using shift techniques despite of its severity. The resultant curves are smooth and continuous. Unfortunately, strong dependence of the spectroscopic parameters of the X(1)Sigma(g) (+) state on the shift parameter is observed. This finding rises serious controversies regarding validity of employing shift techniques for solving the intruder state problem in MS-MRPT. Various alternative approaches of removing intruder states (e.g., modification of the basis set or changing the active space) are tested. None of these conventional techniques is able to fully avoid the quasidegeneracies. We believe that the MS-MRPT calculations on the three lowest A(g) states of manganese dimer constitute a perfect benchmark case for studying the behavior of MRPT in extreme situations.

Choosing a proper complete active space in calculations for transition metal dimers: ground state of Mn2 revisited.

The potential energy curve of the ground state of Mn(2) has been studied using a systematic sequence of complete active spaces. Deficiencies of the routinely used active space, built from atomic 4s and 3d orbitals, has been identified and discussed. It is shown that an additional sigma(g) orbital, originating from the atomic virtual 4p(z) orbitals, is essential for a proper description of static correlation in the (1)Sigma(g)(+) state of Mn(2). The calculated spectroscopic parameters of the (1)Sigma(g)(+) state agree well with available experimental data. The calculated equilibrium bond lengths are located between 3.24 and 3.50 A, the harmonic vibrational frequencies, between 44 and 72 cm(-1), and the dissociation energies, between 0.05 and 0.09 eV. An urgent need for an accurate gas-phase experimental study of spectroscopic constants of Mn(2) is highlighted.