Researcher profile

Asaf Farhi

Asaf Farhi contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
11works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

11 published item(s)

preprint2022arXiv

Purcell factors and Forster resonance energy transfer in proximity to helical structures

Both spontaneous emission and resonant energy transfer can be enhanced significantly when the emitter is placed in the vicinity of metallic or crystal structures. This enhancement can be described using the electromagnetic Green tensor and is determined by the dominant surface modes of the structure. Here we use the eigenpermittivity formalism to derive the spontaneous emission and FRET rates in the quasistatic regime in a two-constituent medium with an anisotropic inclusion. We then apply our results to a helical structure supporting synchronous vibrations and evaluate the contribution of these modes, which are associated with a strong and delocalized response. We show that this contribution can result in large Purcell factors and long-range FRET, which oscillates with the helix pitch. These findings may have implications in understanding and controlling the interactions of molecules close to helical structures such as the microtubules.

preprint2020arXiv

Coupling electrodynamic fields to vibrational modes in helical structures

Helical structures like alpha helices, DNA, and microtubules have profound importance in biology. It has been suggested that these periodic arrangements of constituent units could support collective excitations similarly to crystalline solids. Here, we examine the interaction between such constructs and oscillating dipoles, and evaluate the role of the helicity in the coupling between electrodynamic fields and vibrations. Based on a vibrational and eigenfunction analyses we discover a group of modes of coherent oscillations that give rise to a strong and delocalized response, selectivity in frequency, and typical interaction range. To describe the field scattering due to the structure vibrations we consider an anisotropic permittivity with a helical periodicity, which applies to all vibration types and close dipole locations. This new type of resonances identified here may help explain the role of electrodynamic fields in the diverse functionality of cytoskeletal microtubules in the cellular environment.

preprint2019arXiv

Three-dimensional-subwavelength field localization, time reversal of sources, and infinite-asymptotic degeneracy in spherical structures

High-resolution field localization in three dimensions is one of the main challenges in optics and has immense importance in fields such as chemistry, biology, and medicine. In order to generate the time reversed signal of a monochromatic source \emph{discrete} sources that are modulated according to the wave amplitude on a spherical envelope are required, rendering it applicable only in acoustics. Here we approach these challenges by introducing a spherical layer with a resonant permittivity, which naturally generates the spatially \emph{continuous} time-reversed signal of an atomic and molecular multipole transition at the origin. We start by utilizing a spherical layer with a resonant TM $l=1$ permittivity situated in a uniform medium to generate a free-space-subwavelength focal spot at the origin. We remove the degeneracy of the eigenfunctions of the composite medium by situating a point current source (or polarization) directed parallel to the spherical layer, which generates a focal spot at the origin \emph{independently} of its location. The free-space focal spot has a full width at half maximum of $0.4λ$ in the lateral axes and $0.58λ$ in the axial axis, which is tighter by a factor of $\sqrt{2}$ in each dimension in excitation-collection mode, overcoming the $λ/2$ far-field resolution limit in three dimensions. We then explore two directions to localize electric field with deep-subwavelength resolution in three dimensions using this setup. In addition, we show that spherical structures exhibit a new type of degeneracy in which an infinite number of eigenvalues asymptotically coalesce. This high degeneracy results in a variety of optical phenomena such as strong scattering and enhancement of absorption and emission from an atom or molecule by orders of magnitude compared with a standard resonance.

preprint2016arXiv

A novel method for calculating relative free energy of similar molecules in two environments

Calculating relative free energies is a topic of substantial interest and has many applications including solvation and binding free energies, which are used in computational drug discovery. However, there remain the challenges of accuracy, simple implementation, robustness and efficiency, which prevent the calculations from being automated and limit their use in computational drug discovery. Here we present an exact and complete decoupling analysis in which the partition functions of the compared systems decompose into the partition functions of the common and different subsystems. This decoupling analysis is applicable to submolecules with coupled degrees of freedom and to any potential function, enabling to remove less terms in the transformation. Then we show mathematically, in the context of partition function decoupling, that the two compared systems can be simulated separately, eliminating the need to design a composite system. We demonstrate the decoupling analysis and the separate transformations in an MD calculation of relative free energy for a general force field which is in agreement with experiments. We present a unified soft core technique that will ensure the monotonicity of the numerically integrated function. Finally, we show that when the systems have rugged energy landscape they can be equilibrated without introducing another sampling dimension. The concepts presented in the article have implications in accuracy, efficiency, simplicity and robustness of free energy calculations.

preprint2016arXiv

Eigenstates of the full Maxwell equations for a two-constituent composite medium and their application to a calculation of the local electric field of a time dependent point electric dipole in a flat-slabs microstructure

An exact calculation of the local electric field ${\bf E}({\bf r})$ is described for the case of a time dependent point electric dipole ${\bf p}e^{-iωt}$ in the top layer of an $ε_2$, $ε_1$, $ε_2$ three parallel slabs composite structure, where the $ε_1$ layer has a finite thickness $2d$ but the $ε_2$ layers are infinitely thick. For this purpose we first calculate all the eigenstates of the full Maxwell equations for the case where $μ=1$ everywhere in the system. The eigenvalues appear as special, non-physical values of $ε_1$ when $ε_2$ is given. These eigenstates are then used to develop an exact expansion for the physical values of ${\bf E}({\bf r})$ in the system characterized by physical values of $ε_1(ω)$ and $ε_2(ω)$. Results are compared with those of a previous calculation of the local field of a time dependent point charge in the quasi-static regime. Numerical results are shown for the local electric field in practically important configurations where attaining an optical image with sub-wavelength resolution has practical significance.

preprint2015arXiv

Calculation of molecular free energies in classical potentials

Free energies of molecules can be calculated by quantum computations or by normal mode classical calculations. However, the first can be computationally impractical for large molecules and the second is based on the assumption of harmonic dynamics. We present a novel, accurate and complete calculation of molecular free energies in standard classical potentials. In this method we transform the molecule by relaxing potential terms which depend on the coordinates of a group of atoms in that molecule and calculate the free energy difference associated with the transformation. Then, since the transformed molecule can be treated as non interacting systems, the free energy associated with these atoms is analytically or numerically calculated. This two-step calculation can be applied to calculate free energies of molecules or free energy difference between (possibly large) molecules in a general environment. We suggest the potential application of free energy calculation of chemical reactions in classical molecular simulations.

preprint2015arXiv

Free energy calculation of a molecule by removing VDW and Coulomb interactions in a transformation and treating the molecule as non interacting systems

Free energy calculations in molecular simulations have a variety of applications including determining the strength of molecular processes such as solvation and binding. It has been recently shown that when removing the VDW and Coulomb potential terms of a group of atoms in a molecule by performing a transformation, the molecule can be treated as non interacting systems in the free energy calculation. This treatment is applicable both when the molecule is in vacuum and in liquid and enables a very simple calculation of the free energies associated with the potentials that depend on the relative spherical coordinates of these atoms. Here we demonstrate the method in the free energy calculation of a Methanethiol molecule and compare the results to these obtained by MD simulations in vacuum and in water. The comparison shows agreement between the results and faster computation when using the method by factors varying between 5000 and 10^(12) for the same computational resources.

preprint2014arXiv

Exact analysis of a Veselago lens in the quasi-static regime

The resolution of conventional optical lenses is limited by the wavelength. Materials with negative refractive index have been shown to enable the generation of an enhanced resolution image where both propagating and non-propagating waves are employed. We analyze such a Veselago lens by exploiting some exact one dimensional integral expressions for the quasi-static electric potential of a point charge in that system. Those were recently obtained by expanding that potential in the quasi-static eigenfunctions of a three-flat-slabs composite structure. Numerical evaluations of those integrals, using realistic values for physical parameters like the electric permittivities of the constituent slabs and their thickness, reveal some surprising effects: E.g., the maximum concentration of the electric field occurs not at the geometric optics foci but at the interfaces between the negative permittivity slab and the positive permittivity slabs. The analysis provides simple computational guides for designing such structures in order to achieve enhanced resolution of an optical image.

preprint2013arXiv

A general, efficient and robust method to calculate free energy difference between systems

Calculating free energy differences is a topic of substantial interest and has many applications including molecular docking and hydration, solvation, and binding free energies which is used in computational drug discovery. However, in equilibrium methods the compared molecules are required to have large phase space overlap, which is usually not satisfied for two random systems, and there remain the challenges of robustness and automation. Here a highly efficient and robust method, that enables a wide range of comparisons, will be introduced, demonstrated and compared. In this method instead of transforming between one system into the other to perform the calculation each system is transformed into its replica with the different long range energy terms relaxed, which is inherently correlated with the original one, in order to eliminate the partition function difference arising from these terms. Then, since each transformed system can be treated as non interacting systems, the remaining difference in the (originally highly complex) partition function will cancel out.

preprint2013arXiv

Temperature Integration: an efficient procedure for calculation of free energy differences

We propose a method, Temperature Integration, which allows an efficient calculation of free energy differences between two systems of interest, with the same degrees of freedom, which may have rough energy landscapes. The method is based on calculating, for each single system, the difference between the values of lnZ at two temperatures, using a Parallel Tempering procedure. If our two systems of interest have the same phase space volume, they have the same values of lnZ at high-T, and we can obtain the free energy difference between them, using the two single-system calculations described above. If the phase space volume of a system is known, our method can be used to calculate its absolute (versus relative) free energy as well. We apply our method and demonstrate its efficiency on a toy model of hard rods on a 1-dimensional ring.

preprint2012arXiv

Comparing secondary structures of RNA and calculating the free energy of an interior loop using a novel method for calculating free energy

The thesis consists of two projects. In the first project, we present a software that analyses RNA secondary structures and compares them. The goal of this software is to find the differences between two secondary structures (experimental or predicted) in order to improve or compare algorithms for predicting secondary structures. Then, a comparison between secondary structures predicted by the Vienna package to those found experimentally is presented and cases in which there exists a difference between the prediction and the experimental structure are identified. As the differences originate mainly from faces and hydrogen bonds that are not allowed by the Vienna package, it is suggested that prediction may be improved by integrating them into the software. In the second project we calculate the free energy of an interior loop using Monte-Carlo simulation. We first present a semi-coarse grained model for interior loops of RNA, and the energy model for the different interactions. We then introduce the Monte-Carlo simulations and the method of Parallel Tempering which enables good sampling of configuration space by simulating a system simultaneously at several temperatures. Next we present Thermodynamic Integration, which is a method for calculating free energy differences. Then, We introduce a method that calculates the free energy significantly faster since we need to use only one parameter, $T$. To implement this method, we had to reach a regime in which the partition functions of the two systems are equal, which isn't satisfied for systems that have different entropies in the high temperature limit, so a solution to this problem had to be found. Free energy values calculated for various interior loops are shown, and may, if verified or with a more realistic modelling, be integrated into the Vienna package and supply an alternative to the experiments done.