Scientific Publications

Peer-reviewed work from my time in computational chemistry and nanoplasmonics, plus patents. Full list on Google Scholar.

2026

plasmonX: an Open-Source Code for Nanoplasmonics

This paper introduces plasmonX, an open-source software package that implements the family of atomistic electromagnetic models (ωFQ, ωFQFμ and their multiscale extensions) developed for simulating plasmonic nanostructures. The code allows users to compute the optical response of metal and graphene-based nanosystems at the atomistic level, in both the frequency and time domains, and to couple them with quantum-mechanical descriptions of nearby molecules. It is designed to make state-of-the-art nanoplasmonics modeling accessible and reproducible for the wider community.


2025

The Electric Field Morphology of Plasmonic Picocavities

Picocavities — plasmonic gaps containing atomic-scale defects — confine light below the nanometer scale, but how these defects shape the local electric field had been hard to study with traditional continuum models. Using the atomistic ωFQFμ approach, validated against quantum-mechanical calculations, the authors simulate gold picocavities made of thousands of atoms and map their field morphology in detail. The results reveal strong field inhomogeneities induced by individual atomic defects, with implications for molecular sensing, spectroscopy, and plasmon-driven chemistry, where field gradients play a key role.


2024

Atomistic Multiscale Modeling of Colloidal Plasmonic Nanoparticles

Real-world plasmonic nanoparticles operate in solution, yet most atomistic models treat them in vacuum. This work presents ωFQFμ/FQ, a multiscale classical model that couples an atomistic description of the nanoparticle’s plasmonic response with a polarizable force-field description of the surrounding solvent, making it possible to simulate real-size colloidal nanoparticles of generic chemical nature. Validated against reference quantum-mechanical calculations, the model accurately captures solvent-induced shifts of the plasmon resonance at a fraction of the computational cost.


Real-Time Formulation of Atomistic Electromagnetic Models for Plasmonics

Atomistic electromagnetic models of plasmonic nanoparticles are usually formulated in the frequency domain, which limits the study of time-resolved phenomena and ultrafast laser excitation. This paper reformulates the ωFQ and ωFQFμ models in the real-time domain, propagating the induced charges and dipoles of each atom as the system responds to an external pulse. The approach reproduces frequency-domain results while opening the way to simulating time-dependent plasmon dynamics, nonlinear regimes, and arbitrary excitation schemes in large nanostructures.


2023

QM/Classical Modeling of Surface Enhanced Raman Scattering Based on Atomistic Electromagnetic Models

Surface-enhanced Raman scattering (SERS) lets researchers detect molecules near plasmonic nanostructures, but simulating the effect requires treating both the molecule and the metal at high accuracy. This work combines a quantum-mechanical description of the adsorbed molecule with the atomistic classical models ωFQ and ωFQFμ for the plasmonic substrate, yielding a multiscale method that works for both noble-metal and graphene-based nanostructures. The approach captures the electromagnetic enhancement of Raman signals and provides an interpretable, computationally affordable tool for studying molecule–nanostructure interactions.


Fully Atomistic Modeling of Plasmonic Bimetallic Nanoparticles: Nanoalloys and Core-Shell Systems

Mixing two metals in a single nanoparticle — as a random alloy or a core-shell structure — offers a powerful way to tune plasmonic properties, but modeling such systems atomistically is challenging. This paper extends the ωFQFμ model to bimetallic gold–silver nanoparticles, generalizing its formulation so that each atom carries parameters of its own chemical species. Tested against reference quantum-mechanical calculations, the approach reproduces how composition and atomic arrangement affect the plasmonic response, enabling the in silico exploration of multimetallic nanostructure design.


2022

Do We Really Need Quantum Mechanics to Describe Plasmonic Properties of Metal Nanostructures?

Plasmonics is ultimately a quantum phenomenon, so it is natural to ask whether expensive quantum-mechanical calculations are indispensable to model metal nanostructures. This paper systematically compares the classical atomistic models ωFQ and ωFQFμ against fully quantum time-dependent DFT calculations for metal nanoparticles of decreasing size. The classical models reproduce the quantum reference remarkably well down to very small structures, provided quantum tunneling effects at sub-nanometer junctions are accounted for — showing that accurate, affordable classical atomistic modeling is a viable alternative in most plasmonic regimes.


Assessing the Quality of QM/MM Approaches to Describe Vacuo-to-Water Solvatochromic Shifts

When a molecule moves from vacuum into water, its absorption spectrum shifts — an effect that computational models should reproduce to be predictive in realistic environments. This study benchmarks a range of QM/MM embedding schemes, from simple nonpolarizable force fields to fully polarizable approaches, against experimental vacuo-to-water solvatochromic shifts for a set of chromophores. No single model emerges as uniformly best: the accuracy depends strongly on the nature of the electronic transition, providing practical guidance on which embedding strategy to choose for a given problem.


In Silico Design of Graphene Plasmonic Hot-Spots

Plasmonic hot-spots — tiny regions of enormously enhanced electric field — are the key to ultrasensitive molecular detection, and graphene could offer a cheaper, more versatile platform than noble metals. Using atomistic electromagnetic simulations, this work explores engineered graphene nanostructures, inspired by patterns known to enable single-molecule detection on metal substrates, and focuses on realistic features such as edge defects and grain boundaries. The calculations show that highly localized hot-spots with enhancement factors comparable to noble metals can form near these defects, outlining design rules for graphene-based plasmonic sensors.


2018

Modulation of the Nonlinear Optical Properties of Dibenzo[hi,st]ovalene by Peripheral Substituents

Dibenzo[hi,st]ovalene (DBOV) is a nanographene molecule with strong red emission, high photostability, and optical gain — properties that make it attractive for lasers and strong light–matter coupling in microcavities. This study investigates how different peripheral chemical substituents affect the stimulated emission and nonlinear optical behavior of three DBOV derivatives using ultrafast spectroscopy. The results show that the substituent pattern modulates intermolecular interactions and thus the nonlinear response, offering a chemical handle to tune these nanographenes for photonic applications.


Patents & Theses

Atomistic Models for Plasmonics: from Metal Nanostructures to Molecular Plasmonics

This thesis develops the ωMM family of classical atomistic electromagnetic models for simulating the optical properties of plasmonic materials with quantum-level accuracy at classical-level cost. Starting from the ωFQ model for Drude-like metals, it extends the framework to noble metals by including interband transitions, generalizes it to bimetallic and alloyed nanostructures, and reformulates it in the real-time domain for time-dependent phenomena. Finally, it couples the models to polarizable environments and quantum-mechanical molecular descriptions, enabling applications such as refractive-index sensing and surface-enhanced Raman scattering.


In Silico Process for Identifying the Design of a Nanostructured Plasmonic Sensor

This patent protects a fully computational workflow for designing plasmonic sensors with single-molecule sensitivity. Starting from a target analyte and a library of nanostructured plasmonic substrates, the method maps the analyte–substrate interaction potential at the atomistic level, computes the resulting spectroscopic signal and its enhancement factor, and compares it against a predefined sensitivity threshold. The result is a screening procedure that identifies promising sensor designs before any fabrication, reducing experimental trial and error.