Virtual ultrasound machine operating in a GHz to MHz frequency range for particle-based biomedical simulations

Deze paper introduceert een nieuw deeltjesgebaseerd virtueel ultrasone apparaat dat, dankzij een innovatieve drukoplosser en stabilisatieschema, ultrasone golven in het MHz- tot GHz-bereik effectief simuleert om toepassingen zoals microbel-microdrukkingsmechanismen voor geneesmiddellevering te modelleren.

De kern

Stel je voor dat je een virtuele geluidsmachine hebt die je in je computer kunt programmeren. Deze machine kan geluidsgolven (ultrasone geluiden) nabootsen die door water en kleine deeltjes reizen, precies zoals dat gebeurt in het menselijk lichaam.

Dit is wat de onderzoekers in dit artikel hebben gedaan. Hier is de uitleg in simpele taal, met een paar handige vergelijkingen:

1. Het Probleem: De "Tijdsdilemma"

Ultrasone geluiden (zoals die gebruikt worden voor echo's of medicijnen) bewegen razendsnel. Ze trillen heel snel (miljoenen tot miljarden keer per seconde).

• Het dilemma: Als je probeert dit te simuleren op een computer, heb je een groot probleem. De trillingen gaan zo snel, maar de vloeistof (zoals water) is ook een beetje "stroperig" (viskeus).
• De analogie: Stel je voor dat je een slak (de viskeuze vloeistof) en een supersnelle raket (de geluidsgolf) in dezelfde race zet. Als je de computer probeert te laten rekenen met de snelheid van de raket, moet je de slak ook in duizenden kleine stukjes hakken om hem op tijd te laten bewegen. Dat kost je computer onmogelijk veel tijd en energie.
• De oude oplossing: Wetenschappers gebruikten vaak twee verschillende programma's: één voor de snelle golven en één voor de trage vloeistof, en probeerden ze aan elkaar te plakken. Dat was als het proberen om een Lego-blokje en een stukje klei aan elkaar te lijmen; het zat niet goed en was heel lastig.

2. De Oplossing: De "Slimme Deeltjes" (usSDPD)

De onderzoekers hebben een nieuw soort computerprogramma bedacht, genaamd usSDPD.

• Hoe het werkt: In plaats van de vloeistof te zien als een gladde soep, zien ze het als een verzameling van duizenden kleine balletjes (deeltjes).
• De truc: Ze hebben een slimme rekenmethode bedacht (een "impliciete oplosser") die de computer toestaat om grotere stappen te maken in de tijd.
• De analogie: Stel je voor dat je een film bekijkt. Normaal moet je elke frame bekijken (heel traag). Met hun nieuwe methode mag de computer "sprongen" maken in de film, maar dan wel zo slim dat hij de snelle geluidsgolven niet mist. Hierdoor gaat de simulatie 40 keer sneller dan voorheen.

3. Het Stabiliteitsprobleem: De "Klap"

Geluidsgolven bestaan uit druk (samendrukken) en zuiging (uitrekken).

• Het probleem: Als de computer probeert de "zuigende" kant van de golf na te bootsen (negatieve druk), barsten de oude methoden vaak. De vloeistof "breekt" in de simulatie alsof er een scheur in het water komt.
• De oplossing: De onderzoekers hebben twee veiligheidsmaatregelen toegevoegd:
  1. Ze hebben de deeltjes iets kleiner gemaakt en dichter op elkaar gezet.
  2. Ze hebben een "kunstmatige druk" toegevoegd die voorkomt dat de deeltjes uit elkaar worden getrokken.
• De analogie: Het is alsof je een elastiekje hebt dat te snel uitrekt en dan knapt. Ze hebben een extra veertje toegevoegd aan het elastiekje, zodat het niet knapt, zelfs niet als je er heel hard aan trekt. Hierdoor kan de machine nu veilig simuleren hoe geluidsgolven door water gaan, zonder dat de simulatie crasht.

4. Wat kunnen ze ermee doen? (De "Microbellen")

Om te laten zien dat het werkt, hebben ze een microbelletje (een heel klein blaasje met gas erin, omhuld door een laagje) in hun virtuele water gegooid.

• Het experiment: Ze stuurden een ultrasone golf door het water.
• Het resultaat: Het microbellenletje begon te bewegen en zweefde naar een specifiek punt in de golf (een "drukknopen"). Dit heet akoestoforese.
• Waarom is dit cool? Dit is precies hoe artsen in de toekomst medicijnen kunnen sturen. Ze kunnen microbellen vullen met medicijn en met geluidsgolven precies naar een tumor sturen, waar ze dan openbarsten en het medicijn afgeven.

5. Waarom is dit belangrijk?

Vroeger was het heel moeilijk om te simuleren hoe geluid werkt op het niveau van cellen en kleine deeltjes.

• De "Virtuele Laboratorium": Met deze nieuwe machine kunnen onderzoekers nu duizenden experimenten doen in de computer voordat ze echt in het lab gaan werken.
• De voordelen:
  • Je hoeft geen dure apparatuur te kopen.
  • Je kunt gevaarlijke situaties veilig testen.
  • Je kunt zien wat er gebeurt op het niveau van cellen, wat je met een echte echo niet kunt zien.

Kortom: De onderzoekers hebben een nieuwe, snellere en stevigere manier gevonden om geluidsgolven in water te simuleren op een heel klein niveau. Het is alsof ze een nieuwe bril hebben ontworpen waardoor we kunnen zien hoe geluidsgolven met cellen en medicijntjes dansen, wat helpt bij het ontwikkelen van betere behandelingen voor ziektes.

Technische samenvatting

Probleemstelling

Ultrasone geluidsgolven (US) spelen een cruciale rol in biomedische toepassingen, variërend van beeldvorming tot therapeutische ingrepen zoals gerichte medicijndrang en tumorablatie. Het simuleren van de interactie tussen ultrasone golven en materie (vooral in zachte weefsels en biologische systemen) is echter een aanzienlijke uitdaging voor de computervakwetenschap.

• Schaalprobleem: Er bestaat een groot tijdsverschil tussen viskeuze en akoestische tijdschalen. Continuummethoden (zoals CFD) kunnen grootschalige golfvoortplanting goed modelleren, maar falen in het oplossen van microscopische interacties. Deeltjesgebaseerde methoden bieden moleculaire resolutie, maar kampen vaak met stabiliteits- en efficiëntieproblemen bij grotere schalen.
• Specifieke beperkingen: Bestaande deeltjesmethoden (zoals standaard DPD of SDPD) lijden vaak onder numerieke "bevriezing" (crystallisatie) bij hoge coarsening-niveaus of vertonen instabiliteiten (trekinstabiliteit) bij negatieve drukken, wat essentieel is voor het modelleren van geluidsgolven (die wisselen tussen compressie en expansie).
• Behoefte: Er is behoefte aan een "virtuele ultrasone machine" die zowel de viscositeit als de geluidssnelheid van water nauwkeurig kan reproduceren op micrometer-schaal (MHz tot GHz frequenties), zonder de complexiteit van hybride frameworks (die mesh- en deeltjesmethoden koppelen).

Methodologie

De auteurs introduceren een nieuwe, deeltjesgebaseerde methode genaamd usSDPD (ultrasound Smoothed Dissipative Particle Dynamics), die is geïmplementeerd in de LAMMPS-simulatiecode. De kern van de methode bestaat uit de volgende innovaties:

1. Implicit Compressible Pressure Solver:

   • Om de strenge tijdstapbeperkingen veroorzaakt door de hoge geluidssnelheid in water te omzeilen, wordt een impliciete solver voor drukkrachten gebruikt.
   • Dit maakt het mogelijk om tijdstappen tot 40 keer groter te nemen dan bij expliciete methoden, terwijl dissipatieve en willekeurige krachten expliciet blijven behandeld.
   • De solver lost de continuïteitsvergelijking iteratief op (via een relaxed Jacobi-methode) om de dichtheid binnen een nauwkeurigheid van 0,1% te houden.

2. Stabilisatie bij Negatieve Drukken (usSDPD):

   • Standaard SDPD lijdt onder "trekinstabiliteit" (tensile instability), waarbij het fluïdum voortijdig breekt onder negatieve druk (expansie).
   • Om dit op te lossen, hebben de auteurs twee aanpassingen gedaan:
     • Vermindering deeltjesdiameter: De deeltjesdiameter is verkleind van $0.5h$ naar $0.3h$ (waarbij $h$ de steunbreedte is), wat het aantal buren verhoogt en deeltjes in het "paar-instabiliteitsregime" brengt.
     • Lennard-Jones Potentiaal: Een extra afstotende Lennard-Jones (LJ) potentiaal is toegevoegd om te voorkomen dat deeltjes fysiek onrealistisch dicht bij elkaar komen (paar-instabiliteit), zonder de systematische eigenschappen te verstoren bij normale drukken.
     • Kunstmatige Druk (Artificial Pressure): Een term naar het voorbeeld van Monaghan is toegevoegd om de grootte van de drukkrachten tussen deeltjes te verminderen bij negatieve drukken. Dit verschuift het punt waarop trekinstabiliteit optreedt naar lagere drukken.

3. Virtuele Ultrasone Machine:

   • De simulatieomgeving bestaat uit een open simulatiebox met twee buffers aan de tegenoverliggende kanten die fungeren als ultrasone transducers.
   • Deze buffers worden gemodelleerd met Open Boundary Molecular Dynamics (OBMD), waardoor deeltjes kunnen worden toegevoegd of verwijderd om drukoscillaties te genereren en een staande golf in het gebied van belang (ROI) te creëren.

4. Toepassing op Microbellen:

   • De methode wordt gekoppeld aan een model van ingekapselde microbellen (EMB's). De buitenste vloeistof gebruikt usSDPD, terwijl de binnenste vloeistof van de bel (die veel comprimeerbaarder is) wordt gemodelleerd met standaard DPD. De membraan wordt beschreven als een isotroop elastisch model.

Belangrijkste Resultaten

• Stabiliteit en Nauwkeurigheid: De usSDPD-methode toont aanzienlijk betere stabiliteit dan standaard SDPD, vooral bij negatieve drukken. De druk-dichtheidsrelatie (EOS) blijft lineair zelfs bij hoge negatieve drukken, wat essentieel is voor het correct modelleren van geluidsgolven.
• Efficiëntie: Door de impliciete solver kan de tijdstap met een factor 40 worden vergroot, wat de rekentijd voor simulaties in het MHz- tot GHz-bereik drastisch verlaagt.
• Acoustoforese Simulatie: De auteurs simuleerden de acoustoforese (migratie door akoestische krachten) van een microbel in een staande golf.
  • De microbel migreerde correct naar de drukantinode (zoals theoretisch verwacht voor comprimeerbare deeltjes).
  • De migratiesnelheid schaalde lineair met de viscositeit en omgekeerd evenredig met het kwadraat van de drukamplitude ($\Delta p^{-2}$), in overeenstemming met theoretische modellen.
  • De simulaties werden uitgevoerd op schalen van 0,1 µm tot 10 µm en frequenties van 39 MHz tot 178 MHz.
• Modulariteit: Het framework toont aan dat verschillende vloeistofmodellen (usSDPD voor water, DPD voor de belinhoud) en structuren (elastische membranen) naadloos kunnen worden gekoppeld zonder complexe interpolatieschema's die nodig zijn bij mesh-gebaseerde methoden.

Bijdragen en Betekenis

• Eerste Deeltjesmethode voor US: Dit is de eerste deeltjesgebaseerde methode die zowel de viscositeit als de geluidssnelheid van water nauwkeurig kan reproduceren op micrometer-schaal, terwijl het elastische artefacten bij negatieve drukken minimaliseert.
• Virtueel Laboratorium: De "virtuele ultrasone machine" biedt een krachtig platform voor in silico experimenten. Het stelt onderzoekers in staat om een breed parameterbereik (frequentie, intensiteit, vloeistofeigenschappen) te verkennen zonder de beperkingen van fysieke apparatuur of de moeilijkheid om hoge frequenties experimenteel te genereren.
• Biomedische Toepassingen: De methode opent nieuwe wegen voor het bestuderen van ultrasone-gedreven fenomenen in complexe biologische systemen, zoals:
  • Gerichte medicijndrang via microbellen.
  • Interacties met celmembranen (sonoporation).
  • Dynamica van rode bloedcellen, biofilms en micro-robots.
• Toekomstperspectief: Hoewel de huidige versie isotherm is en geen bulkviscositeit bevat (wat demping beïnvloedt), vormt het een solide basis voor toekomstige uitbreidingen, zoals het modelleren van niet-reflecterende randvoorwaarden (voor lopende golven) en thermische effecten.

Kortom, dit paper introduceert een robuust en schaalbaar computermiddel dat de kloof overbrugt tussen microscopische deeltjesdynamica en macroscopische akoestische golven, wat essentieel is voor de ontwikkeling van geavanceerde ultrasone medische technologieën.