Preprint
Article

Scattering over Varying Amplification Grating

Altmetrics

Downloads

83

Views

14

Comments

0

A peer-reviewed article of this preprint also exists.

This version is not peer-reviewed

Submitted:

22 January 2024

Posted:

24 January 2024

You are already at the latest version

Alerts
Abstract
The scattering pattern from a narrow absorbing/amplifying grating is investigated. A simple model of a narrow amplifying grating is solved exactly numerically and approximately analytically for the regime where the beam’s wavelength is much shorter than the grating’s wavelength. The main results are: The incident angle divides the scattering pattern into two regimes: below and above the incident angles. The former regime has a weak dependence on the incident angle but has a strong dependence on the scattering one. In this regime, a new grating formula is derived. The opposite occurs in the latter regime, which is very sensitive to the incident angle, but has only weak dependence on the scattering angle. Consequently, at certain incident angles, all the scattering is concentrated in the first regime, i.e., all scattering angles are lower than the incident angle.
Keywords: 
Subject: Physical Sciences  -   Optics and Photonics

1. Introduction

The periodicity of crystals and gratings exhibits universal scattering patterns, which are manifested by the well-known grating formula and grating’s spectral resolution [1].
In general, the grating’s scattering is much more complicated and it depends on the grating’s material, profile, and grooving method [1,2,3]. In particular, plasmon excitations in the metal grating are responsible for anomalous scattering [4,5,6]. Moreover, it was found that even dielectric gratings exhibit anomalous scattering behavior [7,8]. These results were in clear contrast to the traditional belief that narrow dielectric gratings can have only a small effect on the incident light beam [9]. It was shown that even when the dielectric grating has weak modulation, the reflection coefficient can be very high [10,11]. The modulation can be either in the refraction index [12] or in the grating’s grooves [7]. These results show that in principle the reflectance at specific incident angles can be arbitrarily close to 100%. It should be stressed that these angles (with 100% reflectance) have nothing to do with the grating formula.
This anomalous conduct is related to a specific case of Fano resonances [8,13,14,15]. Recently, it was shown that in theory perfect reflection (zero transmission) can occur for arbitrarily narrow dielectric grating layer [16]. The layer thickness can be arbitrarily narrow and the grating modulation amplitude can be arbitrarily low and still the effect of 100% reflectance occurs.
This conduct resembles a similar effect in Quantum Mechanics, where the transmission via an oscillating barrier can be suppressed completely at certain incoming particles’ energies. Even though the quantum phenomenon was discovered after its optical equivalent, it was thoroughly investigated for both scattering scenarios [17,18,19,20,21,22,23] as well as for transition between adjacent quantum wells [24,25,26].
This effect occurs when the particle’s energy approximately equals the oscillations’ quanta. However, this is a relatively simple quantum system. When the system becomes more complex the dynamics are much less intuitive. In particular, it has been shown that if the oscillating potential is a well, which is located in an opaque barrier, then this system exhibits complex behavior. The current is extremely sensitive to the system’s parameter, it can either increase substantially (activation and elevation) or almost vanish (complete suppression) [27,28,29,30].
When the barrier is absent this complex behavior vanishes. The reason is that the system’s behavior is determined by the ratio between the oscillating potential’s amplitude and the wavenumber in its vicinity. If this ratio is imaginary no complex behavior appears.
The barrier’s presence makes sure that the wavenumber is approximately real, and since the potential amplitude is real then the ratio is approximately real as well.
However, if the barrier is absent then the ratio can still be real provided the potential amplitude is imaginary. In that case, some of this complex behavior reappears.
While the physical application of an imaginary barrier is not entirely clear, the optical counterpart, i.e., imaginary wavenumber and imaginary refractive index are very common.
These conclusions can be applied to optical equivalent systems where the wavefunction is replaced with the electromagnetic field, and the oscillating potential is replaced with a spatial dielectric (grating as in Ref. [16] but for complex dielectric coefficient).
Then, similarly, the complex structure’s dynamics can reappear by either putting the grating in a low index of refraction environment or by replacing the dielectric grating with an amplifying/absorbing one. In what follows we will focus on the latter: the scattering pattern of an amplifying/absorbing grating is going to be investigated.
Modulated amplification can be done by changing gradually the pumping over the grating causing a gradual change in the population inversion.
However, it is much simpler to create a uniform amplifying layer and attach to it a modulated absorber. If both layers are considerably narrower than the wavelength, then the two-layer system can be regarded as a single layer with modulated amplification/absorption.

2. The system under study

The system under study is a plane wave’s scattering from a grating. This scenario is presented in Figure 1. A TE-polarized plane-wave beam hits an active grating at an incident angle θ . The grating is located on the x=0 plane. The gratings variations occur in the z-direction, in which case the electromagnetic Maxwell’s wave equation can be written
2 x 2 E + 2 z 2 E + k 2 x , z E = 0
where the electric field is polarized in the y-direction (TE polarization) E = y ^ E x , z , and k 2 x , z satisfies
k 2 x , z = k 0 2 1 + i δ x d ε " + δ ε " cos K z
where the grating’s frequency is the reciprocal of the groove spacing Λ (the pitch), i.e.
K 2 π / Λ .
In (2) we followed Ref. [16] and chose a delta function to represent the narrow grating’s layer. As was explained in [16] (and references therein) a delta function is an excellent approximation for layers, which are shorter than the electromagnetic field’s wavelength. To be more specific, a dielectric with the following properties
ε x , y = ε 0 x , y 0 < x < d 0 e l s e where d is the narrow layer’s thickness, can be approximated with great accuracy by δ x 0 d ε 0 x d x . Since in what follows we focus on the limit of extremely narrow layers, this is an excellent approximation. Moreover, the validity of this approximation is improved when the incident angle increases, since the x-component of the beam’s wavenumber decreases. We will see that it is easier to measure the effect in the large incident angle domain. When the incident angle approaches π/2 the approximation is accurate.

3. Exact Numerical Solution

Let θ be the incident angle, therefore the wavenumber is k = k 0 x ^ cos θ + z ^ sin θ (see Figure 1), and the electromagnetic field can be written in all space as
E x , z / E 0 = exp i k 0 x cos θ + z sin θ + n r n exp i q n x + i p n z x < 0 n t n exp i q n x + i p n z x > 0
where p n k 0 sin θ + K n and q n k 0 2 p n 2 = k 0 1 sin θ + λ n / Λ 2 , i.e., q 0 = k 0 cos θ and q 0 = k 0 cos θ . The boundary criterions on the grating requires
E x = + 0 , z = E x = 0 , z and   x E x , z x = + 0 x E x , z x = 0 = i d k 0 2 ε " + δ ε " cos K z E x = 0 , z
and using the dimensionless parameters
s λ / Λ ,   P n p n / k 0 = sin θ + n s Q n q n / k 0 = 1 sin θ + n s 2 ,   α 1 2 d k 0 ε " ,   and   β = d k 0 1 4 δ ε "
the problem reduces to the following difference equation
Q n + α t n + β t n + 1 + t n 1 = Q 0 δ n
which can be rewritten in a more compact matrix form
M 0 t = v
where
M 0       Q 1 + α β       β Q 0 + α β       β Q 1 + α       t t 1 t 0 t 1 ,   and   v 0 Q 0 0 .
The numerical solution of this set of equations is presented in Figure 2.
Eq.(7) illustrates the fact, which was mentioned in the introduction, that the dynamics is governed by the ratio between oscillations amplitude β and Q n . When this ratio is imaginary (note that β is the imaginary part of the oscillations’ amplitude), the spectrum cannot have a complex structure.
As can be seen from this figure, there are incident angles in which the transmission coefficient of the non-negative, i.e. n 0 , modes vanish. A small change in the incoming angle (~1.4% in Figure 2) reduces the positive modes’ transmission probability from maximum value to zero. On the other hand, the incoming angle has a small effect on the negative modes (see Figure 2). Unlike the positive modes regime, where the differences between mode’s amplitudes is small, on the negative regime, the spectrum consists of ridges and valleys. This spectrum structure and behavior is similar to the quantum system of Refs. [29].
In Figure 3 the grating’s scattering pattern is presented for these two incident angles. When the incident angle is sin θ = 0.729 then the e-m field is scattered only in the sin θ s c < 0.729 angles. Nothing is scattered at angles, which exceed the incident angle sin θ s c > 0.729 . Below the incident angle, the e-m field is scattered at only specific angles.
On the other hand when the incident angle is sin θ = 0.737 the e-m is scattered at all angles even beyond the incident angle sin θ s c > 0.737 . Below the incident angle the scattering pattern is similar – in this regime, the e-m field is scattered at only specific angles.
Since all the n 0 modes vanish simultaneously, one can focus on the central mode n = 0 . In Figure 4 the amplitude of the central mode is presented as a function of the incident angle. As can be seen, beyond a certain angle, the transmission coefficient vanishes for several angles.
When the modulation depth β varies, so does the number of incident angles, in which the non-negative modes vanish. As can be seen in Figure 5, the larger the modulation depth, the larger the number of incident angles, in which the non-negative modes vanish.

4. Approximate Analytical Derivation

The numerical studies shows that this effect occurs when the grating constant is considerably larger than the incident beam’s wavelength. In this regime, Eq.(7) can be simplified to allow for analytical derivation. In this regime,
Q n q n / k 0 cos θ n s tan θ
which follows that Eq.(7) can be simplified to
cos θ + α β n s β tan θ t n + t n + 1 + t n 1 = Q 0 β δ n
The solution of this difference equation can be solved with Bessel functions. Using Eq.(9.1.27) of Ref. [31]
F ν 1 x + F ν + 1 x = 2 ν x F ν x the solution of (11) can be written
t n = A J n + n 0 u n < 0 C J n n 0 u n > 0
where n 0 cos θ + α s tan θ , u 2 β s tan θ , J n u are the Bessel functions of the first kind [31] and the coefficients A and C can be calculated by matching the solution and substituting the solution in (12), i.e.,
A J n 0 u = C J n 0 u ,   and   s β tan θ n 0 C J n 0 u + C J 1 n 0 u + A J 1 + n 0 u = Q 0 β .
The solution therefore reads
A = Q 0 β J n 0 u s β tan θ n 0 J n 0 u J n 0 u + J n 0 u J 1 n 0 u + J n 0 u J 1 + n 0 u
C = Q 0 β J n 0 u s β tan θ n 0 J n 0 u J n 0 u + J n 0 u J 1 n 0 u + J n 0 u J 1 + n 0 u
Now we have an approximate solution of the entire spatial spectrum, from which two main results can be derived: (A) the incident angles in which the scattered e-m is bounded to the lower angles, i.e., the scattered angles are lower than the incident angles. (B) The scattering angles below the incident one.

5. Calculating the Incident Angles for Which Scattering Is Bounded

From (12) and (14) if follows that suppression of high angles scattering occurs for C = 0 , i.e. for
J n 0 ( u ) = 0
and in which case all the n 0 modes vanish as well.
Since this approximation was derived for the limit of s < < 1 then both relations 2 β s tan θ > > 1 and n 0 > > 1 hold. Therefore, using Eq.(9.3.3) of Ref. [31],
J ν ν sec B 2 π ν tan B cos ν tan B B π / 4
where in this case
B arccos cos θ + α 2 β   and   ν n 0 .
For a solution, in which the argument of the cosine function of (16) is small whilst still s < < 1 , the parameter B must be accordingly small. Therefore, we seek a solution around
cos θ 2 β α ,
n 0 B 3 / 3 = m + 3 4 π
in which case, (15) and (16) yields
Solving (19) for β , we finally obtain the condition for high angles suppression
β cos θ + α 2 1 + 1 2 3 s tan θ cos θ + α m + 3 4 π 2 / 3
or equivalently, one can find the incident angles θ m for which high angle scattering does not occur
cos θ m 2 β 1 + 1 2 3 s 2 β 2 β α 2 1 m + 3 4 π 2 / 3 1 α
In Figure 6 the approximation (20) is plotted on top of the exact numerical simulation (Figure 5). As can be seen, this approximation shows high agreement with numerical results.

6. Scattering Angles

The scattered beam takes the following form
E x , z / E 0 = n t n exp i q n x + i p n z = n A J n + n 0 2 β s tan θ exp i k 0 1 sin θ + n s 2 x + sin θ + n s z
where A is the constant (13).
To calculate the scattering angles, one can use the fact that s < < 1 . In this case the summation can be approximated by an integral
E x , z / E 0 s 1 A d η J η 0 η / s 2 β s tan θ exp i k 0 1 sin θ + η 2 x + sin θ + η z
Using the method of stationary phase approximation, this integral can be evaluated at
sin θ s c sin θ = η
Note that (24) is the classical grating formula.
Therefore, the scattering amplitude at the scattering angle θ s c is approximately
A J η 0 sin θ s c sin θ / s 2 β s tan θ 2
Again, since s < < 1 one can use Eq.(9.3.3) of Ref. [31], in which case
B 2 1 cos θ + α n s tan θ 2 β
the angles which yields the maximum value of (25) must satisfy
n 0 B 3 / 3 π / 4 = m π
Substituting (26) in (27) and solving the angle, while keeping in mind that sin θ o u t sin θ = η one finally get the grating formula
sin θ o u t m sin θ β 3 π n 0 m + 1 4 2 / 3 2 + cos θ + α cot θ
As can be seen in Figure 7 the scattered angles θ o u t m for m = 0 , 1 , and 2 shows high agreement with the approximation (28). It should be stressed that the grating formula (28) appears on top of the classical one (24).

7. Physical Realization

The realization of this grating requires simultaneous changes in both the real and imaginary components of the refractive index, however, they do not change in the same manner.
To simplify the model, the harmonic change occurs in the imaginary component of the dielectric constant δ ε " cos K z . Therefore, the real and imaginary parts of the refraction index are
n = 2 1 / 2 δ ε " cos K z / 1 + 1 + δ ε " cos K z 2
n = 2 1 / 2 1 + 1 + δ ε " cos K z 2
These components are plotted in Figure 8 for δ ε " = 5 . As can be seen, for this value, the real part of the refractive index changes between 1 and 1.7, while the imaginary part varies between -1.5 to +1.5. Then, if the layer thickness is d / λ 0.03 then β = d k 0 1 4 δ ε " 0.24 , which is sufficient to measure several angles which shows this effect (according to Figure 5). It should be noted that if we seek the effect at high, i.e., grazing, angles, the layers width can be increased substantially, while keeping the narrow layer approximation still. In this case multiple angles should be detected. Moreover, Most of this work concentrated, for simplicity, on the α = 0 scenario, however, by varying α the effect will appear even for very weak β , in which case, it will be easier to measure the effect. Nevertheless, a full study of this effect for different values of α and different modulations is beyond the scope of this paper.

8. Summary

The scattering from an active grating layer was analyzed. In this model, the imaginary component of the dielectric coefficient varies harmonically. The grating’s layer width is shorter than the optical beam wavelength, while the grating’s wavelength is longer than the beam’s wavelength. A new grating formula is formulated, which is the profile on top of the classical grating formula. An approximate analytical expression for this formula was derived. Besides the scattering orders, a special scattering pattern appears for certain incident angles.
It was found that in some cases, the scattering angles θ o u t cannot exceed the incident angle θ , i.e., for specific θ = θ m the scattered angles obey θ o u t θ .

References

  1. Wilcox, C.H. Scattering Theory for Diffraction Gratings; Springer-Verlag: New York, NY, USA, 1984. [Google Scholar]
  2. Gratings: Theory and Numeric Applications; Popov, E. (Ed.) Presses universitaires de Provence; PUP; 2012. [Google Scholar]
  3. Palmer, C.; Loewen, E. Diffraction Grating Handbook; Newport Corporation, 2005. [Google Scholar]
  4. Piper, J.R.; Fan, S. Total absorption in a Graphene Monolayer in the Optical Regime by Critical Coupling with a Photonic Crystal Guided Resonance. ACS Photonics 2014, 1, 347–353. [Google Scholar] [CrossRef]
  5. Hutley, M.C.; Maystre, D. Total absorption of light by a diffraction grating. Opt. Commun. 1976, 19, 431–436. [Google Scholar] [CrossRef]
  6. Jory, M.J.; Vukusic, P.S.; Sambles, J.R. Development of a prototype gas sensor using surface plasmon resonance on gratings. Sens. Actuators B 1994, 17, 203–209. [Google Scholar] [CrossRef]
  7. Mashev, L.; Popov, E. Zero order anomaly of dielectric coated gratings. Optics. Commun. 1985, 55, 377. [Google Scholar] [CrossRef]
  8. Popov, E.; Mashev, L.; Maystre, D. Theoretical study of the anomalies of coated dielectric gratings. Optica Acta 1986, 33, 607–619. [Google Scholar] [CrossRef]
  9. Burckhardt, C.B. Efficiency of a Dielectric Grating. J. Opt. Soc. Am. 1967, 57, 601. [Google Scholar] [CrossRef]
  10. Tibuleac, S.; Magnusson, R. Reflection and Transmission guided-mode resonance filters. J. Opt. Soc. Am. A 1997, 14, 1617. [Google Scholar] [CrossRef]
  11. Tibuleac, S.; Magnusson, R. Narrow-linewidth bandpass filters with diffractive thin-film layers. Opt. Lett. 2001, 26, 584. [Google Scholar] [CrossRef] [PubMed]
  12. Wang, S.S.; Magnusson, R.; Bagby, J.S.; Moharam, M.G. Guided-mode resonance in planar dielectric-layer diffraction gratings. J. Opt. Soc. Am. A 1990, 7, 1470. [Google Scholar] [CrossRef]
  13. Fano, U. The theory of anomalous diffraction gratings and of quasi-stationary waves on metallic surfaces (Sommerfield’s waves). J. Opt. Soc. Am. 1941, 31, 213–222. [Google Scholar] [CrossRef]
  14. Hessel, A.; Oliner, A.A. A new theory of Wood’s anomalies on optical gratings. Appl. Opt. 1965, 4, 1275–1297. [Google Scholar] [CrossRef]
  15. Fan, S.; Suh, W.; Joannopoulos, J.D. Temporal coupled-mode theory for Fano resonant mode in optical resonators. J. Opt. Soc. Am. A 2003, 20, 569–572. [Google Scholar] [CrossRef] [PubMed]
  16. Granot, E. Derivation of analytical expressions for anomalous reflection in the limit of zero thickness and weakly modulated dielectric grating. J. Opt. Soc. Am. A 2022, 39, 2205–2213. [Google Scholar] [CrossRef] [PubMed]
  17. Cota, E.; Jose, J.V.; Rojas, F. Electronic transmission through time-periodic oscillating barriers. Nanostructured Mater. 1993, 3, 349. [Google Scholar] [CrossRef]
  18. Bagwell, P.F.; Lake, R.K. Resonances in transmission through an oscillating barrier. Phys. Rev. B 1992, 46, 15329. [Google Scholar] [CrossRef] [PubMed]
  19. Li, W.; Reichl, L.E. Floquet scattering through a time-periodic potential. Phys. Rev. B 1999, 60, 15732. [Google Scholar] [CrossRef]
  20. Saraga, D.S.; Sassoli-de-Bianchi, M. On the One-Dimensional Scattering by Time-Periodic Potentials: General Theory and Application to Two Specific Models. Helv. Phys. Acta 1997, 70, 751. [Google Scholar]
  21. Martinez, D.F.; Reichl, L.E. Transmission Properties of the oscillating delta-function potential. Phys. Rev. B 2001, 64. [Google Scholar]
  22. Haavig, D.L.; Reifenberger, R. Dynamic transmission and reflection phenomena for a time-dependent rectangular potential. Phys. Rev. B 1982, 26, 6408. [Google Scholar] [CrossRef]
  23. Cai, W.; Hu, P.; Zheng, T.F.; Yudanin, B. Resonance of the one-dimensional electron transmission above a quantum well with dissipation. Phys. Rev. B 1990, 41, 3513. [Google Scholar] [CrossRef]
  24. Kayanuma, Y. Role of phase coherence in the transition dynamics of a periodically driven two-level system. Phys. Rev. A 1994, 50, 843. [Google Scholar] [CrossRef] [PubMed]
  25. Kayanuma, Y. Stokes phase and geometrical phase in a driven two-level system. Phys. Rev. A 1997, 55, 2495. [Google Scholar] [CrossRef]
  26. del Campo, A.; Garcia-Calderon, G.; Muga, J.G. Quantum transients. Phys.Rep. 2009, 476, 1. [Google Scholar] [CrossRef]
  27. Azbel, M.Y. Eigenstate assisted activation. Phys. Rev. Lett 1992, 68, 98. [Google Scholar] [CrossRef] [PubMed]
  28. Zangwill, G.; Granot, E. Eigenstate suppressed activation. Phys. B: Condens. Matter 2015, 461, 140. [Google Scholar] [CrossRef]
  29. Granot, E. Selected elevation in quantum tunneling. Europhys. Lett. 2003, 61, 817. [Google Scholar] [CrossRef]
  30. Granot, E. The Tunneling Current through Oscillating Resonance and the Sisyphus Effect. Adv. Cond. Matter Phys. 2017, 2017, 2435857. [Google Scholar]
  31. Abramowitz, M.; Stegun, A. Handbook of Mathematical Functions; Dover Publications: New York, NY, USA, 1965. [Google Scholar]
Figure 1. iystem schematic: diffraction by a thin grating. The arrows indicate the beam’s scattering directions. The incident beam is polarized in the y-direction, and the wavenumber is k = k 0 x ^ cos θ + z ^ sin θ .
Figure 1. iystem schematic: diffraction by a thin grating. The arrows indicate the beam’s scattering directions. The incident beam is polarized in the y-direction, and the wavenumber is k = k 0 x ^ cos θ + z ^ sin θ .
Preprints 97040 g001
Figure 2. The transmission coefficient amplitude of the spectral modes for three different incoming angles. The parameters used in this simulation were: β = 0.414 , α = 0   s = 0.003 and 521 modes were used.
Figure 2. The transmission coefficient amplitude of the spectral modes for three different incoming angles. The parameters used in this simulation were: β = 0.414 , α = 0   s = 0.003 and 521 modes were used.
Preprints 97040 g002
Figure 3. 2D image of the scattering pattern. The dashed line corresponds to direction of the incident beam. The darker the color the higher is the light intensity E x , z 2 . The parameters in this case were: s = 0.003 , β = 0.414 , α = 0 and 521 modes were used.
Figure 3. 2D image of the scattering pattern. The dashed line corresponds to direction of the incident beam. The darker the color the higher is the light intensity E x , z 2 . The parameters in this case were: s = 0.003 , β = 0.414 , α = 0 and 521 modes were used.
Preprints 97040 g003
Figure 4. The transmission coefficient of the central, i.e., n=0, mode as a function of the incident angle. In this simulation β = 0.45 , α = 0 , s = 0.01 and 501 modes were used.
Figure 4. The transmission coefficient of the central, i.e., n=0, mode as a function of the incident angle. In this simulation β = 0.45 , α = 0 , s = 0.01 and 501 modes were used.
Preprints 97040 g004
Figure 5. The transmission coefficient amplitude of the central mode as a function of the modulation depth and the incident angle. The darker the color the higher is the square of the coefficient amplitude. The same parameters of Figure 4. The parameters in this simulation were s = 0.01 , and α = 0 and 521 modes were used.
Figure 5. The transmission coefficient amplitude of the central mode as a function of the modulation depth and the incident angle. The darker the color the higher is the square of the coefficient amplitude. The same parameters of Figure 4. The parameters in this simulation were s = 0.01 , and α = 0 and 521 modes were used.
Preprints 97040 g005
Figure 6. Zoom in of Figure 5. The dashed curves represents the approximation (20) for four integers m.
Figure 6. Zoom in of Figure 5. The dashed curves represents the approximation (20) for four integers m.
Preprints 97040 g006
Figure 7. Zoom in of Figure 3 with the same parameters. The solid line represents the incident beam’s direction, and the dashed lines represents the scattered directions according to the approximation of the grating’s formula(28).
Figure 7. Zoom in of Figure 3 with the same parameters. The solid line represents the incident beam’s direction, and the dashed lines represents the scattered directions according to the approximation of the grating’s formula(28).
Preprints 97040 g007
Figure 8. The real (dashed curve) and imaginary (solid curve) parts of the refraction index along the grating (z-axis) according to (30) and (29) respectively.
Figure 8. The real (dashed curve) and imaginary (solid curve) parts of the refraction index along the grating (z-axis) according to (30) and (29) respectively.
Preprints 97040 g008
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2024 MDPI (Basel, Switzerland) unless otherwise stated