Source author record

Nicholas D. M. Hine

Nicholas D. M. Hine appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

19works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2022arXiv

Ab initio study on the atomic and electronic structures of twisted InSe bilayer

The electronic properties of the twisted InSe bilayer are studied by large-scale density functional theory. Spectral Function Unfolding reveals that the electronic structure of the twisted system can be described in terms of a combination of features of the bandstructures of the aligned InSe bilayer with different stacking configurations, enabling predictions of the band gap and the effective mass for holes. The effective mass for holes in the twisted InSe bilayer is shown to be influenced primarily by the interlayer distance. The intralayer and interlayer exciton binding energies are thus calculated based on a model recently developed by Ruiz-Tijerina et al. We apply similar analysis to the trilayer heterostructure InSe/hBN/InSe: its electronic structure is shown to be well-described by the superposition of band structures of two InSe monolayers with a small coupling through the hBN layer.

preprint2021arXiv

Ghost anti-crossings caused by interlayer umklapp hybridization of bands in 2D heterostructures

In two-dimensional heterostructures, crystalline atomic layers with differing lattice parameters can stack directly one on another. The resultant close proximity of atomic lattices with differing periodicity can lead to new phenomena. For umklapp processes, this opens the possibility for interlayer umklapp scattering, where interactions are mediated by the transfer of momenta to or from the lattice in the neighbouring layer. Using angle-resolved photoemission spectroscopy to study a graphene on InSe heterostructure, we present evidence that interlayer umklapp processes can cause hybridization between bands from neighbouring layers in regions of the Brillouin zone where bands from only one layer are expected, despite no evidence for moir/'e-induced replica bands. This phenomenon manifests itself as 'ghost' anti-crossings in the InSe electronic dispersion. Applied to a range of suitable 2DM pairs, this phenomenon of interlayer umklapp hybridization can be used to create strong mixing of their electronic states, giving a new tool for twist-controlled band structure engineering.

preprint2020arXiv

Atomic and electronic structure of two-dimensional Mo(1-x)WxS2 alloys

Alloying enables engineering of the electronic structure of semiconductors for optoelectronic applications. Due to their similar lattice parameters, the two-dimensional semiconducting transition metal dichalcogenides of the MoWSeS group (MX2 where M= Mo or W and X=S or Se) can be grown as high-quality materials with low defect concentrations. Here we investigate the atomic and electronic structure of Mo(1-x)WxS2 alloys using a combination of high-resolution experimental techniques and simulations. Analysis of the Mo and W atomic positions in these alloys, grown by chemical vapour transport, shows that they are randomly distributed, consistent with Monte Carlo simulations that use interaction energies determined from first-principles calculations. Electronic structure parameters are directly determined from angle resolved photoemission spectroscopy measurements. These show that the spin-orbit splitting at the valence band edge increases linearly with W content from MoS2 to WS2, in agreement with linear-scaling density functional theory (LS-DFT) predictions. The spin-orbit splitting at the conduction band edge is predicted to reduce to zero at intermediate compositions. Despite this, polarisation-resolved photoluminescence spectra on monolayer Mo0.5W0.5S2 show significant circular dichroism, indicating that spin-valley locking is retained. These results demonstrate that alloying is an important tool for controlling the electronic structure of MX2 for spintronic and valleytronic applications.

preprint2020arXiv

ONETEP + TOSCAM: uniting dynamical mean field theory and linear-scaling density functional theory

We introduce the unification of dynamical mean field theory (DMFT) and linear-scaling density functional theory (DFT), as recently implemented in ONETEP, a linear-scaling DFT package, and TOSCAM, a DMFT toolbox. This code can account for strongly correlated electronic behavior while simultaneously including the effects of the environment, making it ideally suited for studying complex and heterogeneous systems containing transition metals and lanthanides, such as metalloproteins. We systematically introduce the necessary formalism, which must account for the non-orthogonal basis set used by ONETEP. In order to demonstrate the capabilities of this code, we apply it to carbon monoxide-ligated iron porphyrin and explore the distinctly quantum-mechanical character of the iron $3d$ electrons during the process of photodissociation.

preprint2016arXiv

Band parameters and hybridization in 2D semiconductor heterostructures from photoemission spectroscopy

Combining monolayers of different two-dimensional (2D) semiconductors into heterostructures opens up a wealth of possibilities for novel electronic and optical functionalities. Exploiting them hinges on accurate measurements of the band parameters and orbital hybridization in separate and stacked monolayers, many of which are only available as small samples. The recently introduced technique of angle-resolved photoemission spectroscopy with submicron spatial resolution (μ-ARPES) offers the capability to measure small samples, but the energy resolution obtained for such exfoliated samples to date (~0.5 eV) has been inadequate. Here, we show that by suitable heterostructure sample design the full potential of μ-ARPES can be realized. We focus on MoSe2/WSe2 van der Waals heterostructures, which are 2D analogs of 3D semiconductor heterostructures. We find that in a MoSe2/WSe2 heterobilayer the bands in the K valleys are weakly hybridized, with the conduction and valence band edges originating in the MoSe2 and WSe2 respectively. There is stronger hybridization at the Γ point, but the valence band edge remains at the K points. This is consistent with the recent observation of interlayer excitons where the electron and hole are valley polarized but in opposite layers. We determine the valence band offset to be 300 meV, which combined with photoluminescence measurements implies that the binding energy of interlayer excitons is at least 200 meV, comparable with that of intralayer excitons.

preprint2016arXiv

Supercell convergence of charge-transfer energies in pentacene molecular crystals from constrained DFT

Singlet fission (SF) is a multi-exciton generation process that could be harnessed to improve the efficiency of photovoltaic devices. Experimentally, systems derived from the pentacene molecule have been shown to exhibit ultrafast SF with high yields. Charge-transfer (CT) configurations are likely to play an important role as intermediates in the SF process in these systems. In molecular crystals, electrostatic screening effects and band formation can be significant in lowering the energy of CT states, enhancing their potential to effectively participate in SF. In order to simulate these, it desirable to adopt a computational approach which is acceptably accurate, relatively inexpensive, which and scales well to larger systems, thus enabling the study of screening effects. We propose a novel, electrostatically-corrected constrained Density Functional Theory (cDFT) approach as a low-cost solution to the calculation of CT energies in molecular crystals such as pentacene. Here we consider an implementation in the context of the ONETEP linear-scaling DFT code, but our electrostatic correction method is in principle applicable in combination with any constrained DFT implementation, also outside the linear-scaling framework. Our newly developed method allows us to estimate CT energies in the infinite crystal limit, and with these to validate the accuracy of the cluster approximation.

preprint2015arXiv

Linear-scaling time-dependent density-functional theory (TDDFT) beyond the Tamm-Dancoff approximation: obtaining efficiency and accuracy with in situ optimised local orbitals

We present a solution of the full TDDFT eigenvalue equation in the linear response formalism exhibiting a linear-scaling computational complexity with system size, without relying on the simplifying Tamm-Dancoff approximation (TDA). The implementation relies on representing the occupied and unoccupied subspace with two different sets of in situ optimised localised functions, yielding a very compact and efficient representation of the transition density matrix of the excitation with the accuracy associated with a systematic basis set. The TDDFT eigenvalue equation is solved using a preconditioned conjugate-gradients algorithm that is very memory-efficient. The algorithm is validated on a test molecule and a good agreement with results obtained from standard quantum chemistry packages is found, with the preconditioner yielding a significant improvement in convergence rates. The method developed in this work is then used to reproduce experimental results of the absorption spectrum of bacteriochlorophyll (BChl) in an organic solvent, where it is demonstrated that the TDA fails to reproduce the main features of the low energy spectrum, while the full TDDFT equation yields results in good qualitative agreement with experimental data. Furthermore, the need for explicitly including parts of the solvent into the TDDFT calculations is highlighted, making the treatment of large system sizes necessary that are well within reach of the capabilities of the algorithm introduced here. Finally, the linear-scaling properties of the algorithm are demonstrated by computing the lowest excitation energy of BChl in solution. The largest systems considered in this work are of the same order of magnitude as a variety of pigment-protein complexes, opening up the possibility of studying their properties without having to resort to any semiclassical approximations to parts of the protein environment.

preprint2015arXiv

The elevated Curie temperature and half-metallicity in the ferromagnetic semiconductor La$_{x}$Eu$_{1-x}$O

Here we study the effect of La doping in EuO thin films using SQUID magnetometry, muon spin rotation ($μ$SR), polarized neutron reflectivity (PNR), and density functional theory (DFT). The $μ$SR data shows that the La$_{0.15}$Eu$_{0.85}$O is homogeneously magnetically ordered up to its elevated $T_{\rm C}$. It is concluded that bound magnetic polaron behavior does not explain the increase in $T_{\rm C}$ and an RKKY-like interaction is consistent with the $μ$SR data. The estimation of the magnetic moment by DFT simulations concurs with the results obtained by PNR, showing a reduction of the magnetic moment per La$_{x}$Eu$_{1-x}$O for increasing lanthanum doping. This reduction of the magnetic moment is explained by the reduction of the number of Eu-4$f$ electrons present in all the magnetic interactions in EuO films. Finally, we show that an upwards shift of the Fermi energy with La or Gd doping gives rise to half-metallicity for doping levels as high as 3.2 %.

preprint2015arXiv

Tracing potential energy surfaces of electronic excitations via their transition origins: application to Oxirane

We show that the transition origins of electronic excitations identified by quantified natural transition orbital (QNTO) analysis can be employed to connect potential energy surfaces (PESs) according to their character across a widerange of molecular geometries. This is achieved by locating the switching of transition origins of adiabatic potential surfaces as the geometry changes. The transition vectors for analysing transition origins are provided by linear response time-dependent density functional theory (TDDFT) calculations under the Tamm-Dancoff approximation. We study the photochemical CO ring opening of oxirane as an example and show that the results corroborate the traditional Gomer-Noyes mechanism derived experimentally. The knowledge of specific states for the reaction also agrees well with that given by previous theoretical work using TDDFT surface-hopping dynamics that was validated by high-quality quantum Monte Carlo calculations. We also show that QNTO can be useful for considerably larger and more complex systems: by projecting the excitations to those of a reference oxirane molecule, the approach is able to identify and analyse specific excitations of a trans-2,3-diphenyloxirane molecule.

preprint2013arXiv

Anisotropic charge screening and supercell size convergence of defect formation energies

One of the main sources of error associated with the calculation of defect formation energies using plane-wave Density Functional Theory (DFT) is finite size error resulting from the use of relatively small simulation cells and periodic boundary conditions. Most widely-used methods for correcting this error, such as that of Makov and Payne, assume that the dielectric response of the material is isotropic and can be described using a scalar dielectric constant $ε$. However, this is strictly only valid for cubic crystals, and cannot work in highly-anisotropic cases. Here we introduce a variation of the technique of extrapolation based on the Madelung potential, that allows the calculation of well converged dilute limit defect formation energies in non-cubic systems with highly anisotropic dielectric properties. As an example of the implementation of this technique we study a selection of defects in the ceramic oxide Li$_2$TiO$_3$ which is currently being considered as a lithium battery material and a breeder material for fusion reactors.

preprint2013arXiv

Importance of many body effects in the kernel of hemoglobin for ligand binding

We propose a mechanism for binding of diatomic ligands to heme based on a dynamical orbital selection process. This scenario may be described as bonding determined by local valence fluctuations. We support this model using linear-scaling first-principles calculations, in combination with dynamical mean-field theory, applied to heme, the kernel of the hemoglobin metalloprotein central to human respiration. We find that variations in Hund's exchange coupling induce a reduction of the iron 3d density, with a concomitant increase of valence fluctuations. We discuss the comparison between our computed optical absorption spectra and experimental data, our picture accounting for the observation of optical transitions in the infrared regime, and how the Hund's coupling reduces, by a factor of five, the strong imbalance in the binding energies of heme with CO and O_2 ligands.

preprint2013arXiv

Linear-scaling time-dependent density-functional theory in the linear response formalism

We present an implementation of time-dependent density-functional theory (TDDFT) in the linear response formalism enabling the calculation of low energy optical absorption spectra for large molecules and nanostructures. The method avoids any explicit reference to canonical representations of either occupied or virtual Kohn-Sham states and thus achieves linear-scaling computational effort with system size. In contrast to conventional localised orbital formulations, where a single set of localised functions is used to span the occupied and unoccupied state manifold, we make use of two sets of in situ optimised localised orbitals, one for the occupied and one for the unoccupied space. This double representation approach avoids known problems of spanning the space of unoccupied Kohn-Sham states with a minimal set of localised orbitals optimised for the occupied space, while the in situ optimisation procedure allows for efficient calculations with a minimal number of functions. The method is applied to a number of medium sized organic molecules and a good agreement with traditional TDDFT methods is observed. Furthermore, linear scaling of computational cost with system size is demonstrated on a system of carbon nanotubes.

preprint2013arXiv

Simulations of Nanocrystals Under Pressure: Combining Electronic Enthalpy and Linear-Scaling Density-Functional Theory

We present an implementation in a linear-scaling density-functional theory code of an electronic enthalpy method, which has been found to be natural and efficient for the ab initio calculation of finite systems under hydrostatic pressure. Based on a definition of the system volume as that enclosed within an electronic density isosurface [Phys. Rev. Lett., 94, 145501 (2005)], it supports both geometry optimizations and molecular dynamics simulations. We introduce an approach for calibrating the parameters defining the volume in the context of geometry optimizations and discuss their significance. Results in good agreement with simulations using explicit solvents are obtained, validating our approach. Size-dependent pressure-induced structural transformations and variations in the energy gap of hydrogenated silicon nanocrystals are investigated, including one comparable in size to recent experiments. A detailed analysis of the polyamorphic transformations reveals three types of amorphous structures and their persistence on depressurization is assessed.

preprint2012arXiv

Fermi level pinning can determine polarity in semiconductor nanorods

First-principles calculations of polar semiconductor nanorods reveal that their dipole moments are strongly influenced by Fermi level pinning. The Fermi level for an isolated nanorod is found to coincide with a significant density of electronic surface states at the end surfaces, which are either mid-gap states or band-edge states. These states pin the Fermi level, and therefore fix the potential difference across the rod. We provide evidence that this effect can have a determining influence on the polarity of nanorods, and has consequences for the way a rod responds to changes in its surface chemistry, the scaling of its dipole moment with its size, and the dependence of polarity on its composition.

preprint2012arXiv

Linear-scaling DFT+U with full local orbital optimization

We present an approach to the DFT+U method (Density Functional Theory + Hubbard model) within which the computational effort for calculation of ground state energies and forces scales linearly with system size. We employ a formulation of the Hubbard model using nonorthogonal projector functions to define the localized subspaces, and apply it to a local-orbital DFT method including in situ orbital optimization. The resulting approach thus combines linear-scaling and systematic variational convergence. We demonstrate the scaling of the method by applying it to nickel oxide nano-clusters with sizes exceeding 7,000 atoms.

preprint2012arXiv

Vanadium dioxide : A Peierls-Mott insulator stable against disorder

Vanadium dioxide undergoes a first order metal-insulator transition at 340 K. In this work, we develop and carry out state of the art linear scaling DFT calculations refined with non-local dynamical mean-field theory. We identify a complex mechanism, a Peierls-assisted orbital selection Mott instability, which is responsible for the insulating M$_1$ phase, and furthermore survives a moderate degree of disorder.

preprint2011arXiv

Calculating optical absorption spectra for large systems using linear-scaling density-functional theory

A new method for calculating optical absorption spectra within linear-scaling density-functional theory (LS-DFT) is presented, incorporating a scheme for optimizing a set of localized orbitals to accurately represent unoccupied Kohn-Sham states. Three different schemes are compared and the most promising of these, based on the use of a projection operator, has been implemented in a fully-functional LS-DFT code. The method has been applied to the calculation of optical absorption spectra for the metal-free phthalocyanine molecule and the conjugated polymer poly(para-phenylene). Excellent agreement with results from a traditional DFT code is obtained.

preprint2011arXiv

Factors influencing the distribution of charge in polar nanocrystals

We perform first-principles calculations of wurtzite GaAs nanorods to explore the factors determining charge distributions in polar nanostructures. We show that both the direction and magnitude of the dipole moment $\mathbf{d}$ of a nanorod, and its electic field, depend sensitively on how its surfaces are terminated and do not depend strongly on the spontaneous polarization of the underlying lattice. We identify two physical mechanisms by which $\mathbf{d}$ is controlled by the surface termination, and we show that the excess charge on the nanorod ends is not strongly localized. We discuss the implications of these results for tuning nanocrystal properties, and for their growth and assembly.

preprint2010arXiv

Projector self-consistent DFT+U using non-orthogonal generalized Wannier functions

We present a formulation of the density-functional theory + Hubbard model (DFT+U) method that is self-consistent over the choice of Hubbard projectors used to define the correlated subspaces. In order to overcome the arbitrariness in this choice, we propose the use of non-orthogonal generalized Wannier functions (NGWFs) as projectors for the DFT+U correction. We iteratively refine these NGWF projectors and, hence, the DFT+U functional, such that the correlated subspaces are fully self-consistent with the DFT+U ground-state. We discuss the convergence characteristics of this algorithm and compare ground-state properties thus computed with those calculated using hydrogenic projectors. Our approach is implemented within, but not restricted to, a linear-scaling DFT framework, opening the path to DFT+U calculations on systems of unprecedented size.