Source author record

Pietro Faccioli

Pietro Faccioli 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

25works
11topics
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

25 published item(s)

preprint2022arXiv

RNA folding landscapes from explicit solvent all-atom simulations

Atomically detailed simulations of RNA folding have proven very challenging in view of the difficulties of developing realistic force fields and the intrinsic computational complexity of sampling rare conformational transitions. To tackle both these issues, we extend to RNA an enhanced path sampling method previously successfully applied to proteins. In this scheme, the information about the RNA's native structure is harnessed by a soft history-dependent biasing force, which is added to the atomistic force field, thus promoting the generation of productive folding trajectories. Here, we report on the results of simulations in explicit solvent of RNA molecules from 20 to 47 nucleotides long and increasing topological complexity. From a statistical analysis of the folding pathways we infer that, differently from proteins, the underlying free energy landscape is significantly frustrated, even for relatively small chains with a simple topology. The folding mechanism and the thermodynamics are in agreement with the available experiments and some of the existing coarse-grained models. This scheme provides a fully microscopic characterization of RNA folding, relating the kinetics and dynamics of the transition to the chemistry of the chain and its solvent. Therefore, it provides a transferable framework that sets the stage for future translational applications.

preprint2021arXiv

Antimalarial Artefenomel Inhibits Human SARS-CoV-2 Replication in Cells while Suppressing the Receptor ACE2

The steep climbing of victims caused by the new coronavirus disease 2019 (COVID-19) throughout the planet is sparking an unprecedented effort to identify effective therapeutic regimens to tackle the pandemic. The SARS-CoV-2 virus is known to gain entry into various cell types through the binding of one of its surface proteins (spike) to the host Angiotensin-Converting Enzyme 2 (ACE2). Thus, spike-ACE2 interaction represents a major target for vaccines and antiviral drugs. A novel method has been recently described by some of the authors to pharmacologically downregulate the expression of target proteins at the post-translational level. This technology builds on computational advancements in the simulation of folding mechanisms to rationally block protein expression by targeting folding intermediates, hence hampering the folding process. Here, we report the all-atom simulations of the entire sequence of events underlying the folding pathway of ACE2. Our data revealed the existence of a folding intermediate showing two druggable pockets hidden in the native conformation. Both pockets were targeted by a virtual screening repurposing campaign aimed at quickly identifying drugs capable to decrease the expression of ACE2. We identified four compounds capable of lowering ACE2 expression in Vero cells in a dose-dependent fashion. All these molecules were found to inhibit the entry into cells of a pseudotyped retrovirus exposing the SARS-CoV-2 spike protein. Importantly, the antiviral activity has been tested against live SARS-CoV-2 (MEX-BC2/2020 strain). One of the selected drugs (Artefenomel) could completely prevent cytopathic effects induced by the presence of the virus, thus showing antiviral activity against SARS-CoV-2. Ongoing studies are further evaluating the possibility of repurposing these drugs for the treatment of COVID-19.

preprint2021arXiv

Dominant Reaction Pathways by Quantum Computing

Characterizing thermally activated transitions in high-dimensional rugged energy surfaces is a very challenging task for classical computers. Here, we develop a quantum annealing scheme to solve this problem. First, the task of finding the most probable transition paths in configuration space is reduced to a shortest-path problem defined on a suitable weighted graph. Next, this optimization problem is mapped into finding the ground state of a generalized Ising model. A finite-size scaling analysis suggests this task may be solvable efficiently by a quantum annealing machine. Our approach leverages on the quantized nature of qubits to describe transitions between different system's configurations. Since it does not involve any lattice space discretization, it paves the way towards future biophysical applications of quantum computing based on realistic all-atom models.

preprint2021arXiv

Stochastic Dynamics and Bound States of Heavy Impurities in a Fermi Bath

We investigate the dynamics of heavy impurities embedded in an ultra-cold Fermi gas by using a Generalized Langevin equation. The latter -- derived by means of influence functional theory -- describes the stochastic classical dynamics of the impurities and the quantum nature of the fermionic bath manifests in the emergent interaction between the impurities and in the viscosity tensor. By focusing on the two-impurity case, we predict the existence of bound states, in different conditions of coupling and temperature, and whose life-time can be analytically estimated. Our predictions should be testable using cold-gases platforms within current technology.

preprint2021arXiv

Target search of active agents crossing high energy barriers

Target search by active agents in rugged energy landscapes has remained a challenge because standard enhanced sampling methods do not apply to irreversible dynamics. We overcome this non-equilibrium rare-event problem by developing an algorithm generalizing transition-path sampling to active Brownian dynamics. This method is exemplified and benchmarked for a paradigmatic two-dimensional potential with a high barrier. We find that even in such a simple landscape the structure and kinetics of the ensemble of transition paths change drastically in the presence of activity. Indeed, active Brownian particles reach the target more frequently than passive Brownian particles, following longer and counterintuitive search patterns.

preprint2020arXiv

From prompt to direct J/$ψ$ production: new insights on the $χ_{c1}$ and $χ_{c2}$ polarizations and feed-down contributions from a global-fit analysis of mid-rapidity LHC data

While the prompt J/$ψ$ cross section and polarization have been measured with good precision as a function of transverse momentum, $p_{\rm T}$, those of the directly produced J/$ψ$ are practically unknown, given that the cross sections and polarizations of the $χ_{c1}$ and $χ_{c2}$ mesons, large indirect contributors to J/$ψ$ production, are only known with rather poor accuracy. The lack of precise measurements of the $χ_{cJ}$ polarizations induces large uncertainties in the level of their feed-down contributions to the prompt J/$ψ$ yield, because of the polarization-dependent acceptance corrections. The experimental panorama of charmonium production can be significantly improved through a consistent and model-independent global analysis of existing measurements of J/$ψ$, $ψ$(2S) and $χ_{c}$ cross sections and polarizations, faithfully respecting all the correlations and uncertainties. In particular, it is seen that the $χ_{cJ}$ polarizations and feed-down fractions to J/$ψ$ production have a negligible dependence on the J/$ψ$ $p_{\rm T}$, with average values $λ_\vartheta^{χ_{c1}} = 0.55 \pm 0.23$, $λ_\vartheta^{χ_{c2}} = -0.39 \pm 0.22$, $R^{χ_{c1}} = (18.8 \pm 1.4)\%$ and $R^{χ_{c2}} = (6.5 \pm 0.5)\%$. The analysis also shows that $(67.2 \pm 1.9)\%$ of the prompt J/$ψ$ yield is due to directly-produced mesons, of polarization constrained to remarkably small values, $λ_\vartheta^{{\rm J}/ψ} = 0.04 \pm 0.06$.

preprint2016arXiv

Dimensional Reduction of Markov State Models from Renormalization Group Theory

Renormalization Group (RG) theory provides the theoretical framework to define Effective Theories (ETs), i.e. systematic low-resolution approximations of arbitrary microscopic models. Markov State Models (MSMs) are shown to be rigorous ETs for Molecular Dynamics (MD). Based on this fact, we use Real Space RG to vary the resolution of a MSM and define an algorithm for clustering microstates into macrostates. The result is a lower dimensional stochastic model which, by construction, provides the optimal coarse-grained Markovian representation of the system's relaxation kinetics. To illustrate and validate our theory, we analyze a number of test systems of increasing complexity, ranging from synthetic toy models to two realistic applications, built form all-atom MD simulations. The computational cost of computing the low-dimensional model remains affordable on a desktop computer even for thousands of microstates.

preprint2015arXiv

Heavy quark bound states in a quark-gluon plasma: dissociation and recombination

We present a comprehensive approach to the dynamics of heavy quarks in a quark gluon plasma, including the possibility of bound state formation and dissociation. In this exploratory paper, we restrict ourselves to the case of an Abelian plasma, but the extension of the techniques used to the non Abelian case is straightforward. A chain of well defined approximations leads eventually to a generalized Langevin equation, where the force and the noise terms are determined from a correlation function of the equilibrium plasma, and depend explicitly on the configuration of the heavy quarks. We solve the Langevin equation for various initial conditions, various numbers of heavy quark-antiquark pairs, and various temperatures of the plasma. Results of simulations illustrate various expected phenomena: dissociation of bound states as a result of combined effects of screening of the potential and collisions with the plasma constituent, formation of bound pairs (recombination) that occurs when enough heavy quarks are present in the system.

preprint2014arXiv

Quarkonium production in the LHC era: a polarized perspective

Polarization measurements are usually considered as the most difficult challenge for the QCD description of quarkonium production. In fact, global data fits for the determination of the non-perturbative parameters of bound-state formation traditionally exclude polarization observables and use them as a posteriori verifications of the predictions, with perplexing results. With a change of perspective, we move polarization data to the centre of the study, advocating that they actually provide the strongest fundamental indications about the production mechanisms, even before we explicitly consider perturbative calculations. Considering psi(2S) and Y(3S) measurements from LHC experiments and state-of-the-art NLO short-distance calculations in the framework of non-relativistic QCD factorization (NRQCD), we perform a search for a kinematic domain where the polarizations can be correctly reproduced together with the cross sections, by systematically scanning the phase space and accurately treating the experimental uncertainties. This strategy provides a straightforward solution to the "quarkonium polarization puzzle" and reassuring signs that the theoretical framework is reliable. At the same time, the results expose unexpected hierarchies in the non-perturbative NRQCD parameters, that open new paths towards the understanding of bound-state formation in QCD.

preprint2013arXiv

Folding Pathways of a Knotted Protein with a Realistic Atomistic Force Field

We report on atomistic simulation of the folding of a natively-knotted protein, MJ0366, based on a realistic force field. To the best of our knowledge this is the first reported effort where a realistic force field is used to investigate the folding pathways of a protein with complex native topology. By using the dominant-reaction pathway scheme we collected about 30 successful folding trajectories for the 82-amino acid long trefoil-knotted protein. Despite the dissimilarity of their initial unfolded configuration, these trajectories reach the natively-knotted state through a remarkably similar succession of steps. In particular it is found that knotting occurs essentially through a threading mechanism, involving the passage of the C-terminal through an open region created by the formation of the native beta-sheet at an earlier stage. The dominance of the knotting by threading mechanism is not observed in MJ0366 folding simulations using simplified, native-centric models. This points to a previously underappreciated role of concerted amino acid interactions, including non-native ones, in aiding the appropriate order of contact formation to achieve knotting.

preprint2013arXiv

Minimal physical constraints on the angular distributions of two-body boson decays

The angular distribution of the two-body decay of a boson of unknown properties is strongly constrained by angular momentum conservation and rotation invariance, as well as by the nature of the detected decay particles and of the colliding ones. Knowing the border between the "physical" and "unphysical" parameter domains defined by these "minimal constraints" (excluding specific hypotheses on what is still subject of measurement) is a useful ingredient in the experimental determinations of angular distributions and can provide model-independent criteria for spin characterizations. In particular, analysing the angular decay distribution with the general parametrization for the J = 2 case can provide a model-independent discrimination between the J = 0 and J = 2 hypotheses for a particle produced by two real gluons and decaying into two real photons.

preprint2013arXiv

QCD Topology at Finite Temperature: Statistical Mechanics of Selfdual Dyons

Topological phenomena in gauge theories have long been recognized as the driving force for chiral symmetry breaking and confinement. These phenomena can be conveniently investigated in the semi-classical picture, in which the topological charge is entirely carried by (anti-)self-dual gauge configurations. In such an approach, it has been shown that near the critical temperature, the non-zero expectation value of the Polyakov loop (holonomy) triggers the "Higgsing" of the color group, generating the splitting of instantons into $N_c$ self-dual dyons. A number of lattice simulations have provided some evidence for such dyons, and traced their relation with specific observables, such as the Dirac eigenvalue spectrum. In this work, we formulate a model, based on one-loop partition function and including Coulomb interaction, screening and fermion zee modes. We then perform the first numerical Monte Carlo simulations of a statistical ensemble of self-dual dyons,as a function of their density, quark mass and the number of flavors. We study different dyonic two-point correlation functions and we compute the Dirac spectrum, as a function of the ensemble diluteness and the number of quark flavors.

preprint2012arXiv

Observation of chi_c and chi_b nuclear suppression via dilepton polarization measurements

We demonstrate that it is possible to use the polarization of vector quarkonia, measured from dilepton event samples, as an instrument to study the suppression of chi_c and chi_b in heavy-ion collisions, where a direct determination of signal yields involving the identification of low-energy photons is essentially impossible. A change of the observed J/psi and Upsilon(1S) polarizations from proton-proton to central nucleus-nucleus collisions would directly reflect differences in the nuclear dissociation patterns of S and P states and may provide a strong indication for quarkonium sequential suppression in the quark-gluon plasma.

preprint2012arXiv

Questions and prospects in quarkonium polarization measurements from proton-proton to nucleus-nucleus collisions

Polarization measurements are the best instrument to understand how quark and antiquark combine into the different quarkonium states, but no model has so far succeeded in explaining the measured J/psi and Upsilon polarizations. On the other hand, the experimental data in proton-antiproton and proton-nucleus collisions are inconsistent, incomplete and ambiguous. New analyses will have to properly address often underestimated issues: the existence of azimuthal anisotropies, the dependence on the reference frame, the influence of the experimental acceptance on the comparison with other measurements and with theory. Additionally, a recently developed frame-invariant formalism will provide an alternative and often more immediate physical viewpoint and, at the same time, will help probing systematic effects due to experimental biases. The role of feed-down decays from heavier states, a crucial missing piece in the current experimental knowledge, will have to be investigated. Ultimately, quarkonium polarization measurements will also offer new possibilities in the study of the properties of the quark-gluon plasma.

preprint2011arXiv

Computing the Effective Hamiltonian of Low-Energy Vacuum Gauge Fields

A standard approach to investigate the non-perturbative QCD dynamics is through vacuum models which emphasize the role played by specific gauge field fluctuations, such as instantons, monopoles or vortexes. The effective Hamiltonian describing the dynamics of the low-energy degrees of freedom in such approaches is usually postulated phenomenologically, or obtained through uncontrolled approximations. In a recent paper, we have shown how lattice field theory simulations can be used to rigorously compute the effective Hamiltonian of arbitrary vacuum models by stochastically performing the path integral over all the vacuum field fluctuations which are not explicitly taken into account. In this work, we present the first illustrative application of such an approach to a gauge theory and we use it to compute the instanton size distribution in SU(2) gluon-dynamics in a fully model independent and parameter-free way.

preprint2011arXiv

Determination of chi_c and chi_b polarizations from dilepton angular distributions in radiative decays

The angular distributions of the decay products in the successive decays chi_c (chi_b) to J/psi (Upsilon) gamma and J/psi (Upsilon) to l+l- are calculated as a function of the angular momentum composition of the decaying chi meson and of the multipole structure of the photon radiation, using a formalism independent of production mechanisms and polarization frames. The polarizations of the chi states produced in high energy collisions can be derived from the dilepton decay distributions of the daughter J/psi or Upsilon mesons, with a reduced dependence on the details of the photon reconstruction or simulation. Moreover, this method eliminates the dependence of the polarization measurement on the actual details of the multipole structure of the radiative transition. Problematic points in previous calculations of the chi_c decay angular distributions are identified and clarified.

preprint2011arXiv

Investigating Biological Matter with Theoretical Nuclear Physics Methods

The internal dynamics of strongly interacting systems and that of biomolecules such as proteins display several important analogies, despite the huge difference in their characteristic energy and length scales. For example, in all such systems, collective excitations, cooperative transitions and phase transitions emerge as the result of the interplay of strong correlations with quantum or thermal fluctuations. In view of such an observation, some theoretical methods initially developed in the context of theoretical nuclear physics have been adapted to investigate the dynamics of biomolecules. In this talk, we review some of our recent studies performed along this direction. In particular, we discuss how the path integral formulation of the molecular dynamics allows to overcome some of the long-standing problems and limitations which emerge when simulating the protein folding dynamics at the atomistic level of detail.

preprint2011arXiv

Model-independent constraints on the shape parameters of dilepton angular distributions

The coefficients determining the dilepton decay angular distribution of vector particles obey certain positivity constraints and a rotation-invariant identity. These relations are a direct consequence of the covariance properties of angular momentum eigenstates and are independent of the production mechanism. The Lam-Tung relation can be derived as a particular case, simply recognizing that the Drell-Yan dilepton is always produced transversely polarized with respect to one or more quantization axes. The dilepton angular distribution continues to be characterized by a frame-independent identity also when the Lam-Tung relation is violated. Moreover, the violation can be easily characterized by measuring a one-dimensional distribution depending on one shape coefficient.

preprint2010arXiv

A new approach to quarkonium polarization studies

Significant progress in understanding quarkonium production requires improved polarization measurements, fully considering the intrinsic multidimensionality of the problem. We propose a frame-invariant formalism which minimizes the dependence of the measured result on the experimental acceptance, facilitates the comparison with theoretical calculations, and provides a much needed control over systematic effects due to detector limitations and analysis biases. This formalism is a direct and generic consequence of the rotational invariance of the dilepton decay distribution and is independent of any assumptions specific to particular models of quarkonium production.

preprint2010arXiv

Molecular Dynamics at Low Time Resolution

The internal dynamics of macro-molecular systems is characterized by widely separated time scales, ranging from fraction of ps to ns. In ordinary molecular dynamics simulations, the elementary time step dt used to integrate the equation of motion needs to be chosen much smaller of the shortest time scale, in order not to cut-off important physical effects. We show that, in systems obeying the over-damped Langevin Eq., the fast molecular dynamics which occurs at time scales smaller than dt can be analytically integrated out and gives raise to a time-dependent correction to the diffusion coefficient, which we rigorously compute. The resulting effective Langevin equation describes by construction the same long-time dynamics, but has a lower time resolution power, hence it can be integrated using larger time steps dt. We illustrate and validate this method by studying the diffusion of a point-particle in a one-dimensional toy-model and the denaturation of a protein.

preprint2010arXiv

Rotation-invariant observables in parity-violating decays of vector particles to fermion pairs

The di-fermion angular distribution observed in decays of inclusively produced vector particles is characterized by two frame-independent observables, reflecting the average spin-alignment of the produced particle and the magnitude of parity violation in the decay. The existence of these observables derives from the rotational properties of angular momentum eigenstates and is a completely general result, valid for any J=1 state and independent of the production process. Rotation-invariant formulations of polarization and of the decay parity-asymmetry can provide more significant measurements than the commonly used frame-dependent definitions, also improving the quality of the comparisons between the measurements and the theoretical calculations.

preprint2010arXiv

Rotation-invariant relations in vector meson decays into fermion pairs

The rotational properties of angular momentum eigenstates imply the existence of a frame-independent relation among the parameters of the decay distribution of vector mesons into fermions. This relation is a generalization of the Lam-Tung identity, a result specific to Drell-Yan production in perturbative QCD, here shown to be equivalent to the dynamical condition that the dilepton always originates from a transversely polarized photon.

preprint2010arXiv

Towards the experimental clarification of quarkonium polarization

We highlight issues which are often underestimated in the experimental analyses on quarkonium polarization: the relation between the parameters of the angular distributions and the angular momentum composition of the quarkonium, the importance of the choice of the reference frame, the interplay between observed decay and production kinematics, and the consequent influence of the experimental acceptance on the comparison between experimental measurements and theoretical calculations. Given the puzzles raised by the available experimental results, new measurements must provide more detailed information, such that physical conclusions can be derived without relying on model-dependent assumptions. We describe a frame-invariant formalism which minimizes the dependence of the measurements on the experimental acceptance, facilitates the comparison with theoretical calculations, and probes systematic effects due to experimental biases. This formalism is a direct and generic consequence of the rotational invariance of the dilepton decay distribution and is independent of any assumptions specific to particular models of quarkonium production. The use of this improved approach, which exploits the intrinsic multidimensionality of the problem, will significantly contribute to a faster progress in our understanding of quarkonium production, especially if adopted as a common analysis framework by the LHC experiments, which will soon perform analyses of quarkonium polarization in proton-proton collisions.

preprint2000arXiv

CKM matrix: the `over-consistent' picture of the unitarity triangle

In presenting an up-to-date account of the experimental knowledge of the CKM matrix, special emphasis is placed on the exceptional degree of consistency shown by the current Standard Model determination of the unitarity triangle; some implications in the question of how the theoretical nature of the dominant uncertainties affects the Standard Model predictions are discussed. Recent experimental results include improved determinations of |V_{cb}| by OPAL and CLEO, the confirmation of rare (charmless hadronic and electromagnetic penguin) B decays and preliminary sin(2 beta) measurements by BaBar and Belle (new world average: sin(2 beta) = 0.48 +0.22 -0.24). The updated constraints lead to the Standard Model predictions sin(2 beta) = 0.68 +/- 0.03(exp) +/- 0.04(th), sin(2 alpha) = -0.43 +/- 0.15(exp) +/- 0.20(th), gamma = 56 +/- 5(exp)(+5 -6)(th) deg, Delta m_{Bs} = 16.2 (+2.7 -0.3)(exp)(+1.5 -1.0)(th) ps^{-1}.

preprint1998arXiv

Nucleon Structure Functions and Light-Front Dynamics

We present a quark-parton model to describe polarized and unpolarized nucleon structure functions. The twist-two matrix elements for the QCD evolution analysis of lepton-hadron scattering are calculated within a light-front covariant quark model. The relativistic effects in the three-body wave function are discussed for both the polarized and unpolarized cases. Predictions are given for the polarized gluon distributions as will be seen in future experiments.