Efficient Volumetric Method of Moments for Modeling Plasmonic Thin-Film Solar Cells with Periodic Structures

Metallic nanoparticles (NPs) support localized surface plasmon resonances (LSPRs), which enable to concentrate sunlight at the active layer of solar cells. However, full-wave modeling of the plasmonic solar cells faces great challenges in terms of huge computational workload and bad matrix condition. It is tremendously difficult to accurately and efficiently simulate near-field multiple scattering effects from plasmonic NPs embedded into solar cells. In this work, a preconditioned volume integral equation (VIE) is proposed to model plasmonic organic solar cells (OSCs). The diagonal block preconditioner is applied to different material domains of the device structure. As a result, better convergence and higher computing efficiency are achieved. Moreover, the calculation is further accelerated by two-dimensional periodic Green's functions. Using the proposed method, the dependences of optical absorption on the wavelengths and incident angles are investigated. Angular responses of the plasmonic OSCs show the super-Lambertian absorption on the plasmon resonance but near-Lambertian absorption off the plasmon resonance. The volumetric method of moments and explored physical understanding are of great help to investigate the optical responses of OSCs.Comment: 11 pages, 6 figure

Rigorous coherent-structure theory for falling liquid films: Viscous dispersion effects on bound-state formation and self-organization

We examine the interaction of two-dimensional solitary pulses on falling liquid films. We make use of the second-order model derived by Ruyer-Quil and Manneville [Eur. Phys. J. B 6, 277 (1998); Eur. Phys. J. B 15, 357 (2000); Phys. Fluids 14, 170 (2002)] by combining the long-wave approximation with a weighted residuals technique. The model includes (second-order) viscous dispersion effects which originate from the streamwise momentum equation and tangential stress balance. These effects play a dispersive role that primarily influences the shape of the capillary ripples in front of the solitary pulses. We show that different physical parameters, such as surface tension and viscosity, play a crucial role in the interaction between solitary pulses giving rise eventually to the formation of bound states consisting of two or more pulses separated by well-defined distances and travelling at the same velocity. By developing a rigorous coherent-structure theory, we are able to theoretically predict the pulse-separation distances for which bound states are formed. Viscous dispersion affects the distances at which bound states are observed. We show that the theory is in very good agreement with computations of the second-order model. We also demonstrate that the presence of bound states allows the film free surface to reach a self-organized state that can be statistically described in terms of a gas of solitary waves separated by a typical mean distance and characterized by a typical density

Controlled nanochannel lattice formation utilizing prepatterned substrates

Solid substrates can be endued with self-organized regular stripe patterns of nanoscopic lengthscale by Langmuir-Blodgett transfer of organic monolayers. Here we consider the effect of periodically prepatterned substrates on this process of pattern formation. It leads to a time periodic forcing of the oscillatory behavior at the meniscus. Utilizing higher order synchronization with this forcing, complex periodic patterns of predefined wavelength can be created. The dependence of the synchronization on the amplitude and the wavelength of the wetting contrast is investigated in one and two spatial dimensions and the resulting patterns are discussed. Furthermore, the effect of prepatterned substrates on the pattern selection process is investigated

Vortex-antivortex proliferation from an obstacle in thin film ferromagnets

Magnetization dynamics in thin film ferromagnets can be studied using a dispersive hydrodynamic formulation. The equations describing the magnetodynamics map to a compressible fluid with broken Galilean invariance parametrized by the longitudinal spin density and a magnetic analog of the fluid velocity that define spin-density waves. A direct consequence of these equations is the determination of a magnetic Mach number. Micromagnetic simulations reveal nucleation of nonlinear structures from an impenetrable object realized by an applied magnetic field spot or a defect. In this work, micromagnetic simulations demonstrate vortex-antivortex pair nucleation from an obstacle. Their interaction establishes either ordered or irregular vortex-antivortex complexes. Furthermore, when the magnetic Mach number exceeds unity (supersonic flow), a Mach cone and periodic wavefronts are observed, which can be well-described by solutions of the steady, linearized equations. These results are reminiscent of theoretical and experimental observations in Bose-Einstein condensates, and further supports the analogy between the magnetodynamics of a thin film ferromagnet and compressible fluids. The nucleation of nonlinear structures and vortex-antivortex complexes using this approach enables the study of their interactions and effects on the stability of spin-density waves.Comment: 23 pages, 7 figure

Gradient bounds for a thin film epitaxy equation

We consider a gradient flow modeling the epitaxial growth of thin films with slope selection. The surface height profile satisfies a nonlinear diffusion equation with biharmonic dissipation. We establish optimal local and global wellposedness for initial data with critical regularity. To understand the mechanism of slope selection and the dependence on the dissipation coefficient, we exhibit several lower and upper bounds for the gradient of the solution in physical dimensions $d\le 3$

Continuation for thin film hydrodynamics and related scalar problems

This chapter illustrates how to apply continuation techniques in the analysis of a particular class of nonlinear kinetic equations that describe the time evolution through transport equations for a single scalar field like a densities or interface profiles of various types. We first systematically introduce these equations as gradient dynamics combining mass-conserving and nonmass-conserving fluxes followed by a discussion of nonvariational amendmends and a brief introduction to their analysis by numerical continuation. The approach is first applied to a number of common examples of variational equations, namely, Allen-Cahn- and Cahn-Hilliard-type equations including certain thin-film equations for partially wetting liquids on homogeneous and heterogeneous substrates as well as Swift-Hohenberg and Phase-Field-Crystal equations. Second we consider nonvariational examples as the Kuramoto-Sivashinsky equation, convective Allen-Cahn and Cahn-Hilliard equations and thin-film equations describing stationary sliding drops and a transversal front instability in a dip-coating. Through the different examples we illustrate how to employ the numerical tools provided by the packages auto07p and pde2path to determine steady, stationary and time-periodic solutions in one and two dimensions and the resulting bifurcation diagrams. The incorporation of boundary conditions and integral side conditions is also discussed as well as problem-specific implementation issues

Tuning the Polar States of Ferroelectric Films via Surface Charges and Flexoelectricity

Using the self-consistent Landau-Ginzburg-Devonshire approach we simulate and analyze the spontaneous formation of the domain structure in thin ferroelectric films covered with the surface screening charge of the specific nature (Bardeen-type surface states). Hence we consider the competition between the screening and the domain formation as alternative ways to reduce the electrostatic energy and reveal unusual peculiarities of distributions of polarization, electric and elastic fields conditioned by the surface screening length and the flexocoupling strength. We have established that the critical thickness of the film and its transition temperature to a paraelectric phase strongly depend on the Bardeen screening length, while the flexocoupling affects the polarization rotation and closure domain structure and induces ribbon-like nano-scale domains in the film depth far from the top open surface. Hence the joint action of the surface screening (originating from e.g. the adsorption of ambient ions or surface states) and flexocoupling may remarkably modify polar and electromechanical properties of thin ferroelectric films.Comment: 33 pages, 5 figure