Abstract: A theoretical and experimental study of the work function of few-layer graphene is reported. The influence of the number of layers on the work function is investigated in the presence of a substrate, a molecular dipole layer, and combinations of the two. The work function of few-layer graphene is almost independent of the number of layers with only a difference between monolayer and multilayer graphene of about 60 meV. In the presence of a charge-donating substrate the charge distribution is found to decay exponentially away from the substrate and this is directly reflected in the work function of few-layer graphene. A dipole layer changes the work function only when placed in between the substrate and few-layer graphene through a change of the charge transfer between the two.

Abstract: Polarons are quasiparticles emerging in materials from the interaction of extra charge carriers with the surrounding atomic lattice. They appear in a wide va- riety of compounds and can have a profound impact on their properties, making the concept of a polaron a central and ubiquitous topic in material science. Al- though the concept is known for about 75 years, the origin of polarons is not yet fully elucidated. This thesis focuses on WO 3 as a well-known prototypical system for studying polarons, which inherent polaronic nature is linked to its remark- able electrical and chromic properties. The primary objective of this research is to provide a comprehensive atomistic description and understanding of polaron formation in WO 3 using first-principles density functional theory (DFT) calcula- tions. Additionally, the investigation explores the interactions between polarons and the possibility of bipolaron formation. Following a systematic strategy, we first extensively analyze the dielectric and lattice dynamical properties of WO 3 in both the room-temperature P 2 1 /n and ground-state P 2 1 /c phases. Our specific focus is on characterizing the zone-center phonons, which serve as the founda- tion for identifying the phonon modes involved in the polaron formation and charge localization process. Subsequently, we examine the impact of structural distortions on the electronic structure of WO 3 to elucidate the interplay between structural distortions and electronic properties, thereby laying the groundwork for understanding electron-phonon couplings. By incorporating these critical fac- tors, we address our primary research goals. The most common explanation for the polaron formation is associated with the electrostatic screening of the extra charge by the polarizable lattice. Here, we show that, even in ionic crystals, this is not necessarily the case. We demonstrate that polarons in this compound arise primarily from non-polar atomic distortions. We then unveil that this unexpected behavior originates from the undoing of distortive atomic motions, which lowers the bandgap. As such, we coin the name of anti-distortive polaron and validate its appearance through a simple quantum-dot model, in which charge localization is the result of balancing structural, electronic, and confinement energy costs. Then, we also study the polaron-polaron interaction and present the formation of the antiferromagnetic W 4+ bipolaronic state with relatively large formation energy. Our analysis of the W 4+ bipolaronic distortions on the global structure reveals the same behavior as in experiments where the highly distorted monoclinic phase transforms into a tetragonal phase as a function of doping. Additionally, leveraging our previous findings on asymmetric polaronic distortion and examin- ing different merging orientations, we stabilize the antiferromagnetic W 5+ -W 5+ bipolaronic state with an energy lower than the W 4+ state. This thesis clari- fies the formation of unusual medium-size 2D polarons and bipolarons in WO3,which might be relevant to the whole family of ABO 3 perovskites, to which WO 3 is closely related. The simplicity of the concept provides also obvious guidelines for tracking similar behavior in other families of compounds.

Abstract: Cadmium sulfide (CdS)-based yellow pigments have been used in a number of early 20th century artworks, including The Scream series painted by Edvard Munch. Some of these unique paintings are threatened by the discoloration of these CdS-based yellow oil paints because of the oxidation of the original sulfides to sulfates. The experimental data obtained here prove that moisture and cadmium chloride compounds play a key role in promoting such oxidation. To clarify how these two factors effectively prompt the process, we studied the band alignment between CdS, CdCl2, and Cd-(OH)Cl as well as the radicals center dot OH and H3O center dot by density functional theory (DFT) methods. Our results show that a stack of several layers of Cd-(OH)Cl creates a pocket of positive holes at the Cl-terminated surface and a pocket of electrons at the OH-terminated surface by leading in a difference in ionization energy at both surfaces. The resulting band alignment indicates that Cd-(OH)Cl can indeed play the role of an oxidative catalyst for CdS in a moist environment, thus providing an explanation for the experimental evidence.

Abstract: The influence of decoration with impurities and the composition dependent band gap in 2D materials has been the subject of debate for a long time. Here, by using Density Functional Theory (DFT) calculations, we systematically disclose physical properties of chlorinated phosphorene having the stoichiometry of PmCln. By analyzing the adsorption energy, charge density, migration energy barrier, structural, vibrational, and electronic properties of chlorinated phosphorene, we found that (I) the Cl-P bonds are strong with binding energy Eb =-1.61 eV, decreases with increasing n. (II) Cl atoms on phosphorene have anionic feature, (III) the migration path of Cl on phosphorene is anisotropic with an energy barrier of 0.38 eV, (IV) the phonon band dispersion reveal that chlorinated phosphorenes are stable when r <= 0.25 where r = m/n, (V) chlorinated phosphorenes is found to be a photonic crystal in the frequency range of 280 cm-1 to 325 cm-1, (VI) electronic band structure of chlorinated phosphorenes exhibits quasi-flat bands emerging around the Fermi level with widths in the range of 22 meV to 580 meV, and (VII) Cl adsorption causes a semiconducting to metallic/semi-metallic transition which makes it suitable for application as an electroactive material. To elucidate this application, we investigated the change in binding energy (Eb), specific capacity, and open-circuit voltage as a function of the density of adsorbed Cl. The theoretical storage capacity of the chlorinated phosphorene is found to be 168.19 mA h g-1with a large average voltage (similar to 2.08 V) which is ideal number as a cathode in chloride-ion batteries.

Abstract: Electrocatalysts play a pivotal role in advancing the application of water splitting for hydrogen production. This research unveils the potential of defective biphenylenes as high-efficiency catalysts for the hydrogen evolution reaction. Using first-principles simulations, we systematically investigated the structure, stability, and catalytic performance of defective biphenylenes. Our findings unveil that defect engineering significantly enhances the electrocatalytic activity for hydrogen evolution. Specifically, biphenylene with a double-vacancy defect exhibits an outstanding Gibbs free energy of -0.08 eV, surpassing that of Pt, accompanied by a remarkable exchange current density of -3.08 A cm(-2), also surpassing that of Pt. Furthermore, we find the preference for the Volmer-Heyrovsky mechanism in the hydrogen evolution reaction, with a low energy barrier of 0.80 eV. This research provides a promising avenue for developing novel metal-free electrocatalysts for water splitting with earth-abundant carbon elements, making a significant step toward sustainable hydrogen production.

Abstract: Lattice thermal conductivity (kappa(L)) is of great scientific interest for the development of efficient energy conversion technologies. Therefore, microscopic understanding of phonon transport is critically important for designing functional materials. In our previous study (Roshan et al., ACS Applied Energy Mater. 2021, 5, 882-896), anomalous kappa(L) trends were predicted for rocksalt alkaline-earth chalcogenides (AECs). In the present work, we extended it to alkali halides (AHs) and conducted a thorough investigation to explore the role of atomic mass contrast on lattice dynamics and phonon transport properties of 36 binary compounds (20 AHs + 16 AECs). The calculated spectral and cumulative kappa(L) reveal that low-lying optical phonon modes significantly boost kappa(L) alongside acoustic phonons in materials where the atomic mass ratio approaches unity and cophonocity nears zero. Phonon scattering rates are relatively low for materials with a mass ratio close to one, and the corresponding phonon lifetimes are higher, which enhances kappa(L). Phonon lifetimes play a critical role, outweighing phonon group velocities, in determining the anomalous trends in kappa(L) for both AHs and AECs. To further explore the role of atomic mass contrast in kappa(L), the effect of tensile lattice strain on phonon transport has also been investigated. Under tensile strain, both group velocities and phonon lifetimes decrease in the low frequency range, leading to a decrease in kappa(L). This work provides insights on how atomic mass contrast can tune the contribution of optical phonons to kappa(L) and its implications on scattering rates by either enhancing or suppressing kappa(L). These insights would aid in the selection of elements for designing new functional materials with and without atomic mass contrast to achieve relatively high and low kappa(L) values, respectively.

Abstract: The atomically laminated layered ternary transition-metal borides (the MAB phases) have demonstrated outstanding properties and have been applied in various fields. Understanding their thermal and mechanical properties is critical to determining their applicability in various fields such as high-temperature applications. To achieve this, we conducted first-principles calculations based on density-functional theory and the quasiharmonic approximation to determine the thermal expansion coefficients, Gruneisen parameters, bulk moduli, hardness, thermal conductivity, electron-phonon coupling parameters, and the structural and vibrational properties of MoAlB, WAlB, Cr2AlB2, and Ti2InB2. We found varying degrees of anisotropy in the thermal expansion and mechanical properties in spite of similarities in their crystal structures. MoAlB has a mild degree of anisotropy in its thermal expansion coefficient (TEC), while Cr2AlB2 and WAlB display the highest level of TEC anisotropy. We assessed various empirical models to calculate hardness and thermal conductivity, and correlated the calculated values with the material properties such as elastic moduli, Gruneisen parameter, Debye temperature, and type of bonding. Owing to their higher Gruneisen parameters, implying a greater degree of anharmonicity in lattice vibrations and lower phonon group velocities, MoAlB and WAlB have significantly lower lattice thermal conductivity values than those of Cr2AlB2 and Ti2InB2. The hardness and lattice thermal conductivity of MAB phases can be predicted with high accuracy if one utilizes an appropriate model.

Abstract: Hexagonal boron nitride (hBN) is considered as an ideal semiconductor material for solid-state neutron detector, owing to its large neutron scattering section because of the low atomic number of B and excellent physical properties. Here we study the influence of neutron irradiation on crystal structure and on intermediate energy state (IMES) levels induced by the presence of impurities and defects in hBN. Large-size and thick pyrolytic and hot-pressed hBN (PBN and HBN) samples, which can be directly applied for neutron detector devices, are prepared and bombarded by neutrons with different irradiation fluences. The SEM and TEM are used to observe the sample difference of PBN and HBN. X-ray diffraction and Raman spectroscopy are applied to examine the influence of neutron irradiation on lattice structures along different crystal directions of PBN and HBN samples. Photoluminescence (PL) is employed to study the effect of neutron irradiation on IMESs in these samples. We find that the neutron irradiation does not alter the in-plane lattice structures of both PBN and HBN samples, but it can release the inter-layer tensions induced by sample growth of the PBN samples. Interestingly and surprisingly, the neutron irradiation does not affect the IMES levels responsible for PL generation, where PL is attributed mainly from phonon-assisted radiative electron-hole coupling for both PBN and HBN samples. Furthermore, the results indicate that the neutron irradiation can weaken the effective carrier-phonon coupling and exciton transitions in PBN and HBN samples. Overall, both PBN and HBN samples show some degree of the resistance to neutron irradiation in terms of these basic physical properties. The interesting and important findings from this work can help us to gain an in-depth understanding of the influence of neutron irradiation on basic physical properties of hBN materials. These effects can be taken into account when designing and applying the hBN materials for neutron detectors.

Abstract: Recently, a process has been proposed for generating negatively-curved carbon schwarzites via zeolite-templating (Braun et al., 2018). However, the proposed process leads to atomistic models which are not very symmetric and often rather defective. In the present work, an improved generation approach is developed, by imposing symmetry constraints, which systematically leads to defect-free, hence more stable, schwarzites. The stability of the newly predicted symmetric schwarzites is also compared to that of other carbon nanostructures (in particular carbon nanotubes – CNTs), which could also be accommodated within the same templates. Our results suggest that only a few of these (such as FAU, SBT and SBS) can fit schwarzites more stable than CNTs. Our predictions could help experimentalists in the crucial choice of the template for the challenging synthesis of schwarzites. Furthermore, being highly symmetric and stable phases, the models could also be synthesized by means of other experimental procedures.

Abstract: We investigate the coupling factor phi( mu) that quantifies the magnetic flux phi per magnetic moment mu of a point-like magnetic dipole that couples to a superconducting quantum interference device (SQUID). Representing the dipole by a tiny current-carrying (Amperian) loop, the reciprocity of mutual inductances of SQUID and Amperian loop provides an elegant way of calculating phi(mu)(r,e(mu)) vs. position r and orientation e(mu) of the dipole anywhere in space from the magnetic field B-J(r) produced by a supercurrent circulating in the SQUID loop. We use numerical simulations based on London and Ginzburg-Landau theory to calculate phi (mu) from the supercurrent density distributions in various superconducting loop geometries. We treat the far-field regime ( r greater than or similar to a= inner size of the SQUID loop) with the dipole placed on (oriented along) the symmetry axis of circular or square shaped loops. We compare expressions for phi (mu) from simple filamentary loop models with simulation results for loops with finite width w (outer size A > alpha), thickness d and London penetration depth lambda(L )and show that for thin ( d << alpha ) and narrow (w < alpha) loops the introduction of an effective loop size a(eff) in the filamentary loop-model expressions results in good agreement with simulations. For a dipole placed right in the center of the loop, simulations provide an expression phi(mu)(a,A,d,lambda(L)) that covers a wide parameter range. In the near-field regime (dipole centered at small distance z above one SQUID arm) only coupling to a single strip representing the SQUID arm has to be considered. For this case, we compare simulations with an analytical expression derived for a homogeneous current density distribution, which yields excellent agreement for lambda(L)>w,d . Moreover, we analyze the improvement of phi(mu) provided by the introduction of a narrow constriction in the SQUID arm below the magnetic dipole.

Abstract: Knowledge of how heat flows anisotropically in van der Waals (vdW) materials is crucial for thermal management of emerging 2D materials devices and design of novel anisotropic thermoelectric materials. Despite the importance, anisotropic heat transport in vdW materials is yet to be systematically studied and is often presumably attributed to anisotropic speeds of sound in vdW materials due to soft interlayer bonding relative to covalent in-plane networks of atoms. In this work, we investigate the origins of the anisotropic heat transport in vdW materials, through time-domain thermoreflectance (TDTR) measurements and first-principles calculations of anisotropic thermal conductivity of three different phases of MoTe2. MoTe2 is ideal for the study due to its weak anisotropy in the speeds of sound. We find that even when the speeds of sound are roughly isotropic, the measured thermal conductivity of MoTe2 along the c-axis is 5-8 times lower than that along the in-plane axes. We derive meaningful characteristic heat capacity, phonon group velocity, and relaxation times from our first principles calculations for selected vdW materials (MoTe2, BP, h-BN, and MoS2), to assess the contributions of these factors to the anisotropic heat transport. Interestingly, we find that the main contributor to the heat transport anisotropy in vdW materials is anisotropy in heat capacity of the dominant heat-carrying phonon modes in different directions, which originates from anisotropic optical phonon dispersion and disparity in the frequency of heat-carrying phonons in different directions. The discrepancy in frequency of the heat-carrying phonons also leads to similar to 2 times larger average relaxation times in the cross-plane direction, and partially explains the apparent dependence of the anisotropic heat transport on the anisotropic speeds of sound. This work provides insight into understanding of the anisotropic heat transport in vdW materials.

Abstract: Defect-free graphene is impermeable to all atoms(1-5) and ions(6,7) under ambient conditions. Experiments that can resolve gas flows of a few atoms per hour through micrometre-sized membranes found that monocrystalline graphene is completely impermeable to helium, the smallest atom(2,5). Such membranes were also shown to be impermeable to all ions, including the smallest one, lithium(6,7). By contrast, graphene was reported to be highly permeable to protons, nuclei of hydrogen atoms(8,9). There is no consensus, however, either on the mechanism behind the unexpectedly high proton permeability(10-14) or even on whether it requires defects in graphene's crystal lattice(6,8,15-17). Here, using high-resolution scanning electrochemical cell microscopy, we show that, although proton permeation through mechanically exfoliated monolayers of graphene and hexagonal boron nitride cannot be attributed to any structural defects, nanoscale non-flatness of two-dimensional membranes greatly facilitates proton transport. The spatial distribution of proton currents visualized by scanning electrochemical cell microscopy reveals marked inhomogeneities that are strongly correlated with nanoscale wrinkles and other features where strain is accumulated. Our results highlight nanoscale morphology as an important parameter enabling proton transport through two-dimensional crystals, mostly considered and modelled as flat, and indicate that strain and curvature can be used as additional degrees of freedom to control the proton permeability of two-dimensional materials. A study using high-resolution scanning electrochemical cell microscopy attributes proton permeation through defect-free graphene and hexagonal boron nitride to transport across areas of the structure that are under strain.

Abstract: Hydrogen production through photocatalytic water splitting is being developed swiftly to address the ongoing energy crisis. Over the past decade, with the rise of graphene and other two-dimensional (2D) materials, an increasing number of computational and experimental studies have focused on relevant van der Waals (vdW) semiconductor heterostructures for photocatalytic water splitting. In this review, the fundamental mechanism and distinctive performance of type-II and Z-scheme vdW heterostructure photocatalysts are presented. Accordingly, we have conducted a systematic review of recent studies focusing on candidates for photocatalysts, specifically vdW heterostructures involving 2D transition metal disulfides (TMDs), 2D Janus TMDs, and phosphorenes. The photocatalytic performance of these heterostructures and their suitability in theoretical scenarios are discussed based on their electronic and optoelectronic properties, particularly in terms of band structures, photoexcited carrier dynamics, and light absorption. In addition, various approaches for tuning the performance of these potential photocatalysts are illustrated. This strategic framework for constructing and modulating 2D heterostructure photocatalysts is expected to provide inspiration for addressing possible challenges in future studies.

Abstract: The EPS Antwerp Young Minds (AYM) invited Prof. Nir Navon (Yale University) to hold a colloquium for the physics department. For an audience of students and researchers, Prof. Navon presented recent advances in ultracold quantum matter and research from his own lab. His experimental work paves the way to make toy models used by theorists a reality. We sat down afterwards to discuss ultracold physics, box traps and setting up a lab from scratch.