Ab Initio Calculations Research Articles

Oxygen-Donor Metalloligands Induce Slow Magnetization Relaxation in Zero Field for a Cobalt(II) Complex with {CoO4} Motif.

Most 3d metal-based single-molecule magnets (SMMs) use N-ligands or ligands with even softer donors to impart a particular coordination geometry and increase the zero-field splitting parameter |D|, while complexes with hard O-donor ligands showing slow magnetization relaxation are rare. Here, we report that a diamagnetic NiII complex of a tetradentate ligand featuring two N-heterocyclic carbene and two alkoxide-O donors, [LO,ONi], can serve as a {O,O'}-chelating metalloligand to give a trinuclear complex [(LO,ONi)Co(LO,ONi)](OTf)2 (2) with an elongated tetrahedral {CoIIO4} core, D = -74.3 cm-1, and a spin reversal barrier Ueff = 86.9 cm-1 in the absence of an external dc field. The influence of diamagnetic NiII on the electronic structure of the {CoO4} unit in comparison to [Co(OPh)4]2- (A) has been probed with multireference ab initio calculations. These reveal a contrapolarizing effect of the NiII, which forms stronger metal-alkoxide bonds than the central CoII, inducing a change in ligand field splitting and a 5-fold increase in the magnetic anisotropy in 2 compared to A, with an easy magnetization axis along the Ni-Co-Ni vector. This demonstrates a strategy to enhance the SMM properties of 3d metal complexes with hard O-donors by modulating the ligand field character via the coordination of diamagnetic ions and the benefit of robust metalloligands in that regard.

Complete separation of benzene-cyclohexene-cyclohexane mixtures via temperature-dependent molecular sieving by a flexible chain-like coordination polymer

The separation and purification of C6 cyclic hydrocarbons (benzene, cyclohexene, cyclohexane) represent a critically important but energy intensive process. Developing adsorptive separation technique to replace thermally driven distillation processes holds great promise to significantly reduce energy consumption. Here we report a flexible one-dimensional coordination polymer as an efficient adsorbent to discriminate ternary C6 cyclic hydrocarbons via an ideal molecular sieving mechanism. The compound undergoes fully reversible structural transformation associated with removal/re-coordination of water molecules and between activated and hydrocarbon-loaded forms. It exhibits distinct temperature- and adsorbate-dependent adsorption behavior which facilitates the complete separation of benzene, cyclohexene and cyclohexane from their binary and ternary mixtures, with the record-high uptake ratios for C6H6/C6H12 and C6H10/C6H12 in vapor phase and highest binary and ternary selectivities in liquid phase. In situ infrared spectroscopic analysis and ab initio calculations provide insight into the host-guest interactions and their effect on the preferential adsorption and structural transformation.

Photoinduced Electronic and Spin Topological Phase Transitions in Monolayer Bismuth.

Ultrathin bismuth exhibits rich physics including strong spin-orbit coupling, ferroelectricity, nontrivial topology, and light-induced structural dynamics. We use abinitio calculations to show that light can induce structural transitions to four transient phases in bismuth monolayers. These light-induced phases exhibit nontrivial topological character, which we illustrate using the recently introduced concept of spin bands and spin-resolved Wilson loops. Specifically, we find that the topology changes via the closing of the electron and spin band gaps during photoinduced structural phase transitions, leading to distinct edge states. Our study provides strategies to tailor electronic and spin topology via ultrafast control of photoexcited carriers and associated structural dynamics.

Electronic structure of CeCu4In from band structure calculations and X-ray photoelectron spectroscopy

We present a combined new experimental and theoretical study of the electronic structure for the heavy fermion CeCu4In based on X-ray photoelectron spectroscopy (XPS) data and ab-initio band structure calculations. The compound crystallizes in the orthorhombic CeCu4.38In1.62 type of structure (space group Pnnm). Below the Fermi energy the total density of states contains mainly 3d–states of Cu atoms that hybridized with the Ce 4f electronic states. The Ce core-level XPS spectra point to a stable trivalent configuration of Ce atoms in CeCu4In, consistently with the magnetic susceptibility data. For more detailed information about electronic states the fully relativistic band structure was calculated within the density functional theory (DFT) for the first time. Based on these calculations we present calculated photoemission spectra, which very well reproduce the experimental ones. The Fermi level is located at a deep decrease in the density of electronic states and reaches a value of 3.69 states/(eV f.u.). This value is much smaller than the one obtained experimentally, which indicates the importance of many-body effects, which are not properly taken into account in ab-initio calculations. The valence band is formed mainly by Cu(3d) electrons with a small contribution of Ce(4f) electrons in the immediate vicinity of the Fermi level.

Translational eigenstates of He@C60 from four-dimensional abinitio potential energy surfaces interpolated using Gaussian process regression.

We investigate the endofullerene system 3He@C60 with a four-dimensional potential energy surface (PES) to include the three He translational degrees of freedom and C60 cage radius. We compare second order Møller-Plesset perturbation theory (MP2), spin component scaled-MP2, scaled opposite spin-MP2, random phase approximation (RPA)@Perdew, Burke, and Ernzerhof (PBE), and corrected Hartree-Fock-RPA to calibrate and gain confidence in the choice of electronic structure method. Due to the high cost of these calculations, the PES is interpolated using Gaussian Process Regression (GPR), owing to its effectiveness with sparse training data. The PES is split into a two-dimensional radial surface, to which corrections are applied to achieve an overall four-dimensional surface. The nuclear Hamiltonian is diagonalized to generate the in-cage translational/vibrational eigenstates. The degeneracy of the three-dimensional harmonic oscillator energies with principal quantum number n is lifted due to the anharmonicity in the radial potential. The (2l + 1)-fold degeneracy of the angular momentum states is also weakly lifted, due to the angular dependence in the potential. We calculate the fundamental frequency to range between 96 and 110 cm-1 depending on the electronic structure method used. Error bars of the eigenstate energies were calculated from the GPR and are on the order of ∼±1.5 cm-1. Wavefunctions are also compared by considering their overlap and Hellinger distance to the one-dimensional empirical potential. As with the energies, the two abinitio methods MP2 and RPA@PBE show the best agreement. While MP2 has better agreement than RPA@PBE, due to its higher computational efficiency and comparable performance, we recommend RPA as an alternative electronic structure method of choice to MP2 for these systems.

Quantum-Chemistry Study of the Hydrolysis Reaction Profile in Borate Networks: A Benchmark.

This investigation involved an ab initio and Density Functional Theory (DFT) analysis of the hydrolysis mechanism and energetics in a borate network. The focus was on understanding how water molecules interact with and disrupt the borate network, an area where the experimental data are scarce and unreliable. The modeled system consisted of two boron atoms, bridging oxygen atoms, and varying numbers of water molecules. This setup allows for an exploration of hydrolysis under different environmental conditions, including the presence of OH- or H+ ions to simulate basic or acidic environments, respectively. Our investigation utilized both ab initio calculations at the MP2 and CCSD(T) levels and DFT with a range of exchange-correlation functionals. The findings indicate that the borate network is significantly more susceptible to hydrolysis in a basic environment, with respect to an acidic or to a neutral pH setting. The inclusion of explicit water molecules in the calculations can significantly affect the results, depending on the nature of the transition state. In fact, some transition states exhibited closed-ring configurations involving water and the boron-oxygen-boron network; in these cases, there were indeed more water molecules corresponding to lower energy barriers for the reaction, suggesting a crucial role of water in stabilizing the transition states. This study provides valuable insights into the hydrolysis process of borate networks, offering a detailed comparison between different computational approaches. The results demonstrate that the functionals B3LYP, PBE0, and wB97Xd closely approximated the reference MP2 and CCSD(T) calculated reaction pathways, both qualitatively in terms of the mechanism, and quantitatively in terms of the differences in the reaction barriers within the 0.1-0.2 eV interval for the most plausible reaction pathways. In addition, CAM-B3LYP also yielded acceptable results in all cases except for the most complicated pathway. These findings are useful for guiding further computational studies, including those employing machine learning approaches, and experimental investigations requiring accurate reference data for hydrolysis reactions in borate networks.

Discovery of natural agents against Staphylococcus aureus based on EIIC by protein modeling, virtual screening and molecular dynamics

In this study, a reasonable 3D structure of enzyme IIC (EIIC) was established using ab initio methods. Curcumin, baicalein and phloretin were screened out from 127,695 natural products based on EIIC protein in Staphylococcus aureus (S. aureus) by using the virtual screening technology. The antibacterial assays results confirmed that these three compounds could inhibit S. aureus cell growth. The results of molecular docking, molecular dynamics and molecular mechanics/generalized born surface area calculation indicated that curcumin, baicalein and phloretin mainly bound to Glu275, His231, and Phe359 of EIIC through hydrogen bonds and van der Waals forces to obstruct the transport of glucose-based carbohydrates, and thereby inhibited the carbohydrate metabolism of S. aureus. The steered molecular dynamics analysis revealed the existence of transfer channel in EIIC and illustrated that curcumin, baicalein and phloretin could dissociate from EIIC when pulling force on them. Furthermore, alanine scanning and hydrogen bond occupancy analyses found that Glu275 played an important role in stabilizing EIIC-baicalein/curcumin/phloretin complexes, and 275Glu@OE1 was a stable hydrogen bond acceptor. Therefore, Glu275 may be a hot spot residue that can maintain the stable binding of agents to EIIC, which may be a direction for the development of new drugs against S. aureus.

Quantum tomography of molecules using ultrafast electron diffraction.

We propose a quantum tomography (QT) approach to retrieve the temporally evolving reduced density matrix in electronic state basis, where the populations and coherence between the ground state and excited state are reconstructed from the ultrafast electron diffraction signal. In order to showcase the capability of the proposed QT approach, we simulate the nuclear wavepacket dynamics and ultrafast electron diffraction of photoexcited pyrrole molecules using the abinitio quantum chemical CASSCF method. From the simulated time-resolved diffraction data, we retrieve the evolving density matrix in a crude diabatic representation basis and reveal the symmetry of the excited pyrrole wavepacket. Our QT approach opens the route to make a quantum version of "molecular movie" that covers the electronic degree of freedom and equips ultrafast electron diffraction with the power to reveal the coherence between electronic states, relaxation, and dynamics of population transfer.

Dipolar Coupling as a Mechanism for Fine Control of Magnetic States in ErCOT-Alkyl Molecular Magnets.

The design of molecular magnets has progressed greatly by taking advantage of the ability to impart successive perturbations and control vibronic transitions in 4fn systems through the careful manipulation of the crystal field. Herein, we control the orientation and rigidity of two dinuclear ErCOT-based molecular magnets: the inversion-symmetric bridged [ErCOT(μ-Me)(THF)]2 (2) and the nearly linear Li[(ErCOT)2(μ-Me)3] (3). The conserved anisotropy of the ErCOT synthetic unit facilitates the direction of the arrangement of its magnetic anisotropy for the purposes of generating controlled internal magnetic fields, improving control of the energetics and transition probabilities of the electronic angular momentum states with exchange biasing via dipolar coupling. This control is evidenced through the introduction of a second thermal barrier to relaxation operant at low temperatures that is twice as large in 3 as in 2. This barrier acts to suppress through-barrier relaxation by protecting the ground state from interacting with stray local fields while operating at an energy scale an order of magnitude smaller than the crystal field term. These properties are highlighted when contrasted against the mononuclear structure ErCOT(Bn)(THF)2 (1), in which quantum tunneling of the magnetization processes dominate, as demonstrated by magnetometry and ab initio computational methods. Furthermore, far-infrared magnetospectroscopy measurements reveal that the increased rigidity imparted by successive removal of solvent ligands when adding bridging methyl groups, along with the increased excited state purity, severely limits local spin-vibrational interactions that facilitate magnetic relaxation, manifesting as longer relaxation times in 3 relative to those in 2 as temperature is increased.

The Quenching Mechanism and Suppressing Method of Single-Photon Emitters from Point Defects in GaN.

The defect-based single-photon emitters (SPEs) in gallium nitride (GaN) have attracted considerable research interest due to their high emission rate, narrow line width, and room-temperature operation. However, the quenching effect greatly restricts the applications of these SPEs, and the origin of the quenching mechanism is still unclear. Here, based on systematic ab initio calculations, we reveal a possible quenching mechanism originating from the transformation between two different structures of the defect-pair NGaVN in wurtzite GaN. Our results indicate that the defect-pair NGaVN possesses two stable detect-structures A and B, where the structure B has a small zero phonon line (ZPL) and long lifetime. The transformation barrier from structures A to B is only 0.097 eV. Thus, structure A can easily transform to structure B under laser illumination due to thermal fluctuations, causing a quenching phenomenon. Our work also predicts that the barrier energy between defect structures A and B could be effectively adjusted through tuning the triaxial compressive strain of the crystal structure. This provides an effective method to suppress the quenching effect of defect-pair NGaVN in GaN, paving the way for practical applications of SPEs.

ElasTool v3.0: Efficient computational and visualization toolkit for elastic and mechanical properties of materials

Efficient computation and visualization of elastic and mechanical properties are crucial in the selection of materials and the design of new materials. The ElasTool v3.0 toolkit marks a significant advancement in the computational analysis and visualization of elastic and mechanical properties of materials, essential in material selection and design. This enhanced version extends beyond standard calculations like elastic tensor, Young's modulus, bulk modulus, and Poisson's ratio. It introduces capabilities for computing minimum thermal conductivity, linear compressibility, rendering the Christoffel equation, and elastic energy density. Notably, it integrates advanced visualization tools, including compatibility with Plotly and Elate web platforms for interactive web-based property exploration. A key feature of ElasTool v3.0 is the implementation of second-order elastic constants (SOECs) for tubular 2D-based nanostructures and nanotubes. Leveraging high-efficiency strain-matrix sets (OHESS), the toolkit now facilitates efficient computation of elastic constants and mechanical properties at both zero and finite temperatures for 1D, 2D, and 3D dimensions. ElasTool is openly accessible on GitHub: https://github.com/elastool. Program summaryProgram Title: ElasTool v3.0CPC Library link to program files:https://doi.org/10.17632/9hpgmspfmk.1Licensing provisions: GNU General Public License, version 3Programming language: Python 3Nature of problem: A computational toolkit that can efficiently compute the elastic and mechanical properties of materials using ab initio methods and provide advanced visualization capabilities is crucial to accurate material analysis.Solution method:ElasTool addresses this need by incorporating high-efficiency strain-matrix sets that require only small strain sets with ab initio methods such as density functional theory and molecular dynamics. This package is versatile and adept at calculating elastic constants for both 2D and 3D structures, including tubular 2D-based nanostructures and nanotubes and 1D nanoribbons at zero and finite temperatures. Beyond basic elastic moduli, ElasTool now boasts advanced visualization capabilities and computes additional critical properties such as minimum thermal conductivity, elastic energy density, and linear compressibility. This enhancement transforms ElasTool into a comprehensive tool for accelerated material design and modeling, and advanced visualization, offering both efficiency and depth in material property analysis.Additional comments including restrictions and unusual features:ElasTool is currently integrated with the Vienna Ab initio Simulation Package (VASP). Adapting to other electronic structures is straightforward.ElasTool offers serial and parallel execution. Depending on the system size, the number of strains, the nature of the calculations (zero or finite temperature), and the computing architecture, the computing time can vary from a minute to as long as 15 hours.

On the performance of composite schemes in determining equilibrium molecular structures.

Determination of equilibrium molecular structures is an essential ingredient in predicting spectroscopic parameters that help in identifying molecular carriers of microwave transitions. Here, the performance of two different ab initio composite approaches for obtaining equilibrium structures, "energy scheme" and "geometry scheme," is explored and compared to semi-experimental equilibrium structures. This study is performed for a set of 11 molecules which includes diatomics, linear triatomics, and a few non-linear molecules. The ab initio calculations were performed using three tiers of composite chemical recipes. The current results show that as the overall rigor of calculation is increased, the semi-experimental and the ab initio numbers agree to within 0.0003 Å for all molecules in the test set. The composite approach based on correcting the potential energy surface (energy scheme) and the one based on correcting the geometry directly (geometry scheme) show excellent agreement with each other. This work represents a step toward development of efficient and highly accurate procedures for computing ab initio equilibrium structures.

Reducing anisotropy of rhombohedral Bi-rich phase for high-performance Ag-alloyed Sn–Bi low-temperature solders

The rhombohedral structure is with anisotropic properties, leading to brittle mechanical properties. However, it is one of the major constituent phases, the (Bi) phase, in the Sn–Bi-based low-temperature solders, which is demanding in the electronic industry for reducing the soldering temperature and alleviating the warpage in high-density fine-pitch electronic packaging. Herein, we proposed a high-performance Ag-alloyed Sn–Bi (SBA) low-temperature solder by reducing the fraction and, more importantly, the anisotropy of the (Bi) phase. We found that the (Bi) phase exhibited considerable solubility for Ag and Sn, even after thermal aging, without any nanoprecipitates as confirmed with both CALPHAD-type thermodynamic modeling and high-resolution TEM. The solution effects of Ag and Sn in the (Bi) phase is elaborated based on in situ nano-indentation, EPMA compositional analyses, EBSD characterizations, and ab initio calculations. The presence of Ag and Sn in the (Bi) phases resulted in both solid solution hardening and softening effects, depending on the Schmid factor of the grain. As a result, in (Bi) polycrystals, alloying could effectively mitigate the anisotropic mechanical properties of the (Bi) phase and facilitated a more uniform dispersion of forces exerted among (Bi) grains in all orientations. Consequently, the designed SBA solder exhibited superior mechanical properties, e.g., an ultimate tensile strength (UTS) of 55 MPa, elongation of 46%, and toughness of 20.6 × 106 J/m3 at a strain rate of 0.0005 s−1. This understanding opens the door for enhanced rhombohedral phase in alloy design for various applications.

Strain-induced magnetic anisotropy of multi-domain epitaxial EuPd2 thin films

Europium intermetallic compounds show a variety of different ground states and anomalous physical properties due to the interactions between the localized 4f electrons and the delocalized electronic states. Europium is also the most reactive of the rare earth metals which might be the reason why very few works are concerned with the properties of Eu-based thin films. Here we address the low-temperature magnetic properties of ferromagnetic EuPd2 thin films prepared by molecular beam epitaxy. The epitaxial (111)-oriented thin films grow on MgO (100) with eight different domain orientations. We analyze the low-temperature magnetic hysteresis behavior by means of micromagnetic simulations taking the multi-domain morphology explicitly into account and quantify the magnetic crystal anisotropy contribution. By ab initio calculations we trace back the microscopic origin of the magnetic anisotropy to thin film-induced uniform biaxial strain.

Tracking electron motion within and outside of Floquet bands from attosecond pulse trains in time-resolved ARPES

Floquet engineering has recently emerged as a technique for controlling material properties with light. Floquet phases can be probed with time- and angle-resolved photoelectron spectroscopy (Tr-ARPES), providing direct access to the laser-dressed electronic bands. Applications of Tr-ARPES to date focused on observing the Floquet-Bloch bands themselves, and their build-up and dephasing on sub-laser-cycle timescales. However, momentum and energy resolved sub-laser-cycle dynamics between Floquet bands have not been analyzed. Given that Floquet theory strictly applies in time-periodic conditions, the notion of resolving sub-laser-cycle dynamics between Floquet states seems contradictory—it requires probe pulse durations below a laser cycle that inherently cannot discern the time-periodic nature of the light-matter system. Here we propose to employ attosecond pulse train probes with the same temporal periodicity as the Floquet-dressing pump pulse, allowing both attosecond sub-laser-cycle resolution and a proper projection of Tr-ARPES spectra on the Floquet–Bloch bands. We formulate and employ this approach in ab-initio calculations in light-driven graphene. Our calculations predict significant sub-laser-cycle dynamics occurring within the Floquet phase with the majority of electrons moving within and in-between Floquet bands, and a small portion residing and moving outside of them in what we denote as ‘non-Floquet’ bands. We establish that non-Floquet bands arise from the pump laser envelope that induces non-adiabatic electronic excitations during the pulse turn-on and turn-off. By performing calculations in systems with poly-chromatic pumps we also show that Floquet states are not formed on a sub-laser-cycle level. This work indicates that the Floquet-Bloch states are generally not a complete basis set for sub-laser-cycle dynamics in steady-state phases of matter.

The Effective Fragment Molecular Orbital Method: Achieving High Scalability and Accuracy for Large Systems.

The effective fragment molecular orbital (EFMO) method has been developed to predict the total energy of a very large molecular system accurately (with respect to the underlying quantum mechanical method) and efficiently by taking advantage of the locality of strong chemical interactions and employing a two-level hierarchical parallelism. The accuracy of the EFMO method is partly attributed to the accurate and robust intermolecular interaction prediction between distant fragments, in particular, the many-body polarization and dispersion effects, which require the generation of static and dynamic polarizability tensors by solving the coupled perturbed Hartree-Fock (CPHF) and time-dependent HF (TDHF) equations, respectively. Solving the CPHF and TDHF equations is the main EFMO computational bottleneck due to the inefficient (serial) and I/O-intensive implementation of the CPHF and TDHF solvers. In this work, the efficiency and scalability of the EFMO method are significantly improved with a new CPU memory-based implementation for solving the CPHF and TDHF equations that are parallelized by either message passing interface (MPI) or hybrid MPI/OpenMP. The accuracy of the EFMO method is demonstrated for both covalently bonded systems and noncovalently bound molecular clusters by systematically examining the effects of basis sets and a key distance-related cutoff parameter, Rcut. Rcut determines whether a fragment pair (dimer) is treated by the chosen ab initio method or calculated using the effective fragment potential (EFP) method (separated dimers). Decreasing the value of Rcut increases the number of separated (EFP) dimers, thereby decreasing the computational effort. It is demonstrated that excellent accuracy (<1 kcal/mol error per fragment) can be achieved when using a sufficiently large basis set with diffuse functions coupled with a small Rcut value. With the new parallel implementation, the total EFMO wall time is substantially reduced, especially with a high number of MPI ranks. Given a sufficient workload, nearly ideal strong scaling is achieved for the CPHF and TDHF parts of the calculation. For the first time, EFMO calculations with the inclusion of long-range polarization and dispersion interactions on a hydrated mesoporous silica nanoparticle with explicit water solvent molecules (more than 15k atoms) are achieved on a massively parallel supercomputer using nearly 1000 physical nodes. In addition, EFMO calculations on the carbinolamine formation step of an amine-catalyzed aldol reaction at the nanoscale with explicit solvent effects are presented.

Wandering principal optical axes in van der Waals triclinic materials

Nature is abundant in material platforms with anisotropic permittivities arising from symmetry reduction that feature a variety of extraordinary optical effects. Principal optical axes are essential characteristics for these effects that define light-matter interaction. Their orientation – an orthogonal Cartesian basis that diagonalizes the permittivity tensor, is often assumed stationary. Here, we show that the low-symmetry triclinic crystalline structure of van der Waals rhenium disulfide and rhenium diselenide is characterized by wandering principal optical axes in the space-wavelength domain with above π/2 degree of rotation for in-plane components. In turn, this leads to wavelength-switchable propagation directions of their waveguide modes. The physical origin of wandering principal optical axes is explained using a multi-exciton phenomenological model and ab initio calculations. We envision that the wandering principal optical axes of the investigated low-symmetry triclinic van der Waals crystals offer a platform for unexplored anisotropic phenomena and nanophotonic applications.

Transient Photoinduced Pb2+ Disproportionation for Exciton Self-Trapping and Broadband Emission in Low-Dimensional Lead Halide Perovskites.

Low-dimensional lead halide perovskites with broadband emission hold great promise for single-component white-light-emitting (WLE) devices. The origin of their broadband emission has been commonly attributed to self-trapped excitons (STEs) composed of localized electronic polarization with a distorted lattice. Unfortunately, the exact electronic and structural nature of the STE species in these WLE materials remains elusive, hindering the rational design of high-efficiency WLE materials. In this study, by combining ultrafast transient absorption spectroscopy and ab initio calculations, we uncover surprisingly similar STE features in two prototypical low dimensional WLE perovskite single crystals: 1D (DMEDA)PbBr4 and 2D (EDBE)PbBr4, despite of their different dimensionalities. Photoexcited excitons rapidly localize to intrinsic STEs within ∼250 fs, contributing to the white light emission. Crucially, STEs in both systems exhibit characteristic absorption features akin to those of Pb+ and Pb3+. Further atomic level theoretical simulations confirm photoexcited electrons and holes are localized on the Pb2+ site to form Pb+- and Pb3+-like species, resembling transient photoinduced Pb2+ disproportionation. This study provides conclusive evidence on the key excited state species for exciton self-trapping and broadband emission in low dimensional lead halide WLE perovskites and paves the way for the rational design of high-efficiency WLE materials.

Tuning (Anti)Aromaticity: Variations on the [8]-Circulene Framework.

Optoelectronic properties of organic molecules are underpinned by delocalisation and delocalisability of π-electrons. These properties are sensitive to small changes in electron count, whether achieved by heteroatom substitution or redox chemistry. One measure of the delocalisability of π-electrons is the current induced by an external magnetic field, which is diagnostic of (anti)aromaticity. The ab initio ipsocentric method is used here to model diverse ring-current patterns in the family of [8]-circulenes based on tetracyclopenta[def,jkl,pqr,vwx]tetraphenylene (TCPTP), in different charge states, with disjoint hetero-atom substitution, and with CC units systematically replaced by BN pairs. Maps calculated at the CHF/CTOCD-DZ2/6-31G** level reveal that these modifications of the TCPTP framework access the full range of possibilities for current from concentric global circulations (typically counter rotating) to full (non-aromatic) localisation. In the ipsocentric approach, induced current density is partitioned into robust orbital contributions that obey selection rules based on orbital symmetry, energy and nodal character. The selection rules are applied here to interpret current-density and exploit insights gained from simpler models to suggest design strategies for fine-tuning of π-delocalisability (aromaticity and antiaromaticity) in macrocyclic frameworks.

From proton transfer to ionic liquids: How isolated domains of ion pairs grow to form extended network in diphenyl phosphate-bis(2-ethylhexyl) amine mixtures

HypothesisDiphenyl phosphate (DPP)/bis(2-ethylhexyl) amine (BEEA) liquid mixtures are expected to show peculiar self-assembly behavior due to their, respectively, acidic and basic nature which can trigger acid-base reaction. ExperimentThe properties of DPP/BEEA mixtures of different compositions are explored by a combined theoretical (Ab initio calculations) and experimental (rheology, X-ray and light scattering) approach. Findings(i) a proton transfer from DPP to BEEA takes place with formation of cationic and anionic species in liquid phase (ion pair formation); (ii) due to the strong and long-range electrostatic interactions, each of the two charged species is preferentially surrounded by the other one, in a picture resembling the structure of ionic materials; (iii) this gives rise to a striking viscosity increase takes place as DPP concentration increases; (iv) the composition dependencies of all measured parameters show clear deviations around DPP molar fraction of 0.2, indicating unequivocally the formation of isolated domains, formed by DPP-BEEA pairs which coalescence/percolate into a unique, extended, network at higher DPP concentrations.