[DFTB-Plus-User] Mulliken analysis along the trajectory

Lyu, Dongyu dlyu at constructor.university
Wed Feb 28 18:20:38 CET 2024


Dear Bálint,

Thank you very much!
And are the last three columns from the xyz output velocities?
With the Electrostatic Potential option, I can only get the ESP values for fixed points along the trajectory, is that correct? In another words, it's impossible to change the points according to the atomic positions during the simulation?

Best regards,
Dongyu

获取 Outlook for iOS<https://aka.ms/o0ukef>
________________________________
发件人: DFTB-Plus-User <dftb-plus-user-bounces at mailman.zfn.uni-bremen.de> 代表 Bálint Aradi <aradi at uni-bremen.de>
发送时间: Wednesday, February 28, 2024 5:17:49 PM
收件人: User list for DFTB+ related questions <dftb-plus-user at mailman.zfn.uni-bremen.de>
主题: Re: [DFTB-Plus-User] Mulliken analysis along the trajectory

[You don't often get email from aradi at uni-bremen.de. Learn why this is important at https://aka.ms/LearnAboutSenderIdentification ]

Dear Dongyu,

sorry, to be more precise, the xyz output contains the mulliken
population, not the gross charge, so you would have to substract the nr.
of electrons in the neutral state.

As for the ESP, we have this quantity, but we do not print it currently.
It should not be compplicated, but we did not implement it yet. As a
workaround for the moment, you use the ElectrostaticPotential option,
and specify the atomic positions as Points.

Best regards,

Bálint


On 28.02.24 16:54, Lyu, Dongyu wrote:
> Dear Bálint,
>
> Thanks for your reply.
> For the charges, why the printed charges in the xyz output file for the
> last frame are different from the atomic gross charges written in the
> detailed.out? Are they the same property?
> Regarding EPS, Is it possible to extract the value of the electrostatic
> potential at the position of the atom for each frame during the simulation?
>
> Best regards,
> Dongyu
>
>
> ------------------------------------------------------------------------
> *From:* DFTB-Plus-User
> <dftb-plus-user-bounces at mailman.zfn.uni-bremen.de> on behalf of Bálint
> Aradi <aradi at uni-bremen.de>
> *Sent:* Wednesday, February 28, 2024 4:14 PM
> *To:* dftb-plus-user at mailman.zfn.uni-bremen.de
> <dftb-plus-user at mailman.zfn.uni-bremen.de>
> *Subject:* Re: [DFTB-Plus-User] Mulliken analysis along the trajectory
> [You don't often get email from aradi at uni-bremen.de. Learn why this is
> important at https://aka.ms/LearnAboutSenderIdentification
> <https://aka.ms/LearnAboutSenderIdentification> ]
>
> Dear Donyu,
>
> you can find the atomic gross charges in the xyz output file, just after
> the coordinates, for every atom. This is a convention several programs
> understand.
>
> What do you mean by ESP? The value of the electrostatic potential at the
> position of the atom?
>
> Best regards,
>
> Bálint
>
> On 27.02.24 16:30, Lyu, Dongyu wrote:
>> Dear DFTB-plus-User group,
>>
>> I am trying to perform the ground state dynamics of a molecule in gas
>> phase and run the Mulliken Analysis. In the output files I can find the
>> atomic gross charges in detailed.out , which contains the detailed
>> information for the last frame only. Is it possible to do this analysis
>> along the trajectory and print out the Mulliken charges for each step
>> during the simulation? What should I change in the input file?
>> Also, is it possible to print out the ESP of each atom during the
>> simulation?
>> Here is the input file I used:
>> Geometry = GenFormat {
>>   <<< "TMP.gen"
>> }
>>
>> Driver = VelocityVerlet{
>>   TimeStep [fs] = 1.0
>>   Thermostat = NoseHoover {
>>     Temperature [Kelvin] = 300
>>     CouplingStrength [cm^-1] = 3200
>>   }
>>   Steps = 10000
>>   MovedAtoms = 1:-1
>>   MDRestartFrequency = 1
>> }
>>
>> Hamiltonian = DFTB {
>>   Charge = 4
>>   SCC = Yes
>>   SCCTolerance = 1e-7
>>   MaxAngularMomentum {
>>     C = "p"
>>     N = "p"
>>     H = "s"
>>   }
>>   SlaterKosterFiles = Type2FileNames {
>>     Prefix = "/people/FMO/slko/slko/"
>>     Separator = "-"
>>     Suffix = ".skf"
>>   }
>>   ThirdOrderFull = Yes
>>   HubbardDerivs {
>>     C = -0.1492
>>     N = -0.1535
>>     H = -0.1857
>>   }
>> HCorrection = Damping {
>>     Exponent = 4.0
>>   }
>>
>>   Filling = Fermi {
>>     Temperature [ Kelvin ] = 300.0
>>   }
>> }
>>
>> Analysis = {
>>   MullikenAnalysis = Yes
>>   WriteNetCharges = Yes
>>   CalculateForces = Yes
>> }
>> Options = {
>>   WriteDetailedOut = Yes
>> }
>>
>> ParserOptions = {
>> IgnoreUnprocessedNodes = Yes
>> }
>>
>> I would appreciate it if you could help!
>>
>> Best regards,
>> Donyu
>>
>>
>> _______________________________________________
>> DFTB-Plus-User mailing list
>> DFTB-Plus-User at mailman.zfn.uni-bremen.de
>> https://eur01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fmailman.zfn.uni-bremen.de%2Fcgi-bin%2Fmailman%2Flistinfo%2Fdftb-plus-user&data=05%7C02%7Cdlyu%40constructor.university%7C8931aebf372942f4d4de08dc3878d38f%7Cf78e973e5c0b4ab8bbd79887c95a8ebd%7C0%7C0%7C638447338801663791%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C0%7C%7C%7C&sdata=lO5nfr3HeRb0bx7W%2FryObzlqY35WtoKhiYgMzmkMYTs%3D&reserved=0<https://mailman.zfn.uni-bremen.de/cgi-bin/mailman/listinfo/dftb-plus-user> <https://eur01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fmailman.zfn.uni-bremen.de%2Fcgi-bin%2Fmailman%2Flistinfo%2Fdftb-plus-user&data=05%7C02%7Cdlyu%40constructor.university%7C8931aebf372942f4d4de08dc3878d38f%7Cf78e973e5c0b4ab8bbd79887c95a8ebd%7C0%7C0%7C638447338801677932%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C0%7C%7C%7C&sdata=ol2WPHIWCfrhFf1aEg6f9a1NAWkFXF2uxLtqz7wKM4A%3D&reserved=0<https://mailman.zfn.uni-bremen.de/cgi-bin/mailman/listinfo/dftb-plus-user>>
>
> --
> Dr. Bálint Aradi
> Bremen Center for Computational Materials Science, University of Bremen
> https://eur01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.bccms.uni-bremen.de%2Fcms%2Fpeople%2Fb-aradi%2F&data=05%7C02%7Cdlyu%40constructor.university%7C8931aebf372942f4d4de08dc3878d38f%7Cf78e973e5c0b4ab8bbd79887c95a8ebd%7C0%7C0%7C638447338801686568%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C0%7C%7C%7C&sdata=U%2BHBm%2FnQzPOfkcpK2g9xFNTDXNe6nhmyEsKp1ryQl1Y%3D&reserved=0<http://www.bccms.uni-bremen.de/cms/people/b-aradi/> <https://eur01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.bccms.uni-bremen.de%2Fcms%2Fpeople%2Fb-aradi%2F&data=05%7C02%7Cdlyu%40constructor.university%7C8931aebf372942f4d4de08dc3878d38f%7Cf78e973e5c0b4ab8bbd79887c95a8ebd%7C0%7C0%7C638447338801693797%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C0%7C%7C%7C&sdata=GAzFEr3PxMnxmMpTm5T5QFQftx91b9PO15b26VXvoMw%3D&reserved=0<http://www.bccms.uni-bremen.de/cms/people/b-aradi/>>
>
>
>
> _______________________________________________
> DFTB-Plus-User mailing list
> DFTB-Plus-User at mailman.zfn.uni-bremen.de
> https://eur01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fmailman.zfn.uni-bremen.de%2Fcgi-bin%2Fmailman%2Flistinfo%2Fdftb-plus-user&data=05%7C02%7Cdlyu%40constructor.university%7C8931aebf372942f4d4de08dc3878d38f%7Cf78e973e5c0b4ab8bbd79887c95a8ebd%7C0%7C0%7C638447338801700468%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C0%7C%7C%7C&sdata=YEeeaRL4Rvp%2Bijz3PfmUYJRNY1rBkmh56ubAj5rLjtc%3D&reserved=0<https://mailman.zfn.uni-bremen.de/cgi-bin/mailman/listinfo/dftb-plus-user>

--
Dr. Bálint Aradi
Bremen Center for Computational Materials Science, University of Bremen
https://eur01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.bccms.uni-bremen.de%2Fcms%2Fpeople%2Fb-aradi%2F&data=05%7C02%7Cdlyu%40constructor.university%7C8931aebf372942f4d4de08dc3878d38f%7Cf78e973e5c0b4ab8bbd79887c95a8ebd%7C0%7C0%7C638447338801707305%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C0%7C%7C%7C&sdata=e%2BmzVvJ99lfDVdf9wrHDfbE5BWz%2B34PrNOTZQR76L5k%3D&reserved=0<http://www.bccms.uni-bremen.de/cms/people/b-aradi/>


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.zfn.uni-bremen.de/pipermail/dftb-plus-user/attachments/20240228/0bdd28eb/attachment.htm>


More information about the DFTB-Plus-User mailing list