A new shape option — SUPELLIPS — allows automatic generation of targets in the "superellipsoid" family of shapes.

]]>I would like to surround two interacting spherical nanoparticles, each of a layer of silica; can i do this using the function ''CONELLIPS'', if it is possible how can i do that, if not which function can i use.

Thanks for your help! ]]>

Please i need to know if it is possible to have the imaginary part of the scattered wave with DDSCAT. Please if it is possible how can i have it.

thanks for your help ]]>

I have two questions:

1. I'm using the "nSphere" target generation to produce a target without overlapping spheres. I'd like to include nearfield calculations to examine if there's any nearfield interaction between the spheres. For a single sphere, such as the ellipsoid nearfield example, the nearfield calculation extends the computational volume beyond an individual sphere, as specified in the ddscat.par file. However, I'm confused how this works for nSpheres, where dipoles are only placed at sphere locations. Would a nearfield calculation with fractional extensions of 0.5 in all directions expand the computational volume around each sphere or would it somehow expand the computational volume around the entire TUC? Assuming its the former, if I wanted to include sphere-sphere nearfield interactions, should I specify the fractional extension to include nearest neighbor spheres? Or, am I misinterpreting these things?

2. I've been trying to extract the Qrad and Qtorque from calculations. When I include everything (i.e. "DOTORQUE"), the qtable still doesn't include those values. What setting should I double check to make ddscat write the Qrad and Qtorque?

Thank you for the help!

-Matt

I am trying to calculate the Qabs of a target immersed in water with wavelength ranged from 300 nm to 700 nm. Since the refractive index will change with the wavelength frequency, how should I set the refractive index in DDSCAT.par? It seems I could only set a special number instead of variables changing with wavelength. Thanks very much! ]]>

I want to investigate the properties using shape parameter ONIONSHEL for three layers with the spherical shape.

So, as given in DDSCAT manual on page no. 21, if SHPAR1 =1, then we have to define NCOMP = 3.

Now, problem is that how to set the thickness or radius of first, second and third layers with shape parameter SHPAR1, SHPAR2, and SHPAR3?

Please give the example with setting the values SHPAR1, SHPAR2, and SHPAR3 to run the ddscat.par file.

I wish you will reply soon and explain it with depth.

Thanking You

Your Sincerely

Bhatia P. ]]>

Hello, I am new and learning to use DDSCAT 7.3. I am at first trying to reproduce results from Dr. Flatau’s 2-spheres model that are in contact with each other (APPLIED OPTICS / Vol. 32, No. 18 / 20 June 1993 ). I have tried to reproduce the data, but with no luck. First of all, my S11 and S22 values are the same. Also, when I plot and compare the points to the literature values, the forward scattering looks like it reproduces S11, whereas the backscattering looks more like S22 character. Please let me know if you have suggestions for me.

Here is my input file, and header of my shape.dat file.

Thank you,

Yohanna White

*

shape.dat:

>TARELL ellipsoidal grain; AX,AY,AZ= 32.0000 32.0000 32.0000

34512 = NAT

1.000000 0.000000 0.000000 = A_1 vector

0.000000 1.000000 0.000000 = A_2 vector

1.000000 1.000000 1.000000 = lattice spacings (d_x,d_y,d_z)/d

-15.50000 0.50000 0.50000 = lattice offset x0(1-3) = (x_TF,y_TF,z_TF)/d for dipole 0 0 0

JA IX IY IZ ICOMP(x,y,z)

1 -2 -4 -16 1 1 1

2 -1 -4 -16 1 1 1

3 0 -4 -16 1 1 1

4 1 -4 -16 1 1 1

5 -3 -3 -16 1 1 1

ddscat.par:

' ========== Parameter file for v7.3 ==================='

' Preliminaries '

'NOTORQ' = CMTORQ*6 (DOTORQ, NOTORQ) — either do or skip torque calculations

'PBCGS2' = CMDSOL*6 (PBCGS2, PBCGST, GPBICG, QMRCCG, PETRKP) — CCG method

'FFTMKL' = CMETHD*6 (GPFAFT, FFTMKL) — FFT method

'GKDLDR' = CALPHA*6 (GKDLDR, LATTDR, FLTRCD) — DDA method

'NOTBIN' = CBINFLAG (ALLBIN, ORIBIN, NOTBIN)

' Initial Memory Allocation '

100 100 100 = dimensioning allowance for target generation

' Target Geometry and Composition '

'FROM_FILE' = CSHAPE*9 shape directive

no SHPAR parameters needed

2 = NCOMP = number of dielectric materials

'm1.33_0.001' = name of file containing RI info

'm1.33_0.001' = name of file containing RI info

' Additional Nearfield calculation? '

0 = NRFLD (=0 to skip nearfield calc., =1 to calculate nearfield E)

0.0 0.0 0.0 0.0 0.0 0.0 (fract. extens. of calc. vol. in -x,+x,-y,+y,-z,+z)

' Error tolerance '

1.00e-5 = TOL = MAX ALLOWED (NORM OF |G>=AC|E>-ACA|X>)/(NORM OF AC|E>)

' maximum number of iterations allowed '

300 = MXITER

' Interaction cutoff parameter for PBC calculations '

1.00e-5 = GAMMA (1e-2 is normal, 3e-3 for greater accuracy)

' Angular resolution for calculation of <cos>, etc. '

0.5 = ETASCA (number of angles is proportional to [(3+x)/ETASCA]^2 )

' Wavelengths (micron) '

6.283 6.2833 1 'LIN' = wavelengths (first, last, how many, how=LIN,INV,LOG,TAB)

' Refractive index of ambient medium'

1.000 = NAMBIENT

' Effective Radii (micron) '

10.000 10.000 1 'LIN' = eff. radii (first, last, how many, how=LIN,INV,LOG)

' Define Incident Polarizations '

(0,0) (1.,0.) (0.,0.) = Polarization state e01 (k along x axis)

2 = IORTH (=1 to do only pol. state e01; =2 to also do orth. pol. state)

' Specify which output files to write '

0 = IWRKSC (=0 to suppress, =1 to write ".sca" file for each target orient.

' Specify Target Rotations '

0. 0. 1 = BETAMI, BETAMX, NBETA (beta=rotation around a1)

0. 0. 1 = THETMI, THET MX, NTHETA (theta=angle between a1 and k)

0. 0. 1 = PHIMIN, PHIMAX, NPHI (phi=rotation angle of a1 around k)

' Specify first IWAV, IRAD, IORI (normally 0 0 0) '

0 0 0 = first IWAV, first IRAD, first IORI (0 0 0 to begin fresh)

' Select Elements of S_ij Matrix to Print '

2 = NSMELTS = number of elements of S_ij to print (not more than 9)

11 22 = indices ij of elements to print

' Specify Scattered Directions '

'LFRAME' = CMDFRM (LFRAME, TFRAME for Lab Frame or Target Frame)

2 = NPLANES = number of scattering planes

0. 0. 180. 5 = phi, thetan_min, thetan_max, dtheta (in degrees) for plane 1

90. 0. 180. 5 = phi, thetan_min, thetan_max, dtheta (in degrees) for plane 2

I have been trying to run DDSCAT for some custom shapes (FROM_FILE), with the shape.dat files generated from a MATLAB code or from the nanoHUB DDA Convert tool using Blender shapes. These are always created with the shape center at (0, 0, 0). The extinction spectra appear correct, and the E-field intensities generally make sense, but I run into issues somewhere between running ddpostprocess and using MayaVi or Paraview to visualize.

Despite the shape being centered, the E-field intensity map is always off-center and the shape outlines are thus cut off at the edges since they lay outside the computational volume (as ddpostprocess tells me when I try to select tracks from one edge of my shape to the other). I tried to visualize the target.out VTR files on the same set of axes in Paraview, and it looks like the entire shape is somehow thrown off-center since it matches up with the E-field plot. In MayaVi, however, I can only seem to visualize each separately, and while the E-field is still off-center, the target shape appears complete and centered in the outline box. This doesn't seem to change regardless of what I use to create the shape.dat file or what I use to visualize the ddpostprocess output.

While I know I can get a more complete image by changing the fractional extension of calculation volume in ddscat.par, my main concern is that the tracks I specify in ddpostprocess.par will not be directly through the center of the target as desired, since the entire shape is offset by some amount. Does anyone have any insight on how to resolve this issue? Would the output point towards some error in shape generation, in writing ddscat.par, or in writing ddpostprocess.par?

Thanks,

Rahil U.

FATAL ERROR IN PROCEDURE: TARCYL

NAT.GT.MXNAT

]]>I am using DDSCAT 7.3 to calculate the orientationally-averaged values of *Q*_{abs}, *Q*_{sca} and <cos*θ*> of randomly-oriented soot aggregates interacting with monochromatic plane waves with specified polarization states. According to * User Guide for the Discrete Dipole Approximation Code DDSCAT 7.3*, DDSCAT allows the user to specify a general elliptical polarization state for the incident radiation, by specifying the (complex) polarization vector ${\hat e_{01}}$ (see §24). Moreover, if

(1) Assume that the incident plane wave has a specified polarization state represented by the following Jones vector:

(1)\begin{align} J = \left[ {\begin{array}{*{20}{c}} {{E_{0y}}{e^{{\text{ - }}i{\delta _y}}}} \\ {{E_{0z}}{e^{{\text{ - }}i{\delta _z}}}} \end{array}} \right] \Leftrightarrow \left\{ {\begin{array}{*{20}{c}} {J = \frac{1}{{\sqrt {E_{0y}^2 + E_{0z}^2} }}\left[ {\begin{array}{*{20}{c}} {{E_{0y}}} \\ {{E_{0z}}{e^{{\text{ - }}i\delta }}} \end{array}} \right]} \\ {\delta = {\delta _z} - {\delta _y}} \\ {\cos \beta = \frac{{{E_{0y}}}}{{\sqrt {E_{0y}^2 + E_{0z}^2} }}} \\ {s{\text{in}}\beta = \frac{{{E_{0z}}}}{{\sqrt {E_{0y}^2 + E_{0z}^2} }}} \end{array}} \right\} \Rightarrow J = \left[ {\begin{array}{*{20}{c}} {\cos \beta } \\ {s{\text{in}}\beta \cos \delta - is{\text{in}}\beta s{\text{in}}\delta } \end{array}} \right] \end{align}

where *E*_{0y} and *E*_{0z} are the amplitude of the two components of the electric field vector, while *δ*_{y} and *δ*_{z} are the phase of the two components of the electric field vector, **should we simply set ${\hat e_{01}}$ by specifying (0, 0) (cos β, 0) (sinβcosδ, - sinβsinδ) in ddscat.par ?**

(2) Assume that we simply set ${\hat e_{01}} = {\hat y_{LF}}$ and ${\hat e_{02}} = {\hat z_{LF}}$, and calculate the orientationally-averaged values of *Q*_{abs}, *Q*_{sca} and <cos*θ*> in the two orthonormal polarization states, **can we calculate the orientationally-averaged values of Q_{abs}, Q_{sca} and <cosθ> in an arbitrary incident polarization state, just based on the orientationally-averaged values of Q_{abs}, Q_{sca} and <cosθ> in the two orthonormal polarization states ?**

Qext | Qabs | Qsca | g(1)=<cos> | <cos^2> | Qbk | Qpha | |

JO=1: | 9.2037E-01 | 8.1852E-01 | 1.0185E-01 | 7.2936E-02 | 4.0114E-01 | 1.0141E-02 | 5.3381E-01 |

JO=2: | 9.2037E-01 | 8.1852E-01 | 1.0185E-01 | 7.2936E-02 | 4.0114E-01 | 1.0141E-02 | 5.3380E-01 |

I am wondering **if the data in the lines of JO=1 and JO=2 in the red box are the calculation results of Q_{abs}, Q_{sca} and <cosθ> in the incident polarization states ${\hat e_{01}}$ and ${\hat e_{02}}$, respectively ?**

It is highly appreciated if you would instruct me on these questions, and I am looking forward to your reply.

Sincerely yours

Ya-fei Wang

Thanks in advance!

]]>I have modelled a group of nanospheres grouped together almost in a straight line. Can I model the scattering properties such that results obtained from the first nanosphere would become an input for the second nanosphere and so on for the whole network?

Thanks,

Hirak ]]>

I'm running the code in Linux Mint platform. The computer has 16GB RAM, 8 cores & I also installed 'OpenMPI.' Please kindly help me.

Thanks

]]>I am doing the near-field calculations, as given in user guide of DDSCAT 7.3 on 29, 29.1 & 29.2 sections. when, we calculate the near-field calculation after choosing important parameter i.e. 1 = NRFLD (=0 to skip nearfield calc., =1 for nearfield E, =2 for E and B) and 2 = IORTH (=1 to do only pol. state e01; =2 to also do orth. pol. state), and 1 = IWRKSC (=0 to suppress, =1 to write ".sca" file for each target orient.; After doing these changing in ddscat.par file and run in the ddscat. we got the extra file like w000r000k000.E1 & w000r000k000.E2 and so on….w000r000k000.E200. Till now, i followed these steps.

Now, my Question is how to use .E2, .E2, .E3 and so on .E200 files for further calculation ?

Q1. Is it single point calculation.

Q2. If it is not single point calculation, then how to add these file into E1+E2…….=En

Q3. is it run directly through ddpostprocess or we have to make ddpostprocess.par file. if yes, then How to make and which parameter are use in ddpostprocess.par file. and how to run after preparing the ddpostprocess.par file.

I wish you will reply soon.

Thanking You ]]>

The testing problem is the scattering problem of a sphere: the geometry parameters of the target is

'ELLIPSOID'

80 80 80

Then I set OMP_NUM_THREADS=12 and run the executable ddscat.

Here is the problem: it seems the program is not running in openmp mode. The command "top d1" shows only about 200%, instead of 1200%, of CPU is used. No matter what value is given to OMP_NUM_THREADS, the usage of CPU is always about 200%.

Did I miss anything in compiling, or running the program ?

Thank you very much!

]]>To do this I put the ddscat program, and ddscat.par file in a new directory. When I run the program it generates the mtable, qtable, qtable2 and target.out files, but not any .avg, .fml or .sca files. The qtable extinction data is empty. When executing, the program appears to make it through the 100 iterations, but then closes out after. The problem seems to be related to the aeff and/or the shape parameters. When I make the aeff too small or the shape parameters too large, then I get the problem as I described. I have been trying to keep the shape parameters large enough to have ~10000 dipoles. But if I set the aeff to match my experimental samples, then the program doesn't generate extinction data (it does seem to work for much larger aeff sizes).

Below is an example of a .par file that I have been using for a gold rod with dimensions L=37 nm (including end caps) and d=7 nm. Therefore the aeff for this geometry is 0.0068 microns. I tried to make the shape parameters a multiple of 2 (60 x 14), to get over 10000 dipoles.

I have tried this on the Windows precompiled version, as well as on a Linux system with the latest version. I would be tremendously grateful if you could give me any insight on what I might be doing wrong.

' ========== Parameter file for v7.3 ==================='

' Preliminaries '

'NOTORQ' = CMDTRQ*6 (DOTORQ, NOTORQ) — either do or skip torque calculations

'PBCGS2' = CMDSOL*6 (PBCGS2, PBCGST, GPBICG, QMRCCG, PETRKP) — CCG method

'GPFAFT' = CMETHD*6 (GPFAFT, FFTMKL) — FFT method

'GKDLDR' = CALPHA*6 (GKDLDR, LATTDR, FLTRCD) — DDA method

'NOTBIN' = CBINFLAG (NOTBIN, ORIBIN, ALLBIN) — binary output?

' Initial Memory Allocation '

100 100 100 = dimensioning allowance for target generation

' Target Geometry and Composition '

'CYLNDRCAP' = CSHAPE*9 shape directive

60 14 = shape parameters 1 - 2

1 = NCOMP = number of dielectric materials

'../diel/Au_evap' = file with refractive index 1

' Additional Nearfield calculation? '

0 = NRFLD (=0 to skip nearfield calc., =1 to calculate nearfield E)

0.0 0.0 0.0 0.0 0.0 0.0 (fract. extens. of calc. vol. in -x,+x,-y,+y,-z,+z)

' Error Tolerance '

1.00e-5 = TOL = MAX ALLOWED (NORM OF |G>=AC|E>-ACA|X>)/(NORM OF AC|E>)

' Maximum number of iterations '

100 = MXITER

' Integration limiter for PBC calculations '

1.00e-2 = GAMMA (1e-2 is normal, 3e-3 for greater accuracy)

' Angular resolution for calculation of <cos>, etc. '

0.5 = ETASCA (number of angles is proportional to [(3+x)/ETASCA]^2 )

' Wavelengths (micron) '

0.650 0.650 1 'LIN' = wavelengths (1st,last,howmany,how=LIN,INV,LOG,TAB)

' Refractive index of ambient medium '

1.333 = NAMBIENT

' Effective Radii (micron) '

0.0068 0.0068 1 'LIN' = eff. radii (1st,last,howmany,how=LIN,INV,LOG,TAB)

' Define Incident Polarizations '

(0,0) (1.,0.) (0.,0.) = Polarization state e01 (k along x axis)

1 = IORTH (=1 to do only pol. state e01; =2 to also do orth. pol. state)

' Specify which output files to write '

0 = IWRKSC (=0 to suppress, =1 to write ".sca" file for each target orient.

' Specify Target Rotations '

0. 0. 1 = BETAMI, BETAMX, NBETA (beta=rotation around a1)

0. 0. 1 = THETMI, THETMX, NTHETA (theta=angle between a1 and k)

0. 0. 1 = PHIMIN, PHIMAX, NPHI (phi=rotation angle of a1 around k)

' Specify first IWAV, IRAD, IORI (normally 0 0 0) '

0 0 0 = first IWAV, first IRAD, first IORI (0 0 0 to begin fresh)

' Select Elements of S_ij Matrix to Print '

9 = NSMELTS = number of elements of S_ij to print (not more than 9)

11 12 21 22 31 33 44 34 43 = indices ij of elements to print

' Specify Scattered Directions '

'LFRAME' = CMDFRM (LFRAME, TFRAME for Lab Frame or Target Frame)

1 = NPLANES = number of scattering planes

0. 0. 180. 1 = phi, theta_min, theta_max (deg) for plane A

I want calculate nearfield of 2 shape , I see 2 different errors for each shape,one error when run ddpostprossess is " readnf sanity failure: inconsistent j1= 245786 and nat0= 279320"

and another error when run ddscat.par of another shape :"

Program received signal SIGSEGV: Segmentation fault - invalid memory reference.

Backtrace for this error:

#0 0x7F354D819777

#1 0x7F354D819D7E

#2 0x7F354D16BD3F

#3 0x7F354D1B7E2C

#4 0x401DB9 in ddscat_

#5 0x40E2B5 in MAIN__ at DDSCAT_cpp.f90:? "

I dont know meaning of 2 errors?

Thanks advance. ]]>

DDA solution for polarization considering a gold nanorod shows an interesting pattern. The well-documented pattern (i.e. maximum @90°and minimum @180°) was obtained when the transverse wavelength was considered and beta, theta, and phi were set to zero. Using the same settings, the polarization patterns obtained for the longitudinal wavelength show patterns which deviate from their physical picture.

Please suggest the required modification that I need to make for obtaining the correct polarization pattern.

Thanks and regards,

Hirak ]]>

I want to calculate the polarization of two spheres.

One sphere has 1.33+0.01i refractive index and another sphere has 1.9+0.45i.

So I set ddscat.par such as :

2 = NCOMP = number of dielectric materials

'../diel/m1.33_0.01'

'../diel/m1.90_0.45'

Then, m1.33+0.01i is number 1 and m1.90++0.45i is nuber 2. Is it right?

And I set BAM file such as :

x y z a IC1 IC2 IC3 theta_DF phi_DF beta_DF

0.000000 0.000000 0.000000 0.500000 1 1 1 0.00000 0.00000 0.00000

0.000000 1.000000 0.000000 0.500000 2 2 2 0.00000 0.00000 0.00000

Because a sphere at (0,0,0) has m=1.33+0.01i which is number 1, so I set IC1=IC2=IC3=1.

And also at a sphere at (0,1,0) , I set IC1=IC2=IC3=2 because m=1.90++0.45i is nuber 2.

Is it right?

I can't sure that means of number of IC1, IC2, IC3.

I am a new user of ddscat.

I try to calculate scattering and absorption efficiencies of metal nanoparticles of arbitrary shapes and sizes.

I did a serial calculation with a sphere easily. but when I try a parallel calculation with MPI, I find very high values.

It is the same for all the algorithms

ZBCG2 IT= 1 f.err= NaN

…

cgcommon ckpt 20

iter= 8 frac.err= NaN

…

GPBICG IT= 10 f.err= 6.167E+00

…

…

sqrt(rnorm/bnorm)= 3 174106950333229.

sqrt(rnorm/bnorm)= 4 NaN

Other codes using MPI have been successfully implemented in our system, and the job submission procedures are correctly implemented.

The administrator followed all the instructions in the USer Guide. Did we miss an important information?

Regards

]]>I am using DDA to calculate the optical properties such as absorption, scattering and extinction for "CONELLIPS" shape (Two -layers). I have no idea about CONELLIPS shape, this type of target geometry is in-built or not. I have only changed the target geometry-(ELLIPSOID to CONELLIPS) and the composition-(40 40 40 60 60 60 = shape parameters 1 - 6) and run the code for two different materials. After that, i got the error as given below:

Error: FATAL ERROR IN PROCEDURE: TARCEL

AX < BX

EXECUTION ABORTED ]]>

I am working with ellipsoid Nanomatryushka particles having two or more layers. I need to simulate values for the Muller Matrices for the kind of particles. Can I do it with DDSCAT? Can I model Fano Resonance with it?

If I need to modify the CONELLIPSE section for my use, how can I do it?

Regards,

Hirak ]]>

Thank you

]]>I am work with ddscat7.2. Now I want plot polarization direction in each point.

but I cant compile wxxxryyykzzz.poln .Could you help me.

Thank you in advance.

I study scattering using DDSCAT7.3.

I have some questions about CYLNDRPBC.

I want to make regolith using CYLNDRPBC, so I wrote shape parameter such like ::60.0 1.0 1 1.0 1.0 .

Then, in the output file 'w000r000k000.sca', there are only two theta values(0, 180)

I want to (0~180)theta value.

How can I do?

I am research scholar in IIT Delhi, India. My research work is related to find the optical properties of the Nano-particles of different materials and different geometries. I found that the DDSCAT 7.3 is one of the most prominent code to study the characteristically properties of the different geometries. I am highly interested to find the optical properties for such geometries like ANI_ELL_2 ELLIPSO_2 and SPHROID_2. I have downloaded the user guide for DDSCAT 7.3 and its related examples from site of DDSCAT, But I did not found any example related to ANI_ELL_2, ELLIPSO_2 and SPHROID_2geometries.

Could you please provide the us these Examples (code) with their supporting files. The files will help me to study the optical properties of above geometries, which is very useful to enhance the efficiency of photovoltaics.

I will be highly obliged to your favor.

Thanking you

With Warm Regards

Shivani Bhardwaj

]]>Hi Dr. Draine,

We have decided to use the Fortran version of DDSCAT, rather than the C++

version; yours is much more up to date, reliable, and developed (probably

because DDSCAT seems to be your idea and creation in the first place). We

are using it on a high-performance computer because of the processing

requirements.

The DDSCAT program is incredible; the development of it over the years is

proving very useful for people like us using it in independent research.

We are currently trying to use the program to model simple reflection and

refraction, where we convert the incident light as a Gaussian beam into

Electric field components to simulate a plane wave, then convert the

output E-field back into a wave. We are trying to make our own code to generate

visualization graphics; and we believe that this method will allow us to

use DDSCAT for twisted light rather than regular plane waves.

I was hoping you could speak a little to the meaning of the outputs: I do

believe that the wxxxryyykzzz.En binary output is the one we want to use

for our calculations. Are you defining the near field as such because it is

the sum of the incident field and the scattered field? I am a little bit

confused on this; most things discussing near-field and far-field

radiation seem to indicate that there may be some other important difference between

the two fields.

So… I am hoping to clarify that E_sca is indeed the E field that we will

be able to process in order to create our scattered light waveforms. If you

can clarify this for me, that would be awesome. Otherwise, I'm sure with

the help of our professor my group will be able to figure out a solution.

Thank you very much for your time.

Best,

Josh Graber

We have the latest release of ddscat up and running on the high

performance computing center here; however we need the ASCII files as

ddpostprocess.out is said to contain.

The ddposprocess.out file doesn't show up though, and the ddscat.log

file says

"file = w000r000k000.E1

file is incompatible with present version of subroutine SUBREADNF"

We think that might be our problem. Could you possibly speak to that? I

don't know if it's something wrong with the way it was set up on our

system of if there is some new update or an error in the subroutine, but

anything you could give us would help a lot.

Thanks,

Josh Graber

FATAL ERROR IN PROCEDURE: ZBCG2WP

ITERN>ITERMX

EXECUTION ABORTED

How to solve this problem? Thanks a lot!

]]>»»> ITERN>ITERMX

»»> EXECUTION ABORTED

How to solve this problem?

]]>Hi Suman Dhayaj

We have Implemented FFT near-field. One version of this paper is available here:

ftp://ftp.astro.princeton.edu/draine/papers/pdf/OptEx_20_1247.pdf

But if you download DDSCAT code from

http://www.ddscat.com

you can also look at our implementation. It uses FFT but if do not have many dipoles you can calculate

near field witout FFT technique. We had a simple version of the near field calculation available in the

DDSCAT release before the 2012 (we can point to it if you need).

It is difficult for us to study your implementation because of time constraints on our part.

Piotr

On Fri, Nov 18, 2016 at 11:45 AM, Suman Dhayal <moc.liamg|layahdss#moc.liamg|layahdss> wrote:

Dear Professor Draine, Flatau,

This is Ms Suman Dhayal, an associate in research at Duke University. I need some help with the calculations using discrete-dipole approximation.

I have written the codes for calculations using DDA following the description in your JOSA, 11, 1994 paper. I have used it to do some Raman scattering calculations, and I am interested in doing some near field calculations now. I looked through your codes user's guide and some other papers but perhaps I am missing some information. If you could help me with it, I will be grateful to you.

So my query is over whether it will be correct to just take the unknown polarizations obtained by solving the simultaneous equations (7) from your paper and plug them in the general dipole field formula to compute the field in and around a sphere containing a cluster of spheres? I am attaching the small program I have written to do this if you would like to have a look at it.

It will be a big help if you could say whether I am doing it correctly. If not, some suggestion over where I am going wrong and how I may be able to correct it.

Thank you so much for your help,

Suman

]]>Yours truly,

DDSCAT user

Can we generate two infinite cylinders by using a periodic target?

]]>the following ways

- 7.3.2 fixes a bug that interfered with use of MPI for nearfield (NRFLD=1 or 2) calculations in 7.3.1
- 7.3.2 offers two new target options: regular octahedron, and spherical shell composed of uniaxial material (e.g. graphite)

The changes recommended by Dr. Ruzankin have now been implemented in a new release, DDSCAT 7.3.2, released today (2016 August 25) and available for download from ddscat.org — go to the "Downloads" tab. At this point in time we consider DDSCAT 7.3.2 to be fully compatible with MPI. If users encounter bugs with use under MPI, they should report them here.

]]>What is the largest applicable target size compared with wavelength? Any experience to calculate ~100um diameter target with 532nm wavelength?

]]>I have changed the scattering direction from TFRAME to LFRAME for one of example ddscat.par, CYLNDRPBC. Then, this error occurs;

FATAL ERROR IN PROCEDURE: REAPAR

CAN NOT USE LFRAME WHEN JPBC=1 OR 2

EXECUTION ABORTED

What is 'JPBC' ?

]]>(from DDPOSTPROCESS.f90)

! CXEINC(1-NXYZ,3) = complex incident macroscopic E field at all points

! CXESCA(1-NXYZ,3) = complex radiated macroscopic E field at all points

! CXPOL(1-NXYZ,3) = complex polarization/d^3 at all points

! CXADIA(1-NXYZ,3) = diagonal element of polarizability/d^3 at all pts

Can you please explain how each of these outputs are normalized? Specifically, can you please help me understand how to convert scattered electric field, incident electric field and polarization values extracted in DDSCAT to SI units?

]]>The best example I can give to help troubleshoot my issue is this ellipsoid (ddscat.par file copied at the bottom). I intended to create an ellipsoid with dimensions 7x7x11 nm and approximately 10,000 dipoles. I want the Qext values over a wavelength range of 200-1000 nm. I have tried this file for theta=0, theta=90, theta=0-90 (with 10 angles), as well as 0-90 for beta and theta, all while holding all the other angles constant at 0.. and every time I get the exact same Qext values. Shouldn't I get different extinction values when I am looking at the ellipsoid head on versus when it is rotated 90 degrees so I am looking at the long edge? Do you notice anything I am missing or doing wrong or do you have any suggestions of what I might try to change??

I've tried making this larger (70x70x110) and I've tried with as many as 100,000 dipoles as well, never with any success at recognizing the longer edge. I am using version 7.3.1, but not the most recent version of it.

Thank you!!

' ========== Parameter file for v7.3 ==================='

' Preliminaries '

'NOTORQ' = CMDTRQ*6 (DOTORQ, NOTORQ) — either do or skip torque calculations

'PBCGS2' = CMDSOL*6 (PBCGS2, PBCGST, GPBICG, QMRCCG, PETRKP) — CCG method

'GPFAFT' = CMETHD*6 (GPFAFT, FFTMKL) — FFT method

'GKDLDR' = CALPHA*6 (GKDLDR, LATTDR, FLTRCD) — DDA method

'NOTBIN' = CBINFLAG (NOTBIN, ORIBIN, ALLBIN) — binary output?

' Initial Memory Allocation '

100 100 100 = dimensioning allowance for target generation

' Target Geometry and Composition '

'ELLIPSOID' = CSHAPE*9 shape directive

25 25 39.286 = shape parameters 1 - 3

1 = NCOMP = number of dielectric materials

'../diel/polystyeq' = file with refractive index 1

' Additional Nearfield calculation? '

0 = NRFLD (=0 to skip nearfield calc., =1 to calculate nearfield E)

0.0 0.0 0.0 0.0 0.0 0.0 (fract. extens. of calc. vol. in -x,+x,-y,+y,-z,+z)

' Error Tolerance '

1.00e-5 = TOL = MAX ALLOWED (NORM OF |G>=AC|E>-ACA|X>)/(NORM OF AC|E>)

' Maximum number of iterations '

1000 = MXITER

' Integration limiter for PBC calculations '

1.00e-2 = GAMMA (1e-2 is normal, 3e-3 for greater accuracy)

' Angular resolution for calculation of <cos>, etc. '

0.5 = ETASCA (number of angles is proportional to [(3+x)/ETASCA]^2 )

' Wavelengths (micron) '

200 1000 801 'LIN' = wavelengths (1st,last,howmany,how=LIN,INV,LOG,TAB)

' Refractive index of ambient medium '

1.0000 = NAMBIENT

' Effective Radii (micron) '

4.0691 4.0691 1 'LIN' = eff. radii (1st,last,howmany,how=LIN,INV,LOG,TAB)

' Define Incident Polarizations '

(0,0) (1.,0.) (0.,0.) = Polarization state e01 (k along x axis)

2 = IORTH (=1 to do only pol. state e01; =2 to also do orth. pol. state)

' Specify which output files to write '

0 = IWRKSC (=0 to suppress, =1 to write ".sca" file for each target orient.

' Specify Target Rotations '

0. 0. 1 = BETAMI, BETAMX, NBETA (beta=rotation around a1)

0. 90. 10 = THETMI, THETMX, NTHETA (theta=angle between a1 and k)

0. 0. 1 = PHIMIN, PHIMAX, NPHI (phi=rotation angle of a1 around k)

' Specify first IWAV, IRAD, IORI (normally 0 0 0) '

0 0 0 = first IWAV, first IRAD, first IORI (0 0 0 to begin fresh)

' Select Elements of S_ij Matrix to Print '

9 = NSMELTS = number of elements of S_ij to print (not more than 9)

11 12 21 22 31 33 44 34 43 = indices ij of elements to print

' Specify Scattered Directions '

'LFRAME' = CMDFRM (LFRAME, TFRAME for Lab Frame or Target Frame)

1 = NPLANES = number of scattering planes

0. 0. 180. 1 = phi, theta_min, theta_max (deg) for plane A

In version 7.3.1's code (version 160601), in reapar.f90:

1470 ! check that user has not requested more than 1000 orientations and

1471 ! IWRKSC=1

1472

1473 IF(IWRKSC>0.AND.(NBETA*NTHETA*NPHI)>1000)THEN

1474 CALL ERRMSG('FATAL','REAPAR', &

1475 'error: if IWRKSC=1, NBETA*NTHETA*NPHI must be .le. 1000')

1476 ENDIF

Given the changes to the code in 7.3, this can be removed.

]]>In addition, a problem with parameter input for shape option BSLINPBC has been corrected. ]]>