Skip to content

Optimizing Floquet engineering for non-equilibrium steady states with gradient-based methods

MetadataDetails
Publication Date2023-07-26
JournalSciPost Physics
AuthorsAlberto Castro, Shunsuke Sato
InstitutionsUniversidad de Zaragoza, University of Tsukuba
Citations3
AnalysisFull AI Review Included

Technical Documentation & Analysis: Floquet Engineering in Diamond NV Centers

Section titled “Technical Documentation & Analysis: Floquet Engineering in Diamond NV Centers”

Reference: Castro, A. and Sato, S. A. (2023). Optimizing Floquet engineering for non-equilibrium steady states with gradient-based methods. SciPost Phys. 15, 029.


This research demonstrates a powerful methodology for controlling the properties of quantum materials in non-equilibrium steady states (NESS) using optimized periodic driving fields. The findings are highly relevant for the development of advanced diamond-based quantum technologies, particularly those utilizing the Nitrogen-Vacancy (NV) center.

  • Core Achievement: Successful application of Quantum Optimal Control Theory (QOCT) combined with a Lindblad master equation to optimally tune NESS properties in an open quantum system.
  • Model System: The Nitrogen-Vacancy (NV) center in diamond, driven by time-periodic magnetic fields (Floquet engineering).
  • Optimization Goal: Maximization and minimization of the time-averaged spin component, <<Sz>>, by optimizing the Fourier coefficients of the driving fields.
  • Exotic State Preparation: The technique achieved “exotic” NESSs, demonstrating spin expectation values (e.g., <<Sz>> ≈ 0.38) significantly exceeding the maximum allowed in thermal equilibrium (≈ 0.14).
  • Experimental Relevance: By utilizing the Lindblad equation, the framework naturally incorporates decoherence and dissipation (rate Îł), ensuring the optimized fields are relevant for experimentally realizable quantum control.
  • Material Exploration: This work establishes periodic perturbations as a novel degree of freedom for engineering material functionalities, extending material exploration beyond equilibrium phases.

The following hard data points were extracted from the theoretical model and optimization results concerning the NV center system.

ParameterValueUnitContext
Model Fixed UnitNz = 1DimensionlessNV center Hamiltonian parameter (sets energy scale)
Transverse CouplingNxy = 0.05DimensionlessNV center Hamiltonian parameter
Static Magnetic Field (z-axis)Bz = 0.3DimensionlessNV center Hamiltonian parameter
Driving Field Scaling FactorB0 = 0.1DimensionlessAmplitude scaling for periodic perturbation V(u, t)
Inverse Temperature (Fixed)ÎČ = 31/kBTUsed for NESS calculations
Dissipation Rate ConstantÎł = 0.2DimensionlessStandard rate constant used in the Lindblad equation
Thermal Equilibrium Max <Sz>≈ 0.14DimensionlessMaximum Sz achievable by temperature variation (ÎČ > 0)
Optimized Max <<Sz>>≈ 0.38DimensionlessTime-averaged spin component achieved via QOCT (Îș=4.0)
Amplitude Constraint (Max Sz)Îș = 4.0DimensionlessMaximum bound on Fourier coefficients
Cutoff FrequencyωM = 2.0 NzDimensionlessMaximum frequency component used in the driving field

The optimization of Non-Equilibrium Steady States (NESS) was achieved by integrating Quantum Optimal Control Theory (QOCT) with open quantum system dynamics.

  1. Open Quantum System Modeling: The system dynamics were governed by a Lindblad-type master equation (Eq. 1) incorporating a time-periodic Hamiltonian H(t) and Lindblad operators Vij to account for dissipation and decoherence effects (Markovian approximation).
  2. Periodic Driving Parametrization: The external periodic magnetic fields, gx(t) and gy(t), were defined using truncated Fourier expansions (Eq. 24). The control parameters $u$ were the coefficients of these expansions, allowing for “multicolor periodic driving.”
  3. Optimization Target Definition: The objective function G(u) was defined as the time-average of the expectation value of an observable A (specifically, the Sz spin component) over one Floquet period T.
  4. NESS Calculation in Floquet-Liouville Space: The periodic NESS solution ρu(t) was found by solving the linear homogeneous equation in the extended Floquet-Liouville space, which accounts for the time-dependence of the Lindbladian.
  5. Gradient Derivation: A novel procedure was developed to compute the gradient (∂G/∂uk) of the objective function with respect to the control parameters, circumventing difficulties associated with standard adjoint methods for periodic equations.
  6. Constrained Optimization: The Sequential Least-Squares Quadratic Programming (SLSQP) algorithm was used to perform the gradient-based optimization, subject to constraints on the amplitude of the driving field Fourier coefficients (|uj| ≀ Îș), mimicking experimental limitations.

This research highlights the critical role of high-quality diamond substrates and precision engineering in achieving advanced quantum control. 6CCVD is uniquely positioned to supply the necessary materials and customization services required to replicate and extend this Floquet engineering research into practical devices.

Category6CCVD Solution & Value Proposition
Applicable MaterialsOptical Grade Single Crystal Diamond (SCD): The NV center is the core component. Achieving the low dissipation rates (Îł) required for optimal control necessitates ultra-high purity SCD with long coherence times (T2). 6CCVD supplies high-quality SCD plates, essential for minimizing background defects and maximizing quantum performance.
Custom Dimensions & GeometryPrecision Plates and Wafers: We offer custom SCD plates (up to 10x10mm) and large-area PCD wafers (up to 125mm). This capability supports the integration of complex microwave delivery structures necessary to generate the optimized periodic magnetic fields gx(t) and gy(t).
Surface EngineeringUltra-Smooth Polishing (Ra < 1nm): Surface quality is paramount for minimizing decoherence. 6CCVD guarantees Ra < 1nm polishing on SCD, ensuring minimal surface scattering losses and optimal coupling for microwave control fields.
Metalization ServicesIntegrated Quantum Control Circuits: The periodic driving fields are typically delivered via patterned microwave transmission lines. 6CCVD provides in-house custom metalization services, including Ti, Pt, Au, W, Pd, and Cu, allowing researchers to define and pattern the precise geometries required for QOCT field delivery directly onto the diamond substrate.
Thickness ControlFlexible Substrate Thickness: We offer SCD thicknesses ranging from 0.1”m to 500”m. This flexibility is crucial for researchers optimizing the placement of NV centers relative to surface control structures.
Engineering SupportMaterial Selection for Low Dissipation: The success of the QOCT scheme depends heavily on the dissipation rate (γ). 6CCVD’s in-house PhD team provides expert consultation to select diamond materials (e.g., isotopic purity, surface termination) that align with the low-decoherence assumptions of advanced optimal control models for similar NV Center Quantum Control projects.
Global LogisticsReliable Global Shipping: We ensure prompt and secure global delivery (DDU default, DDP available) of high-value diamond materials, supporting international research collaborations.

For custom specifications or material consultation, visit 6ccvd.com or contact our engineering team directly.

View Original Abstract

Non-equilibrium steady states are created when a periodically driven quantum system is also incoherently interacting with an environment - as it is the case in most realistic situations. The notion of Floquet engineering refers to the manipulation of the properties of systems under periodic perturbations. Although it more frequently refers to the coherent states of isolated systems (or to the transient phase for states that are weakly coupled to the environment), it may sometimes be of more interest to consider the final steady states that are reached after decoherence and dissipation take place. In this work, we demonstrate how those final states can be optimally tuned with respect to a given predefined metric, such as for example the maximization of the temporal average value of some observable, by using multicolor periodic perturbations. We show a computational framework that can be used for that purpose, and exemplify the concept using a simple model for the nitrogen-vacancy center in diamond: the goal in this case is to find the driving periodic magnetic field that maximizes a time-averaged spin component. We show that, for example, this technique permits to prepare states whose spin values are forbidden in thermal equilibrium at any temperature.