Strongly Bound Excitons and Anisotropic Linear Absorption in Monolayer Graphullerene.

Graphullerene is a novel two-dimensional carbon allotrope with unique optoelectronic properties. Despite significant experimental characterization and prior density functional theory calculations, unanswered questions remain as to the nature, energy, and intensity of the electronic and optical excitations. Here, we present first-principles calculations of the quasiparticle band structure, neutral excitations, and absorption spectra of monolayer graphullerene and bulk graphullerite, employing the GW-Bethe-Salpeter equation (GW-BSE) approach. We show that strongly bound excitons dominate the absorption spectrum of monolayer graphullerene with binding energies up to 0.8 eV, while graphullerite exhibits less pronounced excitonic effects. Our calculations also reveal a strong linear polarization anisotropy, reflecting the in-plane structural anisotropy from intermolecular coupling between neighboring C60 units. We further show that the presence of Mg atoms, crucial to the synthesis process, induces structural modifications and polarizability effects, resulting in a ∼1 eV quasiparticle gap renormalization and a reduction in the exciton binding energy to ∼0.6 eV.

Polarization Saturation in Multilayered Interfacial Ferroelectrics.

Van der Waals polytypes of broken inversion and mirror symmetries  have been recently shown to exhibit switchable electric polarization even at the ultimate two-layer thin limit. Their out-of-plane polarization has been found to accumulate in a ladder-like fashion with each successive layer, offering 2D building blocks for the bottom-up construction of 3D ferroelectrics. Here, it is demonstrated experimentally that beyond a critical stack thickness, the accumulated polarization in rhombohedral polytypes of molybdenum disulfide saturates. The underlying saturation mechanism, deciphered via density functional theory and self-consistent Poisson-Schrödinger calculations, point to a purely electronic redistribution involving: 1. Polarization-induced bandgap closure that allows for cross-stack charge transfer and the emergence of free surface charge; 2. Reduction of the polarization saturation value, as well as the critical thickness at which it is obtained, by the presence of free carriers. The resilience of polar layered structures to atomic surface reconstruction, which is essentially unavoidable in polar 3D crystals, potentially allows for the design of new devices with mobile surface charges. The findings, which are of general nature, should be accounted for when designing switching and/or conductive devices based on ferroelectric layered materials.

Optically Triggered Emergent Mesostructures in Monolayer WS2.

The ultrahigh surface area of two-dimensional materials can drive multimodal coupling between optical, electrical, and mechanical properties that leads to emergent dynamical responses not possible in three-dimensional systems. We observed that optical excitation of the WS2 monolayer above the exciton energy creates symmetrically patterned mechanical protrusions which can be controlled by laser intensity and wavelength. This observed photostrictive behavior is attributed to lattice expansion due to the formation of polarons, which are charge carriers dressed by lattice vibrations. Scanning Kelvin probe force microscopy measurements and density functional theory calculations reveal unconventional charge transport properties such as the spatially and optical intensity-dependent conversion in the WS2 monolayer from apparent n- to p-type and the subsequent formation of effective p-n junctions at the boundaries between regions with different defect densities. The strong opto-electrical-mechanical coupling in the WS2 monolayer reveals previously unexplored properties, which can lead to new applications in optically driven ultrathin microactuators.

Surface-Guided Crystallization of Xanthine Derivatives for Optical Metamaterial Applications.

Numerous bio-organisms employ template-assisted crystallization of molecular solids to yield crystal morphologies with unique optical properties that are difficult to reproduce synthetically. Here, a facile procedure is presented to deposit bio-inspired birefringent crystals of xanthine derivatives on a template of single-crystal quartz. Crystalline sheets that are several millimeters in length, several hundred micrometers in width, and 300-600 nm thick, are obtained. The crystal sheets are characterized with a well-defined orientation both in and out of the substrate plane, giving rise to high optical anisotropy in the plane parallel to the quartz surface, with a refractive index difference Δn ≈ 0.25 and a refractive index along the slow axis of n ≈ 1.7. It is further shown that patterning of the crystalline stripes with a tailored periodic grating leads to a thin organic polarization-dependent diffractive meta-surface, opening the door to the fabrication of various optical devices from a platform of small-molecule based organic dielectric crystals.

Anisotropic Interlayer Force Field for Group-VI Transition Metal Dichalcogenides.

An anisotropic interlayer force field that describes the interlayer interactions in homogeneous and heterogeneous interfaces of group-VI transition metal dichalcogenides (MX2, where M = Mo, W, and X = S, Se) is presented. The force field is benchmarked against density functional theory calculations for bilayer systems within the Heyd-Scuseria-Ernzerhof hybrid density functional approximation, augmented by a nonlocal many-body dispersion treatment of long-range correlation. The parametrization yields good agreement with the reference calculations of binding energy curves and sliding potential energy surfaces. It is found to be transferable to transition metal dichalcogenide (TMD) junctions outside of the training set that contain the same atom types. Calculated bulk moduli agree with most previous dispersion-corrected density functional theory predictions, which underestimate the available experimental values. Calculated phonon spectra of the various junctions under consideration demonstrate the importance of appropriately treating the anisotropic nature of the layered interfaces. Considering our previous parametrization for MoS2, the anisotropic interlayer potential enables accurate and efficient large-scale simulations of the dynamical, tribological, and thermal transport properties of a large set of homogeneous and heterogeneous TMD interfaces.

Spin state and magnetic coupling in polynuclear Ni(II) complexes from density functional theory: is there an optimal amount of Fock exchange?

Reliable prediction of the ground-state spin and magnetic coupling constants in transition-metal complexes is a well-known challenge for density functional theory (DFT). One popular strategy for addressing this long-standing issue involves the modification of the fraction of Fock exchange in a hybrid functional. Here we explore the viability of this approach using three polynuclear metal-organic complexes based on a Ni4O4 cubane motif, having different ground state spin values (S = 0, 2, 4) owing to the use of different ligands. We systematically search for an optimum fraction of Fock exchange, across various global, range-separated, and double hybrid functionals. We find that for all functionals tested, at best there only exists a very narrow range of Fock exchange fractions which results in a correct prediction of the ground-state spin for all three complexes. The useful range is functional dependent, but general trends can be identified. Typically, at least two similar systems must be used in order to determine both an upper and lower limit of the optimal range. This is likely owing to conflicting demands of minimizing delocalization errors, which typically requires a higher percentage of Fock exchange, and addressing static correlation, which typically requires a lower one. Furthermore, we find that within the optimal range of Fock exchange, the sign and relative magnitude of Ni-Ni magnetic coupling constants are reasonably well reproduced, but there is still room for quantitative improvement in the prediction. Thus, the prediction of spin state and magnetic coupling in polynuclear complexes remains an ongoing challenge for DFT.

Finite Difference Interpolation for Reduction of Grid-Related Errors in Real-Space Pseudopotential Density Functional Theory.

The real-space pseudopotential approach is a well-known method for large-scale density functional theory (DFT) calculations. One of its main limitations, however, is the introduction of errors associated with the positioning of the underlying real-space grid, a phenomenon usually known as the "egg-box" effect. The effect can be controlled by using a finer grid, but this raises the cost of the calculations or even undermines their feasibility altogether. Therefore, there is ongoing interest in the reduction of the effect per a given real-space grid. Here, we present a finite difference interpolation of electron orbitals as a means of exploiting the high resolution of the pseudopotential to reduce egg-box effects systematically. We implement the method in PARSEC, a finite difference real-space pseudopotential DFT code, and demonstrate error mitigation and improved convergence at a low additional computational cost.

Uncovering Multiple Intrinsic Chiral Phases in (PEA)2PbI4 Halide Perovskites.

Two-dimensional (2D) halide perovskites offer a unique platform for investigating the ground state of materials possessing significant anharmonicity. In contrast to three-dimensional perovskites, their 2D counterparts offer substantially fewer degrees of freedom, resulting in multiple well-defined crystal structures. In this work, we thoroughly investigate the anharmonic ground state of the benchmark (PEA)2PbI4 compound, using complementary information from low-temperature X-ray diffraction (XRD) and photoluminescence spectroscopy, supported by density functional theory calculations. We extrapolate four crystallographic configurations from low-temperature XRD. These configurations imply that the ground state has an intrinsic disorder stemming from two coexisting chiral sublattices, each with a bioriented organic spacer molecule. We further show evidence that these chiral structures form unevenly populated ground states, portraying uneven anharmonicity, where the state population may be tuned by surface effects. Our results uncover a disordered ground state that may induce intrinsic grain boundaries, which cannot be ignored in practical applications.

Electrochemical Doping of Halide Perovskites by Noble Metal Interstitial Cations.

Metal halide perovskites are an attractive class of semiconductors, but it has proven difficult to control their electronic doping by conventional strategies due to screening and compensation by mobile ions or ionic defects. Noble-metal interstitials represent an under-studied class of extrinsic defects that plausibly influence many perovskite-based devices. In this work, doping of metal halide perovskites is studied by electrochemically formed Au+ interstitial ions, combining experimental data on devices with a computational analysis of Au+ interstitial defects based on density functional theory (DFT). Analysis suggests that Au+ cations can be easily formed and migrate through the perovskite bulk via the same sites as iodine interstitials (Ii + ). However, whereas Ii + compensates n-type doping by electron capture, the noble-metal interstitials act as quasi-stable n-dopants. Experimentally, voltage-dependent, dynamic doping by current density-time (J-t), electrochemical impedance, and photoluminescence measurements are characterized. These results provide deeper insight into the potential beneficial and detrimental impacts of metal electrode reactions on long-term performance of perovskite photovoltaic and light-emitting diodes, as well as offer an alternative doping explanation for the valence switching mechanism of halide-perovskite-based neuromorphic and memristive devices.

Electronic Coupling and Electrocatalysis in Redox Active Fused Iron Corroles.

Conjugated arrays composed of corrole macrocycles are increasingly more common, but their chemistry still lags behind that of their porphyrin counterparts. Here, we report on the insertion of iron(III) into a ß,ß-fused corrole dimer and on the electronic effects that this redox active metal center has on the already rich coordination chemistry of [H3tpfc] COT, where COT = cyclo-octatetraene and tpfc = tris(pentafluorophenyl)corrole. Synthetic manipulations were performed for the isolation and full characterization of both the 5-coordinate [FeIIItpfc(py)]2COT and 6-coordinate [FeIIItpfc(py)2]2COT, with one and two axial pyridine ligands per metal, respectively. X-Ray crystallography reveals a dome-shaped structure for [FeIIItpfc(py)]2COT and a perfectly planar geometry which (surprisingly at first) is also characterized by shorter Fe-N (corrole) and Fe-N (pyridine) distances. Computational investigations clarify that the structural phenomena are due to a change in the iron(III) spin state from intermediate (S = 3/2) to low (S = 1/2), and that both the 5- and 6-coordinated complexes are enthalpically favored. Yet, in contrast to iron(III) porphyrins, the formation enthalpy for the coordination of the first pyridine to Fe(III) corrole is more negative than that of the second pyridine coordination. Possible interactions between the two corrole subunits and the chelated iron ions were examined by UV-Vis spectroscopy, electrochemical techniques, and density functional theory (DFT). The large differences in the electronic spectra of the dimer relative to the monomer are concluded to be due to a reduced electronic gap, owing to the extensive electron delocalization through the fusing bridge. A cathodic sweep for the dimer discloses two redox processes, separated by 230 mV. The DFT self-consistent charge density for the neutral and cationic states (1- and 2-electron oxidized) reveals that the holes are localized on the macrocycle. A different picture emerges from the reduction process, where both the electrochemistry and the calculated charge density point toward two consecutive electron transfers with similar energetics, indicative of very weak electron communication between the two redox active iron(III) sites. The binuclear complex was determined to be a much better catalyst for the electrochemical hydrogen evolution reaction (HER) than the analogous mononuclear corrole.

Cumulative polarization in conductive interfacial ferroelectrics.

Ferroelectricity in atomically thin bilayer structures has been recently predicted1 and measured2-4 in two-dimensional materials with hexagonal non-centrosymmetric unit-cells. The crystal symmetry translates lateral shifts between parallel two-dimensional layers to sign changes in their out-of-plane electric polarization, a mechanism termed 'slide-tronics'4. These observations have been restricted to switching between only two polarization states under low charge carrier densities5-12, limiting the practical application of the revealed phenomena13. To overcome these issues, one should explore the nature of polarization in multi-layered van der Waals stacks, how it is governed by intra- and interlayer charge redistribution and to what extent it survives the addition of mobile charge carriers14. To explore these questions, we conduct surface potential measurements of parallel WSe2 and MoS2 multi-layers with aligned and anti-aligned configurations of the polar interfaces. We find evenly spaced, nearly decoupled potential steps, indicating highly confined interfacial electric fields that provide a means to design multi-state 'ladder-ferroelectrics'. Furthermore, we find that the internal polarization remains notable on electrostatic doping of mobile charge carrier densities as high as 1013 cm-2, with substantial in-plane conductivity. Using density functional theory calculations, we trace the extra charge redistribution in real and momentum spaces and identify an eventual doping-induced depolarization mechanism.

DFT exchange: sharing perspectives on the workhorse of quantum chemistry and materials science.

In this paper, the history, present status, and future of density-functional theory (DFT) is informally reviewed and discussed by 70 workers in the field, including molecular scientists, materials scientists, method developers and practitioners. The format of the paper is that of a roundtable discussion, in which the participants express and exchange views on DFT in the form of 302 individual contributions, formulated as responses to a preset list of 26 questions. Supported by a bibliography of 777 entries, the paper represents a broad snapshot of DFT, anno 2022.

Real-Space Crystal Structure Analysis by Low-Dose Focal-Series TEM Imaging of Organic Materials with Near-Atomic Resolution.

Structural analysis of beam-sensitive materials by transmission electron microscopy (TEM) represents a significant challenge, as high-resolution TEM (HRTEM) requires high electron doses that limit its applicability to stable inorganic materials. Beam-sensitive materials, e.g., organic crystals, must be imaged under low dose conditions, leading to problematic contrast interpretation and loss of fine structural details. Here, HRTEM imaging of organic crystalline materials with near-atomic resolution of up to 1.6 Å is described, which enables real-space studies of crystal structures, as well as observation of co-existing polymorphs, crystal defects, and atoms. This is made possible by a low-dose focal-series reconstruction methodology, which provides HRTEM images where contrast reflects true object structure and can be performed on contemporary cryo-EM instruments available to many research institutions. Copper phthalocyanine (CuPc), a perchlorinated analogue of CuPc, and indigo crystalline films are imaged. In the case of indigo crystals, co-existing polymorphs and individual atoms (carbonyl oxygen) can be observed. In the case of CuPc, several polymorphs are observed, including a new one, for which the crystal structure is found based on direct in-focus imaging, accomplishing real-space crystal structure elucidation. Such direct analysis can be transformative for structure studies of organic materials.

Optimal Tuning Perspective of Range-Separated Double Hybrid Functionals.

We study the optimal tuning of the free parameters in range-separated double hybrid functionals, based on enforcing the exact conditions of piecewise linearity and spin constancy. We find that introducing the range separation in both the exchange and the correlation terms allows for the minimization of both fractional charge and fractional spin errors for singlet atoms. The optimal set of parameters is system specific, underlining the importance of the tuning procedure. We test the performance of the resulting optimally tuned functionals for the dissociation curves of diatomic molecules. We find that they recover the correct dissociation curve for the one-electron system, H2+, and improve the dissociation curves of many-electron molecules such as H2 and Li2, but they also yield a nonphysical maximum and only converge to the correct dissociation limit at very large distances.

Single Excitation Energies Obtained from the Ensemble "HOMO-LUMO Gap": Exact Results and Approximations.

In calculations based on density functional theory, the "HOMO-LUMO gap" (difference between the highest occupied and lowest unoccupied molecular orbital energies) is often used as a low-cost, ad hoc approximation for the lowest excitation energy. Here we show that a simple correction based on rigorous ensemble density functional theory makes the HOMO-LUMO gap exact in principle and significantly more accurate in practice. The introduced perturbative ensemble density functional theory approach predicts different and useful values for singlet-singlet and singlet-triplet excitations, using semilocal and hybrid approximations. Excitation energies are similar in quality to time-dependent density functional theory, especially at high fractions of exact exchange. The approach therefore offers an easy-to-implement and low-cost route to robust prediction of molecular excitation energies.

Polymorphism, Structure, and Nucleation of Cholesterol·H2O at Aqueous Interfaces and in Pathological Media: Revisited from a Computational Perspective.

We revisit the important issues of polymorphism, structure, and nucleation of cholesterol·H2O using first-principles calculations based on dispersion-augmented density functional theory. For the lesser known monoclinic polymorph, we obtain a fully extended H-bonded network in a structure akin to that of hexagonal ice. We show that the energy of the monoclinic and triclinic polymorphs is similar, strongly suggesting that kinetic and environmental effects play a significant role in determining polymorph nucleation. Furthermore, we find evidence in support of various O-H···O bonding motifs in both polymorphs that may result in hydroxyl disorder. We have been able to explain, via computation, why a single cholesterol bilayer in hydrated membranes always crystallizes in the monoclinic polymorph. We rationalize what we believe is a single-crystal to single-crystal transformation of the monoclinic form on increased interlayer growth beyond that of a single cholesterol bilayer, interleaved by a water bilayer. We show that the ice-like structure is also relevant to the related cholestanol·2H2O and stigmasterol·H2O crystals. The structure of stigmasterol hydrate both as a trilayer film at the air-water interface and as a macroscopic crystal further assists us in understanding the polymorphic and thermal behavior of cholesterol·H2O. Finally, we posit a possible role for one of the sterol esters in the crystallization of cholesterol·H2O in pathological environments, based on a composite of a crystalline bilayer of cholesteryl palmitate bound epitaxially as a nucleating agent to the monoclinic cholesterol·H2O form.

Theory of Chirality Induced Spin Selectivity: Progress and Challenges.

A critical overview of the theory of the chirality-induced spin selectivity (CISS) effect, that is, phenomena in which the chirality of molecular species imparts significant spin selectivity to various electron processes, is provided. Based on discussions in a recently held workshop, and further work published since, the status of CISS effects-in electron transmission, electron transport, and chemical reactions-is reviewed. For each, a detailed discussion of the state-of-the-art in theoretical understanding is provided and remaining challenges and research opportunities are identified.

Anisotropic Interlayer Force Field for Transition Metal Dichalcogenides: The Case of Molybdenum Disulfide.

An anisotropic interlayer force field that describes the interlayer interactions in molybdenum disulfide (MoS2) is presented. The force field is benchmarked against density functional theory calculations for both bilayer and bulk systems within the Heyd-Scuseria-Ernzerhof hybrid density functional approximation, augmented by a nonlocal many-body dispersion treatment of long-range correlation. The parametrization yields good agreement with the reference calculations of binding energy curves and sliding potential energy surfaces for both bilayer and bulk configurations. Benchmark calculations for the phonon spectra of bulk MoS2 provide good agreement with experimental data, and the calculated bulk modulus falls in the lower part of experimentally measured values. This indicates the accuracy of the interlayer force field near equilibrium. Under external pressures up to 20 GPa, the developed force field provides a good description of compression curves. At higher pressures, deviations from experimental data grow, signifying the validity range of the developed force field.

Band gaps of crystalline solids from Wannier-localization-based optimal tuning of a screened range-separated hybrid functional.

Accurate prediction of fundamental band gaps of crystalline solid-state systems entirely within density functional theory is a long-standing challenge. Here, we present a simple and inexpensive method that achieves this by means of nonempirical optimal tuning of the parameters of a screened range-separated hybrid functional. The tuning involves the enforcement of an ansatz that generalizes the ionization potential theorem to the removal of an electron from an occupied state described by a localized Wannier function in a modestly sized supercell calculation. The method is benchmarked against experiment for a set of systems ranging from narrow band-gap semiconductors to large band-gap insulators, spanning a range of fundamental band gaps from 0.2 to 14.2 electronvolts (eV), and is found to yield quantitative accuracy across the board, with a mean absolute error of ∼0.1 eV and a maximal error of ∼0.2 eV.

Space-Filling Curves for Real-Space Electronic Structure Calculations.

Hamiltonian matrices for Kohn-Sham calculations implemented in real space are often large (millions by millions) but very sparse. This poses challenges and opportunities for iterative eigensolvers, which often require a large number of matrix-vector multiplications. As a consequence, an efficient parallel sparse matrix-vector multiplication algorithm is desired. Here, we investigate the benefits of using Hilbert space-filling curves (SFCs) in domain partitioning. We show that the use of Hilbert SFCs in grid-point partitioning brings better locality of the grid points, improves balance of communication, and reduces communication overhead. We also demonstrate an extension of Hilbert SFCs coupled with blockwise operations. The use of blockwise operations helps exploit the vector-processing units in contemporary computational platforms. We illustrate speedup and scalability improvements for an iterative eigensolver based on the Chebyshev-filtered subspace iteration method. Using blockwise Hilbert SFCs, we solve the Kohn-Sham problem for silicon nanocrystals up to 10 nm in diameter, which contain over 26,000 atoms. We illustrate how the density of states of silicon nanocrystals evolves to the bulk limit, where Van Hove singularities are clearly apparent.