Search

When an organic molecule is placed inside a plasmonic cavity formed by two metallic nanoparticles (MNP) under illumination, the electronic excitations of the molecule couple to the plasmonic electromagnetic modes of the cavity, inducing new hybrid light-matter states called polaritons. Atomistic ab initio methods accurately describe the coupling between MNPs and molecules at the nanometer scale and allow us to analyze how atomistic features influence the interaction. In this work, we study the optical response of a porphine molecule coupled to a silver nanoparticle dimer from first principles, within the linear-response time-dependent density functional theory framework, using the recently developed Python Numeric Atomic Orbitals implementation to compute the optical excitations. The optical spectra show the splitting of the resonances of the plasmonic dimer and the molecule into two distinct polaritons, a characteristic feature of the strong light-matter coupling regime. Our results stress the importance of atomistic features, such as the gap configuration in determining the plasmon-exciton coupling strength and in the emergence of molecule-mediated charge-transfer plasmon (CTP) resonances at lower frequencies. Moreover, we show that the strength of the CTP resonance can be tuned by shifting the alignment of the molecular energy levels with respect to the Fermi level of the MNPs., (© 2024 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).)

We present a first-principles theoretical study of the atomistic footprints in the valence electron energy loss spectroscopy (EELS) of nanometer-size metallic particles. Charge density maps of excited plasmons and EEL spectra for specific electron paths through a nanoparticle (Na 380 atom cluster) are modeled using ab initio calculations within time-dependent density functional theory. Our findings unveil the atomic-scale sensitivity of EELS within this low-energy spectral range. Whereas localized surface plasmons (LSPs) are particularly sensitive to the atomistic structure of the surface probed by the electron beam, confined bulk plasmons (CBPs) reveal quantum size effects within the nanoparticle's volume. Moreover, we prove that classical local dielectric theories mimicking the atomistic structure of the nanoparticles reproduce the LSP trends observed in quantum calculations, but fall short in describing the CBP behavior observed under different electron trajectories.

Formation of carbon propeller-like molecules (CPLMs) from starphenes on a graphene substrate under electron irradiation with about 100% yield is observed in molecular dynamics simulations using the REBO-1990EVC_CH potential and CompuTEM algorithm. A CPLM consists of three carbon atomic chains connected to the central hexagon and is formed as a result of the spontaneous breaking of bonds between zigzag atomic rows in starphene arms after hydrogen removal by electron impacts. In the absence of the substrate, the CPLM yield is slightly decreased due to sticking between forming chains, while the formation time is increased threefold. The increase of the kinetic electron energy from 45 to 80 keV has no effect on the CPLM formation. Density functional theory (DFT) calculations performed show the stability of CPLMs with respect to the formation of new bonds between carbon atoms in the chains. DFT calculations using the accurate hybrid B3LYP functional provide an insight into the electronic structure of these new molecules.

Molecular doping provides a route toward designing new organic compounds with improved performance for optoelectronics. Here, we investigate the p-type doping of crystalline diindenoperylene (DIP) with two recently proposed electron-accepting molecular dopants using many-body perturbation theory. For the pristine DIP crystal, the quasiparticle band structure and the optical absorption spectra are found in agreement with the experimental data. Using the same methodology, we then characterize the optical and electronic properties of the two doped DIP crystals. The bandgap of both doped crystals is narrowed considerably due to the formation of hybridized states at the valence band edge. Moreover, a hybrid unoccupied mid-gap band is created with a host-dopant charge-transfer characteristic, giving rise to broader absorption spectra and a much lower absorption onset as compared to pristine DIP. Our results highlight that the interaction and hybridization with the host environment, including many-body effects, must be carefully considered in order to identify appropriate molecular dopants for a given organic crystal., Competing Interests: The authors declare no competing financial interest., (© 2023 American Chemical Society.)

The method of precise cutting of 2D materials by simultaneous action of a catalyst at the tip of the scanning microscope probe and an electron beam in a high-resolution transmission electron microscope is proposed and studied using atomistic simulations by the example of graphene and a nickel catalyst. Reactive molecular dynamics simulations within the Compu-TEM approach for the description of electron impact effects show that the combination of the nickel catalyst and electron irradiation is crucial for graphene cutting. Cuts with straight edges with widths of about 1-1.5 nm can be obtained. The detailed atomistic mechanism of graphene cutting is investigated via the analysis of statistics on atom ejection and bond reorganization reactions induced by the irradiation. The principal and secondary channels of atom ejection which lead to propagation of the cut are shown to be ejection of two-coordinated atoms at the cut edges bonded to the nickel tip and three-coordinated atoms from the defective graphene structure near the tip. At the same time, the ejection of two-coordinated atoms not bonded to the tip and atoms in chains at the cut edges favors smoothing of free cut edges behind the tip. A considerable difference from the atomistic mechanism of cutting a carbon nanotube via the simultaneous action of electron irradiation and nickel catalyst is discussed. The ab initio calculations performed show a decrease of the binding energy of two-coordinated carbon atoms bonded to the nickel cluster in comparison with the same cut edge in the absence of the cluster confirming that the principal channel of atom ejection is related to the cut propagation.

We present a code modularization approach to design efficient and massively parallel cubic- and linear-scaling solvers for electronic structure calculations using atomic orbitals. The modular implementation of the orbital minimization method, in which linear algebra and parallelization issues are handled via external libraries, is demonstrated in the SIESTA code. The distributed block compressed sparse row (DBCSR) and scalable linear algebra package (ScaLAPACK) libraries are used for algebraic operations with sparse and dense matrices, respectively. The MatrixSwitch and libOMM libraries, recently developed within the Electronic Structure Library, facilitate switching between different matrix formats and implement the energy minimization. We show results comparing the performance of several cubic-scaling algorithms, and also demonstrate the parallel performance of the linear-scaling solvers, and their supremacy over the cubic-scaling solvers for insulating systems with sizes of several hundreds of atoms., Competing Interests: We declare we have no competing interests., (© 2023 The Authors.)

Ab initio calculations are performed to study consecutive reconstruction of a zigzag graphene edge. According to the obtained energy profile along the reaction pathway, the first reconstruction step, formation of the first pentagon-heptagon pair, is the slowest one, while the growth of an already nucleated reconstructed edge domain should occur steadily at a much higher rate. Domains merge into one only in 1/4 of cases when they get in contact, while in the rest of the cases, residual defects are left. Structure, energy, and magnetic properties of these defects are studied. It is found that spontaneous formation of pairs of residual defects (i.e., spontaneous domain nucleation) in the fully reconstructed edge is unlikely at temperatures below 1000 K. Using a kinetic model, we show that the average domain length is several micrometers at room temperature and it decreases exponentially upon increasing the temperature at which the reconstruction takes place.

The synthesis of novel organic prototypes combining different functionalities is key to achieve operational elements for applications in organic electronics. Here we set the stage towards individually addressable magneto-optical transducers by the on-surface synthesis of optically active manganese-phthalocyanine derivatives (MnPc) obtained directly on a metallic substrate. We created these 2D nanostructures under ultra-high vacuum conditions with atomic precision starting from a simple phthalonitrile precursor with reversible photo-induced reactivity in solution. These precursors maintain their integrity after powder sublimation and coordinate with the Mn ions into tetrameric complexes and then transform into MnPcs on Ag(111) after a cyclotetramerization reaction. Using scanning tunnelling microscopy and spectroscopy together with DFT calculations, we identify the isomeric configuration of two bi-stable structures and show that it is possible to switch them reversibly by mechanical manipulation. Moreover, the robust magnetic moment brought by the central Mn ion provides a feasible pathway towards magneto-optical transducer fabrication. This work should trigger further research confirming such magneto-optical effects in MnPcs both on surfaces and in liquid environments.

Modelling the inelastic scattering of electrons in water is fundamental, given their crucial role in biological damage. In Monte Carlo track-structure (MC-TS) codes used to assess biological damage, the energy loss function (ELF), from which cross sections are extracted, is derived from different semi-empirical optical models. Only recently have first ab initio results for the ELF and cross sections in water become available. For benchmarking purpose, in this work, we present ab initio linear-response time-dependent density functional theory calculations of the ELF of liquid water. We calculated the inelastic scattering cross sections, inelastic mean free paths, and electronic stopping power and compared our results with recent calculations and experimental data showing a good agreement. In addition, we provide an in-depth analysis of the contributions of different molecular orbitals, species and orbital angular momenta to the total ELF. Moreover, we present single-differential cross sections computed for each molecular orbital channel, which should prove useful for MC-TS simulations., (© 2022 The Authors.)

Electronic circular dichroism (ECD) is a powerful spectroscopy method for investigating chiral properties at the molecular level. ECD calculations with the commonly used linear-response time-dependent density functional theory (LR-TDDFT) framework can be prohibitively costly for large systems. To alleviate this problem, we present here an ECD implementation within the projector augmented-wave method in a real-time-propagation TDDFT framework in the open-source GPAW code. Our implementation supports both local atomic basis sets and real-space finite-difference representations of wave functions. We benchmark our implementation against an existing LR-TDDFT implementation in GPAW for small chiral molecules. We then demonstrate the efficiency of our local atomic basis set implementation for a large hybrid nanocluster and discuss the chiroptical properties of the cluster.

First-principles electronic structure calculations are now accessible to a very large community of users across many disciplines, thanks to many successful software packages, some of which are described in this special issue. The traditional coding paradigm for such packages is monolithic, i.e., regardless of how modular its internal structure may be, the code is built independently from others, essentially from the compiler up, possibly with the exception of linear-algebra and message-passing libraries. This model has endured and been quite successful for decades. The successful evolution of the electronic structure methodology itself, however, has resulted in an increasing complexity and an ever longer list of features expected within all software packages, which implies a growing amount of replication between different packages, not only in the initial coding but, more importantly, every time a code needs to be re-engineered to adapt to the evolution of computer hardware architecture. The Electronic Structure Library (ESL) was initiated by CECAM (the European Centre for Atomic and Molecular Calculations) to catalyze a paradigm shift away from the monolithic model and promote modularization, with the ambition to extract common tasks from electronic structure codes and redesign them as open-source libraries available to everybody. Such libraries include "heavy-duty" ones that have the potential for a high degree of parallelization and adaptation to novel hardware within them, thereby separating the sophisticated computer science aspects of performance optimization and re-engineering from the computational science done by, e.g., physicists and chemists when implementing new ideas. We envisage that this modular paradigm will improve overall coding efficiency and enable specialists (whether they be computer scientists or computational scientists) to use their skills more effectively and will lead to a more dynamic evolution of software in the community as well as lower barriers to entry for new developers. The model comes with new challenges, though. The building and compilation of a code based on many interdependent libraries (and their versions) is a much more complex task than that of a code delivered in a single self-contained package. Here, we describe the state of the ESL, the different libraries it now contains, the short- and mid-term plans for further libraries, and the way the new challenges are faced. The ESL is a community initiative into which several pre-existing codes and their developers have contributed with their software and efforts, from which several codes are already benefiting, and which remains open to the community.

PySCF is a Python-based general-purpose electronic structure platform that supports first-principles simulations of molecules and solids as well as accelerates the development of new methodology and complex computational workflows. This paper explains the design and philosophy behind PySCF that enables it to meet these twin objectives. With several case studies, we show how users can easily implement their own methods using PySCF as a development environment. We then summarize the capabilities of PySCF for molecular and solid-state simulations. Finally, we describe the growing ecosystem of projects that use PySCF across the domains of quantum chemistry, materials science, machine learning, and quantum information science.

A review of the present status, recent enhancements, and applicability of the Siesta program is presented. Since its debut in the mid-1990s, Siesta's flexibility, efficiency, and free distribution have given advanced materials simulation capabilities to many groups worldwide. The core methodological scheme of Siesta combines finite-support pseudo-atomic orbitals as basis sets, norm-conserving pseudopotentials, and a real-space grid for the representation of charge density and potentials and the computation of their associated matrix elements. Here, we describe the more recent implementations on top of that core scheme, which include full spin-orbit interaction, non-repeated and multiple-contact ballistic electron transport, density functional theory (DFT)+U and hybrid functionals, time-dependent DFT, novel reduced-scaling solvers, density-functional perturbation theory, efficient van der Waals non-local density functionals, and enhanced molecular-dynamics options. In addition, a substantial effort has been made in enhancing interoperability and interfacing with other codes and utilities, such as wannier90 and the second-principles modeling it can be used for, an AiiDA plugin for workflow automatization, interface to Lua for steering Siesta runs, and various post-processing utilities. Siesta has also been engaged in the Electronic Structure Library effort from its inception, which has allowed the sharing of various low-level libraries, as well as data standards and support for them, particularly the PSeudopotential Markup Language definition and library for transferable pseudopotentials, and the interface to the ELectronic Structure Infrastructure library of solvers. Code sharing is made easier by the new open-source licensing model of the program. This review also presents examples of application of the capabilities of the code, as well as a view of on-going and future developments.