@article{21381, abstract = {The lack of long-range electrostatics is a key limitation of modern machine learning interatomic potentials (MLIPs), hindering reliable applications to interfaces, charge-transfer reactions, polar and ionic materials, and biomolecules. In this Perspective, we distill two design principles behind the Latent Ewald Summation framework, which can capture long-range interactions, charges, and electrical response just by learning from standard energy and force training data: (i) use a Coulomb functional form with environment-dependent charges to capture electrostatic interactions, and (ii) avoid explicit training on ambiguous density functional theory partial charges. When both principles are satisfied, substantial flexibility remains: essentially any short-range MLIP can be augmented; charge equilibration schemes can be added when desired; dipoles and Born effective charges can be inferred or fine-tuned; and charge/spin-state embeddings or tensorial targets can be further incorporated. We also discuss current limitations and open challenges. Together, these minimal, physics-guided design rules suggest that incorporating long-range electrostatics into MLIPs is simpler and perhaps more broadly applicable than is commonly assumed.}, author = {Kim, Dongjin and Cheng, Bingqing}, issn = {1089-7690}, journal = {The Journal of Chemical Physics}, number = {6}, publisher = {AIP Publishing}, title = {{Long-range electrostatics for machine learning interatomic potentials is easier than we thought}}, doi = {10.1063/5.0316886}, volume = {164}, year = {2026}, } @article{21408, abstract = {Rational design strategies for self-assembly require a detailed understanding of both the equilibrium state and the assembly kinetics. While the former is starting to be well understood, the latter remains a major theoretical challenge, especially in programmable systems and the so-called semi-addressable regime, where binding is often nondeterministic and the formation of off-target structures negatively influences the assembly. Here, we show that it is possible to simultaneously sculpt the assembly outcome and the assembly kinetics through the underexplored design space of binding energies and particle concentrations. By formulating the assembly process as a complex reaction network, we calculate and optimize the tradeoff between assembly speed and quality and show that parameter optimization can speed up assembly by many orders of magnitude without lowering the yield of the target structure. Although the exact speedup varies from design to design, we find the largest speedups for nondeterministic systems where unoptimized assembly is the slowest, sometimes even making them assemble faster than optimized, fully addressable designs. Therefore, these results not only solve a key challenge in semi-addressable self-assembly but further emphasize the utility of semi-addressability, where designs have the potential to be faster as well as cheaper (fewer particle species) and better (higher yield). More broadly, our results highlight the importance of parameter optimization in programmable self-assembly and provide practical tools for simultaneous optimization of kinetics and yield in a wide range of systems.}, author = {Hübl, Maximilian and Goodrich, Carl Peter}, issn = {1089-7690}, journal = {Journal of Chemical Physics}, number = {8}, publisher = {AIP Publishing}, title = {{Simultaneous optimization of assembly time and yield in programmable self-assembly}}, doi = {10.1063/5.0304731}, volume = {164}, year = {2026}, } @article{21840, abstract = {The transport properties of nanofluidic channels are usually studied under constant (DC) voltage or pressure driving. However, the frequency response under sinusoidal (AC) drivings offers rich insights into the time-dependent transport mechanisms. Inspired by recent electrochemical approaches, we investigate the couplings between ionic and electronic transport under AC driving. We show that conduction electrons of the channel walls participate in ionic current via capacitive electrochemical coupling, defining a critical frequency and length scale where electron-dominated conductivity emerges. We further analyze how electron–ion coupling modifies electro-osmotic flows and demonstrate that fluctuation-induced momentum transfer between the electrolyte and wall electrons produces distinct AC transport signatures, depending on the charge carrier polarity. Altogether, we establish a frequency-dependent transport matrix that couples ionic, electronic, and hydrodynamic flows. These findings establish AC nanofluidic transport as a powerful probe of interfacial phenomena under confinement and suggest new directions for engineering nanofluidic functionalities through electron–electrolyte coupling.}, author = {Coquinot, Baptiste and Lizée, Mathieu and Bocquet, Lydéric and Kavokine, Nikita}, issn = {1089-7690}, journal = {The Journal of Chemical Physics}, number = {13}, publisher = {AIP Publishing}, title = {{Electron–electrolyte coupling in AC transport through nanofluidic channels}}, doi = {10.1063/5.0313352}, volume = {164}, year = {2026}, } @article{19880, abstract = {We investigate quantum transport in a two-dimensional electron system coupled to a chiral molecular potential, demonstrating how molecular chirality and orientation affect charge and spin transport properties. We propose a minimal model for realizing true chiral symmetry breaking on a magnetized surface, with a crucial role played by the tilt angle of the molecular dipole with respect to the surface. For non-zero tilting, we show that the Hall response exhibits clear signatures of chirality-induced effects, in both charge- and spin-resolved observables. Concerning the former, tilted enantiomers produce asymmetric Hall conductances and, even more remarkably, the persistence of this feature in the absence of spin–orbit coupling (SOC) signals how the enantiospecific charge response results from electron scattering off the molecular potential. Concerning spin-resolved observables where SOC plays a relevant role, we reveal that chiral symmetry breaking is crucial in enabling spin-flipping processes.}, author = {Al Hyder, Ragheed and Lemeshko, Mikhail and Cappellaro, Alberto}, issn = {1089-7690}, journal = {The Journal of Chemical Physics}, number = {23}, publisher = {AIP Publishing}, title = {{Quantum transport in the presence of a chiral molecular potential}}, doi = {10.1063/5.0271155}, volume = {162}, year = {2025}, } @article{19276, abstract = {Impurity motion in a many-body environment has been a central issue in the field of low-temperature physics for decades. In bosonic quantum fluids, the onset of a drag force experienced by point-like objects is due to collective environment excitations, driven by the exchange of linear momentum between the impurity and the many-body bath. In this work we consider a rotating impurity, with the aim of exploring how angular momentum is exchanged with the surrounding bosonic environment. In order to elucidate these issues, we employ a quasiparticle approach based on the angulon theory, which allows us to effectively deal with the non-trivial algebra of quantized angular momentum in the presence of a many-body environment. We uncover how impurity dressing by environmental excitations can establish an exchange channel, whose effectiveness crucially depends on the initial state of the impurity. Remarkably, we find that there is a critical value of initial angular momentum, above which this channel effectively freezes.}, author = {Cappellaro, Alberto and Bighin, Giacomo and Cherepanov, Igor and Lemeshko, Mikhail}, issn = {1089-7690}, journal = {Journal of Chemical Physics}, number = {7}, publisher = {AIP Publishing}, title = {{Environment-limited transfer of angular momentum in Bose liquids}}, doi = {10.1063/5.0253451}, volume = {162}, year = {2025}, } @article{19279, abstract = {Recent experimental advances in nanofluidics have allowed to explore ion transport across molecular-scale pores, in particular, for iontronic applications. Two-dimensional nanochannels—in which a single molecular layer of electrolyte is confined between solid walls—constitute a unique platform to investigate fluid and ion transport in extreme confinement, highlighting unconventional transport properties. In this work, we study ionic association in 2D nanochannels, and its consequences on non-linear ionic transport, using both molecular dynamics simulations and analytical theory. We show that under sufficient confinement, ions assemble into pairs or larger clusters in a process analogous to a Kosterlitz–Thouless transition, here modified by the dielectric confinement. We further show that the breaking of pairs results in an electric-field dependent conduction, a mechanism usually known as the second Wien effect. However the 2D nature of the system results in non-universal, temperature-dependent, scaling of the conductivity with electric field, leading to ionic coulomb blockade in some regimes. A 2D generalization of the Onsager theory fully accounts for the non-linear transport. These results suggest ways to exploit electrostatic interactions between ions to build new nanofluidic devices.}, author = {Toquer, Damien and Bocquet, Lydéric and Robin, Paul}, issn = {1089-7690}, journal = {Journal of Chemical Physics}, number = {6}, publisher = {AIP Publishing}, title = {{Ionic association and Wien effect in 2D confined electrolytes}}, doi = {10.1063/5.0241949}, volume = {162}, year = {2025}, } @article{19595, abstract = {We investigate the locality of magnetic response in polycyclic aromatic molecules using a novel deep-learning approach. Our method employs graph neural networks (GNNs) with a graph-of-rings representation to predict nucleus independent chemical shifts (NICS) in the space around the molecule. We train a series of models, each time reducing the size of the largest molecules used in training. The accuracy of prediction remains high (MAE < 0.5 ppm), even when training the model only on molecules with up to four rings, thus providing strong evidence for the locality of magnetic response. To overcome the known problem of generalization of GNNs, we implement a k-hop expansion strategy and succeed in achieving accurate predictions for molecules with up to 15 rings (almost 4 times the size of the largest training example). Our findings have implications for understanding the magnetic response in complex molecules and demonstrate a promising approach to overcoming GNN scalability limitations. Furthermore, the trained models enable rapid characterization, without the need for more expensive DFT calculations.}, author = {Davidson, Yair and Philipp, Aviad and Chakraborty, Sabyasachi and Bronstein, Alexander and Gershoni-Poranne, Renana}, issn = {1089-7690}, journal = {Journal of Chemical Physics}, number = {14}, publisher = {AIP Publishing}, title = {{How local is “local”? Deep learning reveals locality of the induced magnetic field of polycyclic aromatic hydrocarbons}}, doi = {10.1063/5.0257558}, volume = {162}, year = {2025}, } @article{15024, abstract = {Electrostatic correlations between ions dissolved in water are known to impact their transport properties in numerous ways, from conductivity to ion selectivity. The effects of these correlations on the solvent itself remain, however, much less clear. In particular, the addition of salt has been consistently reported to affect the solution’s viscosity, but most modeling attempts fail to reproduce experimental data even at moderate salt concentrations. Here, we use an approach based on stochastic density functional theory, which accurately captures charge fluctuations and correlations. We derive a simple analytical expression for the viscosity correction in concentrated electrolytes, by directly linking it to the liquid’s structure factor. Our prediction compares quantitatively to experimental data at all temperatures and all salt concentrations up to the saturation limit. This universal link between the microscopic structure and viscosity allows us to shed light on the nanoscale dynamics of water and ions under highly concentrated and correlated conditions.}, author = {Robin, Paul}, issn = {1089-7690}, journal = {Journal of Chemical Physics}, number = {6}, publisher = {AIP Publishing}, title = {{Correlation-induced viscous dissipation in concentrated electrolytes}}, doi = {10.1063/5.0188215}, volume = {160}, year = {2024}, } @article{17278, abstract = {An azeotrope is a constant boiling point mixture, and its behavior is important for fluid separation processes. Predicting azeotropes from atomistic simulations is difficult due to the complexities and convergence problems of Monte Carlo and free-energy perturbation techniques. Here, we present a methodology for predicting the azeotropes of binary mixtures, which computes the compositional dependence of chemical potentials from molecular dynamics simulations using the S0 method and employs experimental boiling point and vaporization enthalpy data. Using this methodology, we reproduce the azeotropes, or lack thereof, in five case studies, including ethanol/water, ethanol/isooctane, methanol/water, hydrazine/water, and acetone/chloroform mixtures. We find that it is crucial to use the experimental boiling point and vaporization enthalpy for reliable azeotrope predictions, as empirical force fields are not accurate enough for these quantities. Finally, we use regular solution models to rationalize the azeotropes and reveal that they tend to form when the mixture components have similar boiling points and strong interactions.}, author = {Wang, Xiaoyu and Cheng, Bingqing}, issn = {1089-7690}, journal = {Journal of Chemical Physics}, number = {3}, publisher = {AIP Publishing}, title = {{Integrating molecular dynamics simulations and experimental data for azeotrope predictions in binary mixtures}}, doi = {10.1063/5.0217232}, volume = {161}, year = {2024}, } @article{14321, abstract = {We demonstrate the possibility of a coupling between the magnetization direction of a ferromagnet and the tilting angle of adsorbed achiral molecules. To illustrate the mechanism of the coupling, we analyze a minimal Stoner model that includes Rashba spin–orbit coupling due to the electric field on the surface of the ferromagnet. The proposed mechanism allows us to study magnetic anisotropy of the system with an extended Stoner–Wohlfarth model and argue that adsorbed achiral molecules can change magnetocrystalline anisotropy of the substrate. Our research aims to motivate further experimental studies of the current-free chirality induced spin selectivity effect involving both enantiomers.}, author = {Al Hyder, Ragheed and Cappellaro, Alberto and Lemeshko, Mikhail and Volosniev, Artem}, issn = {1089-7690}, journal = {The Journal of Chemical Physics}, keywords = {Physical and Theoretical Chemistry, General Physics and Astronomy}, number = {10}, publisher = {AIP Publishing}, title = {{Achiral dipoles on a ferromagnet can affect its magnetization direction}}, doi = {10.1063/5.0165806}, volume = {159}, year = {2023}, } @article{14603, abstract = {Computing the solubility of crystals in a solvent using atomistic simulations is notoriously challenging due to the complexities and convergence issues associated with free-energy methods, as well as the slow equilibration in direct-coexistence simulations. This paper introduces a molecular-dynamics workflow that simplifies and robustly computes the solubility of molecular or ionic crystals. This method is considerably more straightforward than the state-of-the-art, as we have streamlined and optimised each step of the process. Specifically, we calculate the chemical potential of the crystal using the gas-phase molecule as a reference state, and employ the S0 method to determine the concentration dependence of the chemical potential of the solute. We use this workflow to predict the solubilities of sodium chloride in water, urea polymorphs in water, and paracetamol polymorphs in both water and ethanol. Our findings indicate that the predicted solubility is sensitive to the chosen potential energy surface. Furthermore, we note that the harmonic approximation often fails for both molecular crystals and gas molecules at or above room temperature, and that the assumption of an ideal solution becomes less valid for highly soluble substances.}, author = {Reinhardt, Aleks and Chew, Pin Yu and Cheng, Bingqing}, issn = {1089-7690}, journal = {Journal of Chemical Physics}, number = {18}, publisher = {AIP Publishing}, title = {{A streamlined molecular-dynamics workflow for computing solubilities of molecular and ionic crystals}}, doi = {10.1063/5.0173341}, volume = {159}, year = {2023}, } @article{12705, abstract = {The elasticity of disordered and polydisperse polymer networks is a fundamental problem of soft matter physics that is still open. Here, we self-assemble polymer networks via simulations of a mixture of bivalent and tri- or tetravalent patchy particles, which result in an exponential strand length distribution analogous to that of experimental randomly cross-linked systems. After assembly, the network connectivity and topology are frozen and the resulting system is characterized. We find that the fractal structure of the network depends on the number density at which the assembly has been carried out, but that systems with the same mean valence and same assembly density have the same structural properties. Moreover, we compute the long-time limit of the mean-squared displacement, also known as the (squared) localization length, of the cross-links and of the middle monomers of the strands, showing that the dynamics of long strands is well described by the tube model. Finally, we find a relation connecting these two localization lengths at high density and connect the cross-link localization length to the shear modulus of the system.}, author = {Sorichetti, Valerio and Ninarello, Andrea and Ruiz-Franco, José and Hugouvieux, Virginie and Zaccarelli, Emanuela and Micheletti, Cristian and Kob, Walter and Rovigatti, Lorenzo}, issn = {1089-7690}, journal = {Journal of Chemical Physics}, number = {7}, publisher = {American Institute of Physics}, title = {{Structure and elasticity of model disordered, polydisperse, and defect-free polymer networks}}, doi = {10.1063/5.0134271}, volume = {158}, year = {2023}, } @article{12249, abstract = {The chemical potential of a component in a solution is defined as the free energy change as the amount of that component changes. Computing this fundamental thermodynamic property from atomistic simulations is notoriously difficult because of the convergence issues involved in free energy methods and finite size effects. This Communication presents the so-called S0 method, which can be used to obtain chemical potentials from static structure factors computed from equilibrium molecular dynamics simulations under the isothermal–isobaric ensemble. This new method is demonstrated on the systems of binary Lennard-Jones particles, urea–water mixtures, a NaCl aqueous solution, and a high-pressure carbon–hydrogen mixture. }, author = {Cheng, Bingqing}, issn = {1089-7690}, journal = {The Journal of Chemical Physics}, keywords = {Physical and Theoretical Chemistry, General Physics and Astronomy}, number = {12}, publisher = {AIP Publishing}, title = {{Computing chemical potentials of solutions from structure factors}}, doi = {10.1063/5.0107059}, volume = {157}, year = {2022}, } @article{11400, abstract = {By varying the concentration of molecules in the cytoplasm or on the membrane, cells can induce the formation of condensates and liquid droplets, similar to phase separation. Their thermodynamics, much studied, depends on the mutual interactions between microscopic constituents. Here, we focus on the kinetics and size control of 2D clusters, forming on membranes. Using molecular dynamics of patchy colloids, we model a system of two species of proteins, giving origin to specific heterotypic bonds. We find that concentrations, together with valence and bond strength, control both the size and the growth time rate of the clusters. In particular, if one species is in large excess, it gradually saturates the binding sites of the other species; the system then becomes kinetically arrested and cluster coarsening slows down or stops, thus yielding effective size selection. This phenomenology is observed both in solid and fluid clusters, which feature additional generic homotypic interactions and are reminiscent of the ones observed on biological membranes.}, author = {Palaia, Ivan and Šarić, Anđela}, issn = {1089-7690}, journal = {The Journal of Chemical Physics}, keywords = {Physical and Theoretical Chemistry, General Physics and Astronomy}, number = {19}, publisher = {AIP Publishing}, title = {{Controlling cluster size in 2D phase-separating binary mixtures with specific interactions}}, doi = {10.1063/5.0087769}, volume = {156}, year = {2022}, } @article{17907, abstract = {Carbyne is a linear allotrope of carbon that is composed of a chain of sp-hybridized carbon atoms. Through appropriate engineering of the chain termination, carbyne can harbor helical states where the π-electron delocalization twists along the axis of the chain. Herein, we present a comprehensive analysis of these helical states at the tight-binding level. We demonstrate that, in general, the molecular orbital coefficients of the helical states trace out an ellipse, in analogy to elliptically polarized light. Helical states can be realized in a model, inspired by the structure of cumulene, which considers a chain terminated by sp2-hybridized atoms oriented at a nontrivial dihedral angle. We provide a complete analytic solution for this model. Additionally, we present a variation of the model that yields perfect helical states that trace out a circle as opposed to an ellipse. Our results provide a deeper understanding of helical states and lay a foundation for more advanced levels of theory.}, author = {Gunasekaran, Suman and Venkataraman, Latha}, issn = {1089-7690}, journal = {The Journal of Chemical Physics}, number = {12}, publisher = {AIP Publishing}, title = {{Tight-binding analysis of helical states in carbyne}}, doi = {10.1063/5.0021146}, volume = {153}, year = {2020}, } @article{9658, abstract = {Macroscopic models of nucleation provide powerful tools for understanding activated phase transition processes. These models do not provide atomistic insights and can thus sometimes lack material-specific descriptions. Here, we provide a comprehensive framework for constructing a continuum picture from an atomistic simulation of homogeneous nucleation. We use this framework to determine the equilibrium shape of the solid nucleus that forms inside bulk liquid for a Lennard-Jones potential. From this shape, we then extract the anisotropy of the solid-liquid interfacial free energy, by performing a reverse Wulff construction in the space of spherical harmonic expansions. We find that the shape of the nucleus is nearly spherical and that its anisotropy can be perfectly described using classical models.}, author = {Cheng, Bingqing and Ceriotti, Michele and Tribello, Gareth A.}, issn = {1089-7690}, journal = {The Journal of Chemical Physics}, number = {4}, publisher = {AIP Publishing}, title = {{Classical nucleation theory predicts the shape of the nucleus in homogeneous solidification}}, doi = {10.1063/1.5134461}, volume = {152}, year = {2020}, } @article{8587, abstract = {Inspired by the possibility to experimentally manipulate and enhance chemical reactivity in helium nanodroplets, we investigate the effective interaction and the resulting correlations between two diatomic molecules immersed in a bath of bosons. By analogy with the bipolaron, we introduce the biangulon quasiparticle describing two rotating molecules that align with respect to each other due to the effective attractive interaction mediated by the excitations of the bath. We study this system in different parameter regimes and apply several theoretical approaches to describe its properties. Using a Born–Oppenheimer approximation, we investigate the dependence of the effective intermolecular interaction on the rotational state of the two molecules. In the strong-coupling regime, a product-state ansatz shows that the molecules tend to have a strong alignment in the ground state. To investigate the system in the weak-coupling regime, we apply a one-phonon excitation variational ansatz, which allows us to access the energy spectrum. In comparison to the angulon quasiparticle, the biangulon shows shifted angulon instabilities and an additional spectral instability, where resonant angular momentum transfer between the molecules and the bath takes place. These features are proposed as an experimentally observable signature for the formation of the biangulon quasiparticle. Finally, by using products of single angulon and bare impurity wave functions as basis states, we introduce a diagonalization scheme that allows us to describe the transition from two separated angulons to a biangulon as a function of the distance between the two molecules.}, author = {Li, Xiang and Yakaboylu, Enderalp and Bighin, Giacomo and Schmidt, Richard and Lemeshko, Mikhail and Deuchert, Andreas}, issn = {1089-7690}, journal = {The Journal of Chemical Physics}, keywords = {Physical and Theoretical Chemistry, General Physics and Astronomy}, number = {16}, publisher = {AIP Publishing}, title = {{Intermolecular forces and correlations mediated by a phonon bath}}, doi = {10.1063/1.5144759}, volume = {152}, year = {2020}, } @article{7422, abstract = {Biochemical reactions often occur at low copy numbers but at once in crowded and diverse environments. Space and stochasticity therefore play an essential role in biochemical networks. Spatial-stochastic simulations have become a prominent tool for understanding how stochasticity at the microscopic level influences the macroscopic behavior of such systems. While particle-based models guarantee the level of detail necessary to accurately describe the microscopic dynamics at very low copy numbers, the algorithms used to simulate them typically imply trade-offs between computational efficiency and biochemical accuracy. eGFRD (enhanced Green’s Function Reaction Dynamics) is an exact algorithm that evades such trade-offs by partitioning the N-particle system into M ≤ N analytically tractable one- and two-particle systems; the analytical solutions (Green’s functions) then are used to implement an event-driven particle-based scheme that allows particles to make large jumps in time and space while retaining access to their state variables at arbitrary simulation times. Here we present “eGFRD2,” a new eGFRD version that implements the principle of eGFRD in all dimensions, thus enabling efficient particle-based simulation of biochemical reaction-diffusion processes in the 3D cytoplasm, on 2D planes representing membranes, and on 1D elongated cylinders representative of, e.g., cytoskeletal tracks or DNA; in 1D, it also incorporates convective motion used to model active transport. We find that, for low particle densities, eGFRD2 is up to 6 orders of magnitude faster than conventional Brownian dynamics. We exemplify the capabilities of eGFRD2 by simulating an idealized model of Pom1 gradient formation, which involves 3D diffusion, active transport on microtubules, and autophosphorylation on the membrane, confirming recent experimental and theoretical results on this system to hold under genuinely stochastic conditions.}, author = {Sokolowski, Thomas R and Paijmans, Joris and Bossen, Laurens and Miedema, Thomas and Wehrens, Martijn and Becker, Nils B. and Kaizu, Kazunari and Takahashi, Koichi and Dogterom, Marileen and ten Wolde, Pieter Rein}, issn = {1089-7690}, journal = {The Journal of Chemical Physics}, number = {5}, publisher = {AIP Publishing}, title = {{eGFRD in all dimensions}}, doi = {10.1063/1.5064867}, volume = {150}, year = {2019}, } @article{9659, abstract = {The curvature dependence of interfacial free energy, which is crucial in quantitatively predicting nucleation kinetics and the stability of bubbles and droplets, is quantified by the Tolman length δ. For solid-liquid interfaces, however, δ has never been computed directly due to various theoretical and practical challenges. Here we perform a direct evaluation of the Tolman length from atomistic simulations of a solid-liquid planar interface in out-of-equilibrium conditions, by first computing the surface tension from the amplitude of thermal capillary fluctuations of a localized version of the Gibbs dividing surface and by then calculating how much the surface energy changes when it is defined relative to the equimolar dividing surface. We computed δ for a model potential, and found a good agreement with the values indirectly inferred from nucleation simulations. The agreement not only validates our approach but also suggests that the nucleation free energy of the system can be perfectly described using classical nucleation theory if the Tolman length is taken into account.}, author = {Cheng, Bingqing and Ceriotti, Michele}, issn = {1089-7690}, journal = {The Journal of Chemical Physics}, number = {23}, publisher = {AIP Publishing}, title = {{Communication: Computing the Tolman length for solid-liquid interfaces}}, doi = {10.1063/1.5038396}, volume = {148}, year = {2018}, } @article{14006, abstract = {We present a theoretical formalism for the calculation of attosecond delays in molecular photoionization. It is shown how delays relevant to one-photon-ionization, also known as Eisenbud-Wigner-Smith delays, can be obtained from the complex dipole matrix elements provided by molecular quantum scattering theory. These results are used to derive formulae for the delays measured by two-photon attosecond interferometry based on an attosecond pulse train and a dressing femtosecond infrared pulse. These effective delays are first expressed in the molecular frame where maximal information about the molecular photoionization dynamics is available. The effects of averaging over the emission direction of the electron and the molecular orientation are introduced analytically. We illustrate this general formalism for the case of two polyatomic molecules. N2O serves as an example of a polar linear molecule characterized by complex photoionization dynamics resulting from the presence of molecular shape resonances. H2O illustrates the case of a non-linear molecule with comparably simple photoionization dynamics resulting from a flat continuum. Our theory establishes the foundation for interpreting measurements of the photoionization dynamics of all molecules by attosecond metrology.}, author = {Baykusheva, Denitsa Rangelova and Wörner, Hans Jakob}, issn = {1089-7690}, journal = {The Journal of Chemical Physics}, keywords = {Physical and Theoretical Chemistry, General Physics and Astronomy}, number = {12}, publisher = {AIP Publishing}, title = {{Theory of attosecond delays in molecular photoionization}}, doi = {10.1063/1.4977933}, volume = {146}, year = {2017}, }