Show Posts

This section allows you to view all posts made by this member. Note that you can only see posts made in areas you currently have access to.


Messages - DSarkar

Pages: [1] 2
1
Hi,

I am trying to run multiple MD simulations at the same time on multiple nodes of our cluster. Could you please let me know how the license checking takes place when I run the simulation? Is the license checked at the beginning of execution of each python script, or does it check continually at some interval, or at the beginning of each module (such as Geometry Optimization, different Molecular Dynamics modules, etc.)? It would help me define the pbs script accordingly. In case it is checked only once before the start of a python file, do you think it may be possible for me to define some "clones" of the atkpython file and separately set the path to each license file corresponding to the particular node (I have obtained different academic license for each node)?

Regards,

Debarghya

2
General Questions / Re: User Defined Tabulated Potential
« on: May 4, 2016, 17:15 »
Hi Julian and Jess,

Thanks for your replies. I ran the attached script and got the results as well. However, I am getting some error while running my script, which I cannot debug. I would be grateful if you could shine some light on it. (I have the tabulated potentials as 60000 points for each pair, 6 in total, and delta_r=0.0001 Angstrom)

Regards,

Debarghya

3
General Questions / User Defined Tabulated Potential
« on: May 4, 2016, 05:47 »
Hi,

What is the starting radius of the user defined tabulated potential? Is it 0 by definition?

Regards,

Debarghya

4
General Questions / Non-equilibrium momentum exchange
« on: April 29, 2016, 19:18 »
Hi,

I am trying to simulate thermal transport in VNL using the non-equilibrium momentum exchange formalism, and had a few clarifications regarding that:
1. In the MD Analyzer, the temperature profile that is shown, is that the average temperature at each point between tmin and tmax?
2. How exactly is the heat rate reported at the end of the log file, calculated? Is it also an average heat rate over all the exchange processes in the simulation? Or is it also averaged over all the simulation steps (so that essentially higher the simulation steps, lower the heat rate)?
3. Do you have a suggestion for a "good" heat flux? Or could you suggest a method for determining it based on maybe material parameters?
4. In a paper http://www.tandfonline.com/doi/pdf/10.1080/0026897031000068578, which also uses the same formalism, it is stated that atoms at the interface are loosely bound by weak harmonic springs to avoid dissolution. Can that be done in VNL?
5. Does this simulation consider atomic vibrational energy, or is it only translational energy that is considered?
6. Any system, irrespective of whether it has an interface, would come to the same temperature throughout if allowed sufficient time to relax. So instead of calling temperature profile at equilibrium, I would prefer calling it temperature profile at steady state, where the heat exchange rate and the time of relaxation are consistently chosen. It seems, therefore, that the choice of this combination is very critical. Do you have any suggestion regarding this, or any rule of thumb to follow based on material used, cross-section area, etc.?
7. Do you have any suggestion about the length of the "Hot" and "Cold" regions? I believe this also has a bearing on the resultant heat flux, since a longer length would possibly give a broader distribution of velocities and therefore, higher flux--is that understanding correct? Apart from that, what else needs to be considered for these regions? Once again, could you give a good number of atoms that you think should be included in these regions?
8. Some papers have set the boundary condition as temperature, and have calculated the heat flux. Is there a way in VNL to do so, i.e. have thermostats connected to some of the atoms, and see the evolution?
9. If I am interested in changing the exchange rate over time, which cannot be currently done directly in VNL, is it possible to carry over the results of one momentum exchange simulation to another? I know I can carry over the bulk configuration and velocities, but at this point I am a bit confused if that should be sufficient for the next momentum exchange simulation I would be looking for. Do you think I need to read some other results from the previous simulation and use that for some calculation in the next momentum exchange step?
10. Any other general comments inspired by these questions are highly welcome.

As a note in case it helps, I am currently looking at a Si-Ge-Si system where the temperature is above the melting point of Ge but below that of Si. I highly appreciate your comments. I have greatly benefited from all your previous replies.

Regards,

Debarghya

5
General Questions / Re: ZBL Potential
« on: April 27, 2016, 19:11 »
Hi Julian,

Thanks a lot for the detailed reply. I shall follow your suggestion to implement the potential. In general, do you think it would be a valid approach if I implemented the ZBL using General1Potential defined in VNL, using the correct parameters as given in the paper?

Regards,

Debarghya

6
General Questions / ZBL Potential
« on: April 19, 2016, 05:30 »
Hi,

I was trying to reproduce the results of this paper: http://www.tandfonline.com/doi/abs/10.1080/10420150.2013.792818#.VxWj-zArKhc. They have used a combination of Coulomb potential, Buckingham potential, and ZBL potential for the interactions. Referring to the tutorial http://docs.quantumwise.com/tutorials/combining_potentials.html, I have some fair idea of implementing the Coulomb and the Buckingham potentials. However, I have a few questions regarding the ZBL potential.
1. As per the description in http://www.quantumwise.com/documents/manuals/latest/ReferenceManual/index.html/ref.tersoffzblpotential.html, and looking through external references, I don't understand what bf and rf signify, and what values to put in for my structure.
2. How to set the outer cut-off radius (referred to as r_cut in the several other potential implementations)?
3. If I set the r_cut for ZBL and r_i for Buckingham, will the range of radii between them be automatically spline interpolated?
4. Any suggestion for the type of Coulomb interaction to use: DSF or SPME?
I greatly appreciate your useful replies.

Regards,

Debarghya

7
Yes, it does. Thank you.

8
Hi,

As described in several tutorials, we can extract the position and velocity of atoms in an ensemble from any time point in a simulation run by opening the "xyz.nc" file using the Movie Tool, and sending the configuration to the required module (builder, script generator, etc.). While this is a very convenient method for many applications, for relatively large system and/or simulation sizes, this method is quite time consuming. So I was wondering if it is possible to do either of the following:
1. Extract the bulk structure using some script at any time step without opening the "xyz.nc" file in the Movie Tool.
2. Open the file in the Movie Tool and run some script to select the bulk structure at, say the time step currently selected there, and send it to one of the modules.
Any suggestion/advice is highly appreciated.

Regards,

Debarghya

9
 Hi Julian,

Thanks for the detailed explanation. This understanding would definitely help me a lot in designing future simulations. I had a question about the effect of changing the coupling constant of the thermostat and barostat. In an attempt to understand the machinery better, I made a few sweeps of the parameters and have attached a small snapshot here due to file size constraints.
My observations are (all of which may not be obvious from this slide):
1. Thermostat time constant has the dominant effect in determining the envelope of KE /temperature variation dynamics.
2. Barostat time constant has its effect on the detailed fluctuations.
3. At higher coupling constants like >25 fs, a time step of 1 fs is sufficient to resolve the dynamics (since 0.5 fs timestep gives the same temperature profile). 4. But at lower coupling constants like 5 fs, there is a significant difference between 1fs timestep and 0.5 fs timestep.
My queries are:
1. Does the discrepancy between observations 3 and 4 above appear because the external forces "dominate" over the natural dynamics of the system and/or the external coupling is so fast that we need higher sampling to resolve it?
2. Based on all observations above, what would be an optimum coupling constant and time step for the system? What would be a standard procedure to determine these parameters properly in general for any system?
3. How can we know quantitatively that the system is relaxed? In other words, what command would give us the force on each atom which can be minimized?

Thanks a lot for your help.

Regards,

Debarghya

10
Hi Julian,

Thanks a lot for the script and the explanation.

Debarghya

11
General Questions / Re: Resource reqirement
« on: March 30, 2016, 16:22 »
Hi Ulrik,

Thanks for letting me know.
With respect to your comment: "If your calculation involves geometry optimization, you can use that the time required for each force evaluation step is approximately the same for a given system." How can we know the time taken for each step? Could you please elaborate a little more?

Regards,

Debarghya

12
Hi Jess.

Thanks a lot for your suggestions. One thing that occurred to me was whether the wavy nature of the system and an oscillatory behavior is present because my time step is not sufficiently small to resolve the highest vibrational frequencies of the atoms. In that case, I would want to find the vibrational frequencies and set the time step, say T=0.01*smallest vibrational frequencies. Do you think this is a valid approach? Also, how can this be done in VNL?

Regards,

Debarghya

13
Hi,

I am trying to simulate the interfacial thermal resistance at Si-Ge interfaces, and was following the tutorial on the Si(100) and Si(110) interface using MD. I have attached a file describing the geometry I am currently looking at. From the NPT stabilization of this and similar other previous simulations, it seems to me that the structure is not fully relaxed. In the NPT simulation, there appears to be large destabilizing force at the beginning, which results in weird situations like a wavy character throughout the Si-Ge-Si structure, premature melting at localized regions, etc. In this simulation I set the maximum number of steps in Geometry Optimization as 100000 (maximum possible in VNL), but it ran for 14 fs, and went onto the NPT step. Could you please help me setup the proper relaxation method for this kind of interface? Is there a suitable tutorial or other resource that you can possibly point me to, so that I can get a better hold on the geometry optimization methods? Thanks a lot, I really appreciate your help.

Regards,

Debarghya


14
General Questions / Resource reqirement
« on: March 25, 2016, 00:06 »
Hi,

Is there a way to estimate the time required to finish a simulation given the computation resources for a certain simulation setup? It would be helpful for me to estimate the time required as I am running on a common queued cluster, so that I can request resources optimally.

Regards,

Debarghya


15
General Questions / Re: Boundary conditions
« on: March 24, 2016, 08:30 »
Thanks for the clarification, Jess.

Pages: [1] 2