Tilnærming til diskret dipol - Discrete dipole approximation
Diskret dipol -tilnærming ( DDA ), også kjent som koblet dipol -tilnærming , er en metode for å beregne spredning av stråling av partikler med vilkårlig form og periodiske strukturer. Gitt et mål for vilkårlig geometri, søker man å beregne dens sprednings- og absorpsjonsegenskaper ved en tilnærming til kontinuummålet med et begrenset utvalg av små polariserbare dipoler . Denne teknikken brukes i et flertall av anvendelser inklusiv nanofotonikk , radar spredning, aerosol fysikk og astro .
Enkle konsepter
Grunnideen til DDA ble introdusert i 1964 av DeVoe som brukte den for å studere de optiske egenskapene til molekylære aggregater; retardasjonseffekter var ikke inkludert, så DeVoes behandling var begrenset til aggregater som var små sammenlignet med bølgelengden. DDA, inkludert retardasjonseffekter, ble foreslått i 1973 av Purcell og Pennypacker som brukte den til å studere interstellare støvkorn. Enkelt sagt, DDA er en tilnærming til kontinuummålet med et begrenset utvalg av polariserbare punkter. Punktene får dipolmomenter som svar på det lokale elektriske feltet. Dipolene interagerer med hverandre via sine elektriske felt, så DDA blir også noen ganger referert til som den koblede dipol -tilnærmingen.
Naturen gir den fysiske inspirasjonen til DDA - i 1909 viste Lorentz at de dielektriske egenskapene til et stoff kan være direkte relatert til polariserbarheten til de enkelte atomene det var sammensatt av, med et spesielt enkelt og eksakt forhold, Clausius -Mossotti -forholdet (eller Lorentz-Lorenz), når atomene er plassert på et kubikkgitter. Vi kan forvente at, akkurat som en kontinuumrepresentasjon av et fast stoff er passende på lengdeskalaer som er store sammenlignet med den interatomiske avstanden, kan en rekke polariserbare punkter nøyaktig tilnærme responsen til et kontinuumsmål på lengdeskalaer som er store sammenlignet med interdipol separasjon.
For et begrenset utvalg av punktdipoler kan spredningsproblemet løses nøyaktig, så den eneste tilnærmingen som er tilstede i DDA er erstatning av kontinuummålet med en rekke N-punktdipoler. Erstatningen krever spesifikasjon av både geometrien (plasseringen av dipolene) og dipolens polariserbarhet. For monokromatiske hendelsesbølger kan den selvkonsistente løsningen for de oscillerende dipolmomentene bli funnet; fra disse beregnes absorpsjons- og spredningstverrsnittene. Hvis DDA -løsninger oppnås for to uavhengige polarisasjoner av hendelsesbølgen, kan den komplette amplituden -spredningsmatrisen bestemmes. Alternativt kan DDA stammer fra volumintegralligning for det elektriske feltet . Dette fremhever at tilnærmingen av punktdipoler tilsvarer den å diskretisere integralligningen, og reduseres dermed med avtagende dipolstørrelse.
Med erkjennelsen av at polariserbarhetene kan være tensorer, kan DDA enkelt påføres anisotrope materialer. Utvidelsen av DDA til å behandle materialer med magnetisk følsomhet uten null er også grei, selv om magnetiske effekter i de fleste applikasjoner er ubetydelige.
Utvidelser
Fremgangsmåten ble forbedret ved Draine , Flatau, og Goodman som påføres hurtig Fourier-transformasjon for å beregne konvolusjon problem som oppstår i DDA som tillot å beregne spredningen av store mål. De distribuerte diskret dipol -tilnærming åpen kildekode DDSCAT. Det er nå flere DDA -implementeringer, utvidelser til periodiske mål og partikler plassert på eller i nærheten av et plant underlag. og sammenligninger med eksakt teknikk ble publisert. Andre aspekter som validitetskriteriene for den diskrete dipoltilnærmingen ble publisert. DDA ble også utvidet til å anvende rektangulære eller kuboide dipoler som er mer effektive for sterkt oblate eller prolate partikler.
Diskrete tilnærmingskoder for dipolen
Det er anmeldelser så vel som publisert sammenligning av eksisterende koder. De fleste kodene gjelder vilkårlig formede inhomogene ikke-magnetiske partikler og partikkelsystemer i fritt rom eller homogent dielektrisk vertsmedium. De beregnede mengdene inkluderer vanligvis Mueller-matriser , integrerte tverrsnitt (utryddelse, absorpsjon og spredning), interne felt og vinkeloppløste spredte felt (fasefunksjon).
Generelle åpen kildekode DDA-koder
Disse kodene bruker vanligvis vanlige rutenett (kubikk eller rektangulær kuboid), konjugert gradientmetode for å løse store systemer med lineære ligninger, og FFT-akselerasjon av matrise-vektorproduktene som bruker konvolusjonsteorem. Kompleksiteten i denne tilnærmingen er nesten lineær i antall dipoler for både tid og minne.
Navn | Forfattere | Referanser | Språk | Oppdatert | Funksjoner |
---|---|---|---|---|---|
DDSCAT | Draine og Flatau | Fortran | 2019 (v. 7.3.3) | Kan også håndtere periodiske partikler og effektivt beregne nær felt . Bruker OpenMP -akselerasjon. | |
DDscat.C ++ | Choliy | C ++ | 2017 (v. 7.3.1) | Versjon av DDSCAT oversatt til C ++ med noen ytterligere forbedringer. | |
LEGG TIL EN | Yurkin, Hoekstra og bidragsytere | C | 2020 (v. 1.4.0) | Implementerer rask og grundig vurdering av et plant underlag, og tillater rektangulære kuboide vokser for svært oblate eller prolaterende partikler. Kan også beregne utslippsforbedring (forfallshastighet) av punktutsendere. Nær-felt beregningen er ikke veldig effektivt. Bruker Message Passing Interface (MPI) parallellisering og kan kjøres på GPU ( OpenCL ). | |
OpenDDA | McDonald | C | 2009 (v. 0.4.1) | Bruker både OpenMP- og MPI -parallellisering. Fokuserer på beregningseffektivitet. | |
DDA-GPU | Kieß | C ++ | 2016 | Kjører på GPU (OpenCL). Algoritmer er delvis basert på ADDA. | |
VIE-FFT | Sha | C/C ++ | 2019 | Beregner også nær felt og materialabsorpsjon. Navnet annerledes, men algoritmene er veldig like de som ble brukt i den vanlige DDA. | |
VoxScatter | Samuel Groth, Polimeridis og White | Matlab | 2020 | Inneholder forkondisjonering av akselerasjon | |
IF-DDA | PC Chaumet, A. Sentenac og D. Sentenac | Fortran og grafisk brukergrensesnitt skrevet i C ++ med Qt | 2020 | Idiotvennlig diskret dipol -tilnærming. Koden er tilgjengelig på github. Kan brukes med openMP og HDF5. |
Spesialiserte DDA -koder
Denne listen inneholder koder som ikke er kvalifisert for forrige seksjon. Årsakene kan omfatte følgende: kildekoden er ikke tilgjengelig, FFT -akselerasjon er fraværende eller redusert, koden fokuserer på spesifikke applikasjoner som ikke tillater enkel beregning av standard spredningsmengder.
Navn | Forfattere | Referanser | Språk | Oppdatert | Funksjoner | |
---|---|---|---|---|---|---|
DDSURF, DDSUB, DDFILM | Schmehl, Nebeker og Zhang | Fortran | 2008 | Streng håndtering av semi-uendelig substrat og endelige filmer (med vilkårlig partikkelplassering), men bare 2D FFT- akselerasjon brukes. | ||
DDMM | Mackowski | Fortran | 2002 | Beregner T-matrise , som deretter kan brukes til å effektivt beregne orienterings-gjennomsnittlige spredningsegenskaper. | ||
CDA | McMahon | Matlab | 2006 | |||
DDA-SI | Loke | Matlab | 2014 (v. 0.2) | Streng håndtering av underlaget, men ingen FFT -akselerasjon brukes. | ||
PyDDA | Python | 2015 | Reimplementering av DDA-SI | |||
e -DDA | Vaschillo og Bigelow | Fortran | 2019 (v. 2.0) | Simulerer elektron-energitapspektroskopi og katodoluminescens. Bygget på DDSCAT 7.1. | ||
DDEELS | Geuquet, Guillaume og Henrard | Fortran | 2013 (v. 2.1) | Simulerer elektron-energitapspektroskopi og katodoluminescens. Håndterer underlaget gjennom tilnærming til bildet, men ingen FFT -akselerasjon brukes. | ||
T-DDA | Edalatpour | Fortran | 2015 | Simulerer nærfeltets strålende varmeoverføring. Beregningsflaskehalsen er direkte matriseinversjon (ingen FFT -akselerasjon brukes). Bruker OpenMP og MPI parallellisering. | ||
IF-DDAM | PC Chaumet, A. Sentenac og D. Sentenac | Fortran og grafisk brukergrensesnitt skrevet i C ++ med Qt | 2021 | Idiotvennlig diskret dipol -tilnærming. Koden er tilgjengelig på github. Kan brukes med openMP og HDF5. Streng håndtering av flerlag (støtte guidet bølge eller plasmon og vilkårlig partikkelplassering). |
Galleri med former
Se også
- Beregningsmessig elektromagnetikk
- Mie teori
- Endelig-forskjell tidsdomene metode
- Metode for øyeblikk (elektromagnetikk)