Hückel Molecular Orbital Analysis for Stability and Instability of Stacked Aromatic and Stacked Antiaromatic Systems

Face-to-face stacking of aromatic compounds leads to stacked antiaromaticity, while that of antiaromatic compounds leads to stacked aromaticity. This is a prediction with a long history; in the late 2000s, the prediction was confirmed by high-precision quantum chemical calculations, and finally, in 2016, a π-conjugated system with stacked aromaticity was synthesized. Several variations have since been reported, but essentially, they are all the same molecule. To realize stacked aromaticity in a completely new and different molecular system and to trigger an extension of the concept of stacked aromaticity, it is important to understand the origin of stacked aromaticity. The Hückel method, which has been successful in giving qualitatively correct results for π-conjugated systems despite its bold assumptions, is well suited for the analysis of stacked aromaticity. We use this method to model the face-to-face stacking systems of benzene and cyclobutadiene molecules and discuss their stacked antiaromaticity and stacked aromaticity on the basis of their π-electron energies. By further developing the discussion, we search for clues to realize stacked aromaticity in synthesizable molecular systems

Frontier Orbital Views of Stacked Aromaticity

Recent studies have theoretically and experimentally demonstrated that antiaromatic molecules with 4n π electrons exhibit stacked aromaticity according to π–π stacking when arranged in a face-to-face manner. However, the mechanism of its occurrence has not been clearly studied. In this study, we investigated the mechanism of stacked aromaticity using cyclobutadiene. When the antiaromatic molecules are stacked in a face-to-face manner, the orbital interactions between the degenerate singly occupied molecular orbitals (SOMOs) of the monomer unit cause a larger energy gap between the degenerate highest-occupied molecular orbitals (HOMOs) and the lowest-unoccupied molecular orbitals (LUMOs) of the dimer. However, the antiaromatic molecules are more stable in less symmetric conformations, mainly because of pseudo-Jahn–Teller distortions. In the case of cyclobutadiene, the two SOMOs of the monomer unit split into HOMO and LUMO because of the bond alternation. When the molecules are stacked in a face-to-face manner, the HOMO–LUMO gap of the dimer is smaller than that of the monomer due to the interactions between the HOMOs and LUMOs of the two monomer units. When the monomer units are within a specific distance of each other, the HOMO and LUMO of the dimer, which correspond to antibonding and bonding between the units, respectively, are interchanged. This alternation of molecular orbitals may result in an increase in the bond strength between the monomer units, exhibiting stacked aromaticity. We demonstrated that it is possible to control the distance exhibited by stacked aromaticity by engineering the HOMO–LUMO gap of the monomer units

Topology Dictates Magnetic and Conductive Properties of a π‑Stacked System: Insight into Possible Coexistence of Magnetic and Conductive Systems

In this paper, conductivity and magnetism in alternant hydrocarbons are discussed based on the topology of π-conjugated networks. In a molecular system with two spin centers, when the spins are separated by an odd-length walk, they interact antiferromagnetically with each other, but when they are separated by an even-length walk, they interact ferromagnetically. The conduction through the pathway connecting the two spins is expected to be effective for the former case, while ineffective for the latter case, but both show almost the same conductance in a magnetic system. This is because in the latter case, a feature in the electron transmission spectrum that causes destructive quantum interference is localized away from the Fermi level of the electrode and in a very narrow energy range, not affecting the zero-bias conductance. This tendency is further accentuated by generating weak coupling between the electrode surfaces and the spins to preserve the radical character of the molecule sandwiched between two electrodes. Although there is a challenge on how to stabilize radical molecules in a confined environment between electrodes, what is presented in this paper would give a clue on how to construct a system where magnetism and conductivity coexist

Exploring Metal Nanocluster Catalysts for Ammonia Synthesis Using Informatics Methods: A Concerted Effort of Bayesian Optimization, Swarm Intelligence, and First-Principles Computation

This paper details the use of computational and informatics methods to design metal nanocluster catalysts for efficient ammonia synthesis. Three main problems are tackled: defining a measure of catalytic activity, choosing the best candidate from a large number of possibilities, and identifying the thermodynamically stable cluster catalyst structure. First-principles calculations, Bayesian optimization, and particle swarm optimization are used to obtain a Ti8 nanocluster as a catalyst candidate. The N2 adsorption structure on Ti8 indicates substantial activation of the N2 molecule, while the NH3 adsorption structure suggests that NH3 is likely to undergo easy desorption. The study also reveals several cluster catalyst candidates that break the general trade-off that surfaces that strongly adsorb reactants also strongly adsorb products

Molecular Understanding of Adhesion of Epoxy Resin to Graphene and Graphene Oxide Surfaces in Terms of Orbital Interactions

The adhesion mechanism of epoxy resin (ER) cured material consisting of diglycidyl ether of bisphenol A (DGEBA) and 4,4′-diaminodiphenyl sulfone (DDS) to pristine graphene and graphene oxide (GO) surfaces is investigated on the basis of first-principles density functional theory (DFT) with dispersion correction. Graphene is often used as a reinforcing filler incorporated into ER polymer matrices. The adhesion strength is significantly improved by using GO obtained by the oxidation of graphene. The interfacial interactions at the ER/graphene and ER/GO interfaces were analyzed to clarify the origin of this adhesion. The contribution of dispersion interaction to the adhesive stress at the two interfaces is almost identical. In contrast, the DFT energy contribution is found to be more significant at the ER/GO interface. Crystal orbital Hamiltonian population (COHP) analysis suggests the existence of hydrogen bonding (H-bonding) between the hydroxyl, epoxide, amine, and sulfonyl groups of the ER cured with DDS and the hydroxyl groups of the GO surface, in addition to the OH−π interaction between the benzene rings of ER and the hydroxyl groups of the GO surface. The H-bond has a large orbital interaction energy, which is found to contribute significantly to the adhesive strength at the ER/GO interface. The overall interaction at the ER/graphene is much weaker due to antibonding type interactions just below the Fermi level. This finding indicates that only dispersion interaction is significant when ER is adsorbed on the graphene surface

Exploring the Optimal Alloy for Nitrogen Activation by Combining Bayesian Optimization with Density Functional Theory Calculations

Binary alloy catalysts have the potential to exhibit higher activity than monometallic catalysts in nitrogen activation reactions. However, owing to the multiple possible combinations of metal elements constituting binary alloys, an exhaustive search for the optimal combination is difficult. In this study, we searched for the optimal binary alloy catalyst for nitrogen activation reactions using a combination of Bayesian optimization and density functional theory calculations. The optimal alloy catalyst proposed by Bayesian optimization had a surface energy of ∼0.2 eV/Å2 and resulted in a low reaction heat for the dissociation of the NN bond. We demonstrated that the search for such binary alloy catalysts using Bayesian optimization is more efficient than random search

Photo-Induced Ring-Opening Reaction of Flav-3-en-2-ol Monitored by Time-Resolved Infrared Spectroscopy

Photo-induced ring-opening reaction from flav-3-en-2-ol to 2-hydroxychalcone has been studied by time-resolved infrared (TR-IR) spectroscopy and quantum chemical calculations. A vibrational band due to the CO stretching modes for intermediate species, enol forms of 2-hydroxychalcone in the electronic ground state, was observed at 1632 cm–1 in the TR-IR spectra after photoexcitation of flav-3-en-2-ol. We also found that the CO stretching modes of the keto forms of 2-hydroxychalcone at 1664 cm–1 appeared immediately after photoexcitation and increased in intensity in synchronization with the depletion of the 1632 cm–1 band. Because the decay of the 1632 cm–1 band and the rise of the 1664 cm–1 band were fitted with bi-exponential model functions with common rate constants 0.5 and 11 μs–1, we propose that two kinds of enol form, single bond cis- (s-cis-) and trans- (s-trans-) enols, transformed into keto forms, cis-2-hydroxychalcone (Cc) and trans-2-hydroxychalcone (Ct), respectively. Quantum chemically calculated IR spectra of related species are consistent with the proposal. The observed temporal behavior of the TR-IR spectra indicates that there were reaction paths to the photogeneration of Cc and Ct within the time resolution of the TR-IR spectrometer (∼0.1 μs) in addition to the reaction paths via the enol forms of 2-hydroxychalcone