# Alexey A. Sokol
Alexey A. Sokol describes himself as a theoretical physicist working on solid state physics, chemistry, and materials science, with a particular interest in localized states and defects and in the development of hybrid QM/MM embedded cluster techniques.[1](https://profiles.ucl.ac.uk/8781-alexey-sokols) • [2](https://orcid.org/0000-0003-0178-1147) His current research applies these methods and semi-classical atomistic models to nanoporous catalysts, metal oxides, and wide-gap semiconductors.[1](https://profiles.ucl.ac.uk/8781-alexey-sokols)
| Fact | Detail |
|---|---|
| Position | Principal Research Associate, Materials Chemistry Section, UCL Chemistry, since 2016[1](https://profiles.ucl.ac.uk/8781-alexey-sokols) |
| Field | Theoretical physicist; localized states and defects, and hybrid QM/MM embedded cluster techniques[2](https://orcid.org/0000-0003-0178-1147) |
| Known for | Hybrid QM/MM embedded cluster methods and the ChemShell software for defect and catalysis modeling[2](https://orcid.org/0000-0003-0178-1147) • [3](https://pubs.rsc.org/en/content/articlehtml/2023/cp/d3cp00648d) |
| Signature work | The 2004 hybrid QM/MM embedding method for ionic surfaces[4](https://doi.org/10.1002/qua.20032) |
| Training | PhD, Physical Chemistry, University of London, on defects in zeolites, supervised by Richard Catlow at the Royal Institution[1](https://profiles.ucl.ac.uk/8781-alexey-sokols) |
| Current focus | Predictive multiscale free energy simulations of hybrid transition metal catalysts (EPSRC award to October 2026)[5](https://gtr.ukri.org/person/9BA9F540-1CD9-45A8-92E5-C3EDCB0BABFF) |
## Training
His doctoral work was carried out at the Royal Institution of Great Britain under Richard Catlow, on the theory of defects in zeolites; his UCL profile records the PhD as completed in 1997.[1](https://profiles.ucl.ac.uk/8781-alexey-sokols) The dissertation, "Defect structures in zeolite crystals", is recorded in UCL Discovery with a publication date of 2000.[6](https://discovery.ucl.ac.uk/id/eprint/10106347)
The thesis identified the major defect species in zeolites: aluminum framework substitutionals, which are the Brønsted acid sites, together with hydroxyl nests, and vicinal disilanols, while peroxide-containing defects were found to be the major Lewis acid sites.[6](https://discovery.ucl.ac.uk/id/eprint/10106347)
## Career
Sokol was a Research Associate at the [Royal Institution](https://www.edgepedia.org/royal-institution)'s Davy Faraday Research Laboratory from 28 February 1997 to 31 July 2007, ORCID listing the role as Post Doctoral Research Assistant.[1](https://profiles.ucl.ac.uk/8781-alexey-sokols) • [2](https://orcid.org/0000-0003-0178-1147) He continued in Catlow's group as a postdoctoral researcher and moved with the group to UCL Chemistry, joining the Department on 1 August 2007.[1](https://profiles.ucl.ac.uk/8781-alexey-sokols) He was promoted to Senior Research Associate in 2010 and to Principal Research Associate in 2016.[1](https://profiles.ucl.ac.uk/8781-alexey-sokols)
## Representative work
According to a 2004 paper in the International Journal of Quantum Chemistry, a hybrid QM/MM technique was implemented in ChemShell for reactions at the surfaces of ionic solids, including reconstructed polar surfaces and interfaces. Its applications centered on electron trapping at the oxygen-terminated polar surfaces of ZnO and related surface F centers, species that proved to be active catalytic centers in methanol synthesis over ZnO.[4](https://doi.org/10.1002/qua.20032) A subsequent Royal Society survey of computational approaches to heterogeneous catalysis describes Sokol's group's use of DFT-based molecular cluster and embedded cluster QM/MM techniques, with the embedding matrix simulated by shell model potentials, across three case studies: alkene epoxidation over the microporous TS-1 catalyst, methanol synthesis on ZnO and Cu/ZnO, and C–H bond activation over Li-doped MgO.[7](https://doi.org/10.1098/rsta.2004.1529)
## ChemShell and the embedded cluster method
**How the method works.** An embedded cluster calculation treats a defect and its immediate surroundings with a quantum mechanical method while representing the rest of the ionic crystal classically. In Py-ChemShell, ionically bonded systems such as transition metal oxides require an ionic embedding procedure in which pseudopotentials are applied to atoms in a boundary region to localize the electron density within the inner QM region, because link atoms are inappropriate at such a boundary.[3](https://pubs.rsc.org/en/content/articlehtml/2023/cp/d3cp00648d) A doctoral thesis from this line of work set out the design goals explicitly: a method and code for embedded cluster calculations of point defects in the bulk and at surfaces of ionic crystals, able to handle charged defects and spectroscopic defect properties.[8](https://discovery.ucl.ac.uk/id/eprint/10106319/1/Development_and_application_of.pdf)
**The software.** ChemShell is a scriptable computational chemistry environment emphasizing multiscale [QM/MM simulation](https://www.edgepedia.org/qm-mm-simulation); it was redeveloped from the ground up as an open-source, Python-based platform (Py-ChemShell) for modeling chemical reactions on surfaces and within microporous solids on massively parallel computing systems.[9](https://doi.org/10.1021/acs.jctc.8b01036) A 2023 review in [Physical Chemistry Chemical Physics](https://www.edgepedia.org/physical-chemistry-chemical-physics) describes ChemShell as a leading software package for QM/MM calculations in catalysis modeling.[3](https://pubs.rsc.org/en/content/articlehtml/2023/cp/d3cp00648d)
## Embedded cluster versus supercell methods
The main alternative for defect calculations is the periodic supercell. Simulating an isolated defect with periodic plane-wave methods usually requires large supercells to avoid interactions between periodic images, and charged defects additionally need corrections for the long-range electrostatic interaction between those images; embedded cluster calculations instead exploit the locality of the defect, representing the environment by an embedding potential.[10](http://dollywood.itp.tuwien.ac.at/%7Eflorian/1.4922260.pdf) A 2015 implementation of density functional embedding theory in [VASP](https://www.edgepedia.org/vasp-software) demonstrated that embedded cluster models can reproduce the electronic structure of point defects in bulk semiconductors.[10](http://dollywood.itp.tuwien.ac.at/%7Eflorian/1.4922260.pdf)
Direct benchmarks exist on both sides. A 2017 study in Theoretical Chemistry Accounts compared cluster and supercell approaches directly, using defects in diamond as the test case.[11](https://doi.org/10.1007/s00214-017-2071-5)
The supercell side carries its own unresolved error problem. A 2009 [Physical Review Letters](https://www.edgepedia.org/physical-review-letters) paper states that, despite numerous attempts, a general scheme to correct finite-size errors in charged-defect supercell calculations was not yet available, and proposes an efficient method based on a rigorous analysis of electrostatics.[12](https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.102.016402) A methodological review goes further, arguing that any analytical error-correction scheme relying on electrostatic considerations is not appropriate for deriving reliable defect formation energies, especially for relaxed geometries, and proposing finite-size scaling instead, demonstrated on III–V semiconductors.[13](https://doi.org/10.1088/0965-0393/17/8/084003)
## Recent developments
The 2023 PCCP review consolidated the state of ChemShell for catalysis modeling and cataloged recent applications of the ionic QM/MM embedded cluster approach: characterizing native point defects in GaN, creating optimized models of rutile TiO2 surfaces, studying oxygen vacancies in TiO2 with DFT and high-level wavefunction methods, improving interatomic potentials for CeO2, investigating vacancies in MnO for CO2 transformation, and studying the defect properties of Cu in ZnO, an industrial methanol-synthesis catalyst.[3](https://pubs.rsc.org/en/content/articlehtml/2023/cp/d3cp00648d) UKRI records an EPSRC award of £834,868 to [University College London](https://www.edgepedia.org/university-college-london) running from 22 April 2024 to October 2026 for "Predictive multiscale free energy simulations of hybrid transition metal catalysts", listing Alexey Sokol.[5](https://gtr.ukri.org/person/9BA9F540-1CD9-45A8-92E5-C3EDCB0BABFF)
## References
1. Alexey Sokols | About | University College London. https://profiles.ucl.ac.uk/8781-alexey-sokols
2. Alexey A. Sokol (0000-0003-0178-1147) – ORCID. https://orcid.org/0000-0003-0178-1147
3. Multiscale QM/MM modelling of catalytic systems with ChemShell, Phys. Chem. Chem. Phys., 2023. https://pubs.rsc.org/en/content/articlehtml/2023/cp/d3cp00648d
4. Hybrid QM/MM embedding approach for the treatment of localized surface states in ionic materials, Int. J. Quantum Chem., 2004. https://doi.org/10.1002/qua.20032
5. Alexey Sokol – UKRI Gateway to Research. https://gtr.ukri.org/person/9BA9F540-1CD9-45A8-92E5-C3EDCB0BABFF
6. Defect structures in zeolite crystals – UCL Discovery. https://discovery.ucl.ac.uk/id/eprint/10106347
7. Computational approaches to the determination of active site structures and reaction mechanisms in heterogeneous catalysts, Phil. Trans. R. Soc. A. https://doi.org/10.1098/rsta.2004.1529
8. Development and Application of Embedded Cluster Methodologies for Defects in Ionic Materials (PhD thesis). https://discovery.ucl.ac.uk/id/eprint/10106319/1/Development_and_application_of.pdf
9. Open-Source, Python-Based Redevelopment of the ChemShell Multiscale QM/MM Environment, J. Chem. Theory Comput. https://doi.org/10.1021/acs.jctc.8b01036
10. Implementation of density functional embedding theory within the projector augmented-wave method and applications to semiconductor defect states, J. Chem. Phys., 2015. http://dollywood.itp.tuwien.ac.at/%7Eflorian/1.4922260.pdf
11. Comparison between cluster and supercell approaches: the case of defects in diamond, Theor. Chem. Acc., 2017. https://doi.org/10.1007/s00214-017-2071-5
12. Fully Ab Initio Finite-Size Corrections for Charged-Defect Supercell Calculations, Phys. Rev. Lett. 102, 016402, 2009. https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.102.016402
13. Density functional theory calculations of defect energies using supercells, Modelling Simul. Mater. Sci. Eng. https://doi.org/10.1088/0965-0393/17/8/084003
---
*Topic: Encyclopedia › Physical world and mathematics › Physical and mathematical scientists › Chemists › Researchers in physical, theoretical, and computational chemistry › Quantum chemistry and electronic structure theory › Density functional theory researchers*
*Initially written Sep 21, 2026 · Reviewed: Sep 25, 2026 · Edited: Sep 25, 2026; Oct 11, 2026; Oct 15, 2026; Oct 25, 2026 · Last review: Sep 25, 2026*
*Copyright 2026 EdgeChat AI, a subsidiary of Biostate AI.*
License: Edgepedia Community License 1.0, https://www.edgepedia.org/edgepedia/license