Source author record

Vadim A. Markel

Vadim A. Markel 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

12works
7topics
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

12 published item(s)

preprint2016arXiv

Can photonic crystals be homogenized in higher bands?

We consider conditions under which photonic crystals (PCs) can be homogenized in the higher photonic bands and, in particular, near the $Γ$-point. By homogenization we mean introducing some effective local parameters $ε_{\rm eff}$ and $μ_{\rm eff}$ that describe reflection, refraction and propagation of electromagnetic waves in the PC adequately. The parameters $ε_{\rm eff}$ and $μ_{\rm eff}$ can be associated with a hypothetical homogeneous effective medium. In particular, if the PC is homogenizable, the dispersion relations and isofrequency lines in the effective medium and in the PC should coincide to some level of approximation. We can view this requirement as a necessary condition of homogenizability. In the vicinity of a $Γ$-point, real isofrequency lines of two-dimensional PCs can be close to mathematical circles, just like in the case of isotropic homogeneous materials. Thus, one may be tempted to conclude that introduction of an effective medium is possible and, at least, the necessary condition of homogenizability holds in this case. We, however, show that this conclusion is incorrect: complex dispersion points must be included into consideration even in the case of strictly non-absorbing materials. By analyzing the complex dispersion relations and the corresponding isofrequency lines, we have found that two-dimensional PCs with $C_4$ and $C_6$ symmetries are not homogenizable in the higher photonic bands. We also draw a distinction between spurious $Γ$-point frequencies that are due to Brillouin-zone folding of Bloch bands and "true" $Γ$-point frequencies that are due to multiple scattering. Understanding of the physically different phenomena that lead to the appearance of spurious and "true" $Γ$-point frequencies is important for the theory of homogenization.

preprint2016arXiv

Nonlinear inverse problem by T-matrix completion. I. Theory

We propose a conceptually new method for solving nonlinear inverse scattering problems (ISPs) such as are commonly encountered in tomographic ultrasound imaging, seismology and other applications. The method is inspired by the theory of nonlocality of physical interactions and utilizes the relevant formalism. We formulate the ISP as a problem whose goal is to determine an unknown interaction potential $V$ from external scattering data. Although we seek a local (diagonally-dominated) $V$ as the solution to the posed problem, we allow $V$ to be nonlocal at the intermediate stages of iterations. This allows us to utilize the one-to-one correspondence between $V$ and the T-matrix of the problem, $T$. Here it is important to realize that not every $T$ corresponds to a diagonal $V$ and we, therefore, relax the usual condition of strict diagonality (locality) of $V$. An iterative algorithm is proposed in which we seek $T$ that is (i) compatible with the measured scattering data and (ii) corresponds to an interaction potential $V$ that is as diagonally-dominated as possible. We refer to this algorithm as to the data-compatible T-matrix completion (DCTMC). This paper is Part I in a two-part series and contains theory only. Numerical examples of image reconstruction in a strongly nonlinear regime are given in Part II. The method described in this paper is particularly well suited for very large data sets that become increasingly available with the use of modern measurement techniques and instrumentation.

preprint2016arXiv

Solution of the inverse scattering problem by T-matrix completion. II. Simulations

This is Part II of the paper series on data-compatible T-matrix completion (DCTMC), which is a method for solving nonlinear inverse problems. Part I of the series contains theory and here we present simulations for inverse scattering of scalar waves. The underlying mathematical model is the scalar wave equation and the object function that is reconstructed is the medium susceptibility. The simulations are relevant to ultrasound tomographic imaging and seismic tomography. It is shown that DCTMC is a viable method for solving strongly nonlinear inverse problems with large data sets. It provides not only the overall shape of the object but the quantitative contrast, which can correspond, for instance, to the variable speed of sound in the imaged medium.

preprint2015arXiv

Nonasymptotic Homogenization of Periodic Electromagnetic Structures: Uncertainty Principles

We show that artificial magnetism of periodic dielectric or metal/dielectric structures has limitations and is subject to at least two "uncertainty principles". First, the stronger the magnetic response (the deviation of the effective permeability tensor from identity), the less accurate ("certain") the predictions of any homogeneous model. Second, if the magnetic response is strong, then homogenization cannot accurately reproduce the transmission and reflection parameters and, simultaneously, power dissipation in the material. These principles are general and not confined to any particular method of homogenization. Our theoretical analysis is supplemented with a numerical example: a hexahedral lattice of cylindrical air holes in a dielectric host. Even though this case is highly isotropic, which might be thought as conducive to homogenization, the uncertainty principles remain valid.

preprint2014arXiv

Inversion of the star transform

We define the star transform as a generalization of the broken ray transform introduced by us in previous work. The advantages of using the star transform include the possibility to reconstruct the absorption and the scattering coefficients of the medium separately and simultaneously (from the same data) and the possibility to utilize scattered radiation which, in the case of the conventional X-ray tomography, is discarded. In this paper, we derive the star transform from physical principles, discuss its mathematical properties and analyze numerical stability of inversion. In particular, it is shown that stable inversion of the star transform can be obtained only for configurations involving odd number of rays. Several computationally-efficient inversion algorithms are derived and tested numerically.

preprint2013arXiv

Theoretical and numerical analysis of current-driven electromagnetic homogenization and the problem of effective medium parameters for finite samples

Reflection and refraction of electromagnetic waves by artificial periodic composites (metamaterials) can be accurately modeled by an effective medium theory only if the boundary of the medium is explicitly taken into account and the two effective parameters of the medium -- the index of refraction and the impedance -- are correctly determined. Theories that consider infinite periodic composites do not satisfy the above condition. As a result, they cannot model reflection and transmission by finite samples with the desired accuracy and are not useful for design of metamaterial-based devices. As an instructive case in point, we consider the "current-driven" homogenization theory, which has recently gained popularity. We apply this theory to the case of one-dimensional periodic medium wherein both exact and homogenization results can be obtained analytically in closed form. We show that, beyond the well-understood zero-cell limit, the current-driven homogenization result is inconsistent with the exact reflection and transmission characteristics of the slab.

preprint2012arXiv

Homogenization of Maxwell's equations in periodic composites

We consider the problem of homogenizing the Maxwell equations for periodic composites. The analysis is based on Bloch-Floquet theory. We calculate explicitly the reflection coefficient for a half-space, and derive and implement a computationally-efficient continued-fraction expansion for the effective permittivity. Our results are illustrated by numerical computations for the case of two-dimensional systems. The homogenization theory of this paper is designed to predict various physically-measurable quantities rather than to simply approximate certain coefficients in a PDE.

preprint2011arXiv

Quantum theory of electric polarization nonlinearity in metal nanofilms

We develop a quantum theory of electron confinement in metal nanofilms. The theory is used to compute the nonlinear response of the film to a static or low-frequency external electric field and to investigate the role of boundary conditions imposed on the metal surface. We find that the sign and magnitude of the nonlinear polarizability depends dramatically on the type of boundary condition used.

preprint2011arXiv

Theoretical and numerical investigation of the size-dependent optical effects in metal nanoparticles

We further develop the theory of quantum finite-size effects in metallic nanoparticles, which was originally formulated by Hache, Ricard and Flytzanis [J. Opt. Soc. Am. B 3, 1647 (1986)] and (in a somewhat corrected form) by Rautian [Sov. Phys. JETP 85, 451 (1997)]. These references consider a metal nanoparticle as a degenerate Fermi gas of conduction electrons in an infinitely-high spherical potential well. This model (referred to as the HRFR model below) yields mathematical expressions for the linear and the third-order nonlinear polarizabilities of a nanoparticle in terms of infinite nested series. These series have not been evaluated numerically so far and, in the case of nonlinear polarizability, they can not be evaluated with the use of conventional computers due to the high computational complexity involved. Rautian has derived a set of remarkable analytical approximations to the series but direct numerical verification of Rautian's approximate formulas remained a formidable challenge. In this work, we derive an expression for the third-order nonlinear polarizability, which is exact within the HRFR model but amenable to numerical implementation. We then evaluate the expressions obtained by us numerically for both linear and nonlinear polarizabilities. We investigate the limits of applicability of Rautian's approximations and find that they are surprizingly accurate in a wide range of physical parameters. We also discuss the limits of small frequencies (comparable or below the Drude relaxation constant) and of large particle sizes (the bulk limit) and show that these limits are problematic for the HRFR model, irrespectively of any additional approximations used. Finally, we compare the HRFR model to the purely classical theory of nonlinear polarization of metal nanoparticles developed by us earlier [Phys. Rev. Lett. 100, 47402 (2008)].

preprint2010arXiv

Comment on "Poynting vector, heating rate, and stored energy in structured materials: A first principles derivation"

In this Comment I argue that Silveirinha's criticism of my earlier work [M.G.Silveirinha, Phys.Rev.B 80, 235120 2009)] is based on an unphysical excitation model which involves an "external current" which overlaps with a continuous medium but is not subject to constitutive relations. When this excitation model is replaced by the conventional and experimentally-relevant model of excitation by external electromagnetic fields, it can be easily shown that Silveirinha's formulas contain the very results he wanted to disprove. This and a few other misconceptions present in Silveirinha's paper are subject of this Comment.

preprint2010arXiv

Inversion formulas for the broken-ray Radon transform

We consider the inverse problem of the broken ray transform (sometimes also referred to as the V-line transform). Explicit image reconstruction formulas are derived and tested numerically. The obtained formulas are generalizations of the filtered backprojection formula of the conventional Radon transform. The advantages of the broken ray transform include the possibility to reconstruct the absorption and the scattering coefficients of the medium simultaneously and the possibility to utilize scattered radiation which, in the case of the conventional X-ray tomography, is typically discarded.

preprint2008arXiv

Single-Scattering Optical Tomography: Simultaneous Reconstruction of Scattering and Absorption

We demonstrate that simultaneous reconstruction of scattering and absorption of a mesoscopic system using angularly-resolved measurements of scattered light intensity is possible. Image reconstruction is realized based on the algebraic inversion of a generalized Radon transform relating the scattering and absorption coefficients of the medium to the measured light intensity and derived using the single-scattering approximation to the radiative transport equation.